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Abstract 



Quantum Information Technology (QIT) promises faster, more secure means of data 
manipulation by making use of the quantum properties of matter. In recent years, 
QIT has experienced a steady transition from the realm of the theoretical to that of 
the practical, and seemingly abstract properties, such as entanglement, can now be 
observed and, to an extent, manipulated. However, current research is still investi- 
gating means to control these elusive quantum properties, which seldom survive the 
transition from the microscopic to the macroscopic. This thesis examines some of 
the more fundamental requirements of a successful quantum computation, namely 
the ability to transmit quantum information with maximum efficiency, and the cre- 
ation of entanglement. I focus specifically on neutron entanglement, showing that 
the spins of two or more distinct neutrons can be measurably entangled by forward- 
scattering from an isotropic medium. This result is of great interest, as inter- neutron 
entanglement has never yet been observed. Furthermore, a mathematical treatment 
of neutron scattering puts forth intriguing questions regarding the interpretation of 
'time' in scattering experiments. Though the problem is by no means new, a unique 
resolution is outstanding. I present a simple treatment based on the Heisenberg 
S'-matrix, from which it emerges that in certain situations the quantum-mechanical 
time parameter appearing in the effective time-evolution operator for the spin sys- 
tem has an intuitive physical interpretation. The final part of the thesis deals with 
quantum information transfer in arrays of permanently coupled dipolar systems. It 
is shown that spin chains with dipolar couplings offer high fidelity long-distance 
state transmission, but transfer times in unmodulated chains are unfeasibly long. 
Possible optimization methods are discussed, concluding with a review of recent 
achievements in this field. Finally, I give a brief summary of each chapter, and 
discuss future extensions of this work. 
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Chapter 1 
Introduction 



In the past four decades, the evolution of computer hardware has borne out a predic- 
tion commonly known as Moore's Law. In its most recent formulation [1], Moore's 
Law decrees that the number of transistors on an integrated circuit doubles roughly 
every two years. In a keynote speech delivered earlier this year at the Intel developer 
forum in Shanghai, the senior vice president of Intel technologies stated the trend 
is set to continue for at least another twenty years. By that time, a microprocessor 
no bigger than a few square nanometres will have come to house as many as 10^^ 
transistors [21 E]- Already in the early 1970s, scientists began to realize the implica- 
tions of this continued miniaturization of computer hardware: microchip transistors 
would soon be so densely packed as to experience non-negligible intrinsic quantum 
mechanical effects. The idea of a computer whose functioning would rely in a signifi- 
cant way on quantum processes aroused scientists' curiosity: what if all the internal 
processing of a computer were to follow quantum laws? Would such a machine have 
any advantage over a classical computer? How could it be built? 

One of the first to address these issues was Richard Feynman. During a conference 
at Massachusetts Institute of Technology (MIT) in 1981, Feynman observed that 
a classical computer could not efficiently simulate the evolution of a quantum sys- 
tem, and proposed a basic model of quantum "simulator" that would overcome this 
limitation. Feynman's ideas were formalized by David Deutsch, who in the seminal 
paper "Quantum Theory, the Church- Turing Principle and the Universal Quantum 
Computer" [1] provided a more comprehensive description of a quantum computer's 
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features and capabilities. In comparing a quantum computer to a classical machine, 
Deutscli stated: 

... [the universal quantum computer] admits a further class of programs 
which evolve computational basis states into linear superpositions of each 
other. 

This is, in fact, the defining feature of a quantum computer: the ability to manipu- 
late data that is not in an eigenstate of the computational basis. Whereas classical 
bits can exist only in the states or 1, quantum bits (or qubits), can be found in 
superpositions of the form a|0) + where a and /3 are real or complex coefficients 
satisfying the condition |ap + = 1. Consequently, a quantum computer of N 
registers can exist in a superposition of 2^ possible states, which greatly increases 
its computing power compared to a classical computer of the same size. Consider, 
for example, a function of many variables, f{xi,X2---Xn), with x„ = [0,1]. This 
function is defined constant if it is single- valued for all permutations of {xi, X2...x„}, 
and balanced if its value is for half of the input domain and 1 for the other half. 
Classically, the characterization of the function might require 2^"""^^ + 1 evaluations. 
With a quantum computer, just one evaluation is enough to give an answer that is 
always correct. The quantum algorithm that achieves this is known as the Deutsch- 
Josza algorithm, and was proposed by Deutsch, in collaboration with Richard Josza, 
in 1992 [5]. 

Two years later, Peter Shor of Bell Labs discovered a remarkable algorithm that 
would allow large integers to be factored very quickly, thus invalidating even the 
most secure public- key cryptosystems available at the time [H]. Arriving as it did 
in a period of intense debate on the individual's right to absolute privacy of com- 
munication, Shor's discovery sparked a huge interest in quantum computation and 
kick-started some truly outstanding progress in the field. The very next year, a 
group at NIST (Colorado) succeeded in making the first quantum logic gate with 
trapped ions, according to a proposal advanced previously by Cirac and Zoller of 
Innsbruck University [71[H]. In 1996 Lov Grover of Bell Labs developed an algorithm 
that would allow the quadratic speed up of a database search P , and 1997 brought 
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the first publications on quantum simulation using liquid-state nuclear magnetic 
resonance (NMR) techniques, by Cory, Chuang et al. of MIT [10]. The following 
year, two independent groups at MIT and Oxford University succeeded in building 
the first two-qubit NMR quantum computers, and executed the Deutsch-Josza al- 
gorithm and Grover's algorithm for a database with four entries [TTl [T2| [T3] . The 
new millennium saw the construction of five- and seven-qubit NMR computers, and 
in 2001 researchers at IBM's Almaden Research Centre and Stanford University im- 
plemented Shor's algorithm to factor the number 15 on a 7-qubit NMR computer 

The appearance of a working quantum computer less than fifteen years after its 
'official' conception was a truly impressive achievement, greatly owed to a backup 
of almost five decades of work in the field of magnetic resonance. However, it was 
soon realized that NMR quantum computing would become unfeasible for systems 
of more than 10 qubits, which is considerably fewer than the number required for 
any truly interesting application. As a result, the focus of the scientific community 
largely shifted to other QC strategies, and some of the most popular proposals 
currently include trapped ions, superconducting Josephson junctions, charge states 
of quantum dots, and atoms trapped in optical potentials, usually referred to as 
optical lattices. The past few years have yielded numerous achievements in all these 
fields, but current research is still grappling with issues such as qubit generation, 
manipulation and stability, as well as reliable data storage and transfer. 

This thesis proposes to examine some of the more fundamental requirements of a suc- 
cessful quantum computation, namely the ability to transmit quantum information 
with maximum efficiency, and the creation of entanglement, a stronger-than-classical 
quantum correlation that many believe embodies 'quantumness' itself. I will begin 
with a discussion of this intriguing phenomenon, describing in a mainly qualitative 
fashion its principal properties and its role in quantum computation and commu- 
nication. Experimental achievements will be reviewed, and possible 'gaps in the 
market' will be identified. Chapter [3] attempts to fill one such gap by proposing 
an original scheme to generate entanglement between distinct neutrons, which has 
never been observed to date. The entanglement arises from the sequential scat- 
tering of a neutron pair from a macroscopic substrate, which is described using a 
fully time-dependent approach based on unitary evolution. This somewhat atypical 
treatment of the scattering problem raises the question of the physical interpreta- 
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tion of the quantum-mechanical time parameter, and a detailed comparison with 
standard scattering theory reveals interesting, and often pleasingly simple results. 
Following a discussion of this point in chapter HI the performance of the protocol will 
be assessed in chapter |5l together with experimental requirements and feasibility in 
light of methods and instrumentation currently at our disposal. Chapter [6] then de- 
scribes an extension of the scheme, which allows for greater experimental flexibility. 
In chapter [71 the focus shifts slightly from entanglement generation to entanglement 
distribution by transmission of arbitrary quantum states. The discussion will hinge 
on quantum communication using spin chains, remarkable yet simple systems iden- 
tified some years ago as potentially faithful carriers of quantum information. I will 
review the original proposal for quantum state transfer through nearest-neighbour- 
coupled Heisenberg spin chains, and present a new model for near-perfect spin chain 
communication based on long-range dipolar interactions in chapter [HI Finally, chap- 
ter [9] is dedicated to a short summary and discussion of the main results of this 
thesis, together with a vision for improvements, extensions and future avenues of 
investigation. 



Chapter 2 

The Role of Entanglement 



This chapter introduces quantum entanglement, a property identified by Schrddinger 
in the mid 1930's as the very embodiment of 'quantumness' itself. After a brief 
historical note, I will illustrate the main properties of entangled states and the 
concept of an 'entanglement measure', a theoretical tool used to quantify entan- 
glement and study its dynamics. Entanglement generation schemes will be reviewed, 
with particular emphasis on long-distance entanglement via ancillary mediators. I 
will discuss detection of entangled states using witnesses, and current experimental 
achievements, concluding with a few personal remarks on the ontological value of 
understanding what Einstein called 'spooky action at a distance '. 



2.1 Introduction 



The act of communication between two parties can be interpreted as a movement of 
information from one to the other. This information somehow hnks the two parties, 
hence a successful communication fulfils the objective of establishing a correlation 
between a sender and a receiver. In the quantum world, this correlation has rather 
peculiar properties, and upon close analysis raises questions that stray from the 
realm of physics into that of philosophy. We are speaking, of course, of entanglement. 

From an operational viewpoint, entanglement is a physical resource, which enables 
faithful communication of quantum information and is believed by many to fuel 



12 



Chapter 2. The Role of Entanglement 



13 



the power of quantum computation ( [l5l [16], [T7] and references therein). Entangled 
quantum systems share a stronger-than-classical correlation, such that a measure- 
ment on one will instantaneously project the other into a state conditional on the 
outcome of that measurement, even if they are separated by a space-like interval. 
In the words of Schrodinger himself [T8] : 

When two systems, of which we know the states by their respective 
representatives, enter into temporary physical interaction due to known 
forces between them, and [...] separate again, then they can no longer be 
described [...] by endowing each of them with a representative of its own 
[...] By the interaction the two representatives have become entangled. 

[...] even though we restrict the disentangling measurements to one sys- 
tem, the representative obtained for the other system is by no means 
independent of the particular choice of observations which we select for 
that purpose and which by the way are entirely arbitrary. It is rather 
discomforting that the theory should allow a system to be steered or 
piloted into one or the other type of state at the experimenter's mercy 
in spite of his having no access to it. 

Schrodinger's discomfort stemmed from the apparent conflict between the predic- 
tions of quantum mechanics and the reigning opinion that a complete physical theory 
should be local and realistic [19]. According to such a theory, (i) a measurement 
made at one point in space should not be influenced by anything outside its past 
light cone; (ii) physical properties should have a well-defined value whether they are 
measured or not. It is believed Einstein once said, "1 like to think that the moon is 
there even if 1 am not looking at it" [SU]- The idea that observations made on one 
party of a entangled pair could instantaneously project the other into a well-defined 
quantum state contradicts both these axioms, suggesting (i) that faster than light 
communication is possible, and (ii) that physical quantities lack objective reality. 
Quantum mechanics was therefore deemed incomplete, and local 'hidden variable' 
(LHV) models were proposed to explain the manifestation of this 'spooky action at a 
distance'. After many years of heated debate, British physicist John Bell concluded 
that empirical evidence alone would resolve the controversy, and in 1966 published 
details of an experiment designed to test the deterministic world view [21] . The key 
result of this publication was a measurable inequality, derived independently of any 
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type of mathematical formalism, which should be respected by any LHV. A more 
general version of Bell's inequality was measured for the first time in the 1980s in 
a series of experiments by Aspect et ai, whose results conclusively precluded the 
existence of local hidden variables, therefore ascribing to quantum mechanics the 
characteristics of a non-local theory [22]. The overwhelming majority of Bell's test 
experiments carried out since have reached similar conclusions [23]; however, there 
remain shadows of uncertainty owing to the possibility of 'loopholes', which arise 
from the need to make assumptions about the experimental setup |24)H 




Entangled states have many interesting properties. It has recently been shown that 
entangled particles may signal to each other more than 10000 times faster than the 
speed of light [2S|- Luckily for special relativity, this does not imply super-luminal 
communication, because the signal cannot be used to send a message. For this to 
be possible, control over a measurement outcome is required - but this outcome 
is always probabilistic. Indeed, the information carried by a maximally entangled 
state relates to the whole rather than the components, which are individually in 
maximally mixed states. This property can be phrased in terms of entropy, in that 
the state of the whole is always more ordered than that of the subsystems [26]. 
The operational capabilities of entangled states can then be expressed as entropic 
inequalities which show, for example, that if two parties share a maximally entangled 
state and a classical channel, each party can arrive at full knowledge of the state 
without any further exchange of qubits [2Z|. 

Entanglement is perhaps best appreciated in the context of quantum computation. 
It is believed that entanglement is, at least partly, responsible for the exponential 
speedup of a quantum process over a classical one. It has been shown that an ar- 
bitrary quantum computation requires exponentially higher resource overheads if 
entanglement is not involved [T^l [IHl • This view is upheld by the fact that many 
quantum algorithms, such as Shor's factorization and Grover's search, rely on the 
production of highly entangled intermediate states [28l[29]. Sophisticated quantum 
search (SQS) algorithms capable of outperforming a classical algorithm without re- 
course to entanglement have been suggested; however, it has also been observed 
that equally efficient classical versions of SQS exist, therefore no true advantage 

^An example is the 'fair sampling' assumption, which relates to optical realizations of the 
experiment, and states that the number of photons reaching a detector is representative of the 
number of photons emitted by the source. 
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over classical computation can be claimed [301 EI] • Simulations of quantum systems 
also rely on entanglement, because not all quantum processes can be reproduced by 
a classical computer. In some exceptional cases, highly entangled quantum compu- 
tations can be efficiently simulated by classical means [321 [331 IM]- However, as a 
general rule, this is not possible unless the amount of entanglement present in the 
system does not exceed k logn, where k > and n is the number of input qubits 

[13 Eg. 

Entanglement is by no means a comprehensively understood phenomenon, however 
despite its more esoteric characteristics it is now well established as a real phys- 
ical resource. Unfortunately, entanglement is often somewhat fragile, and can be 
diminished or destroyed if the entangled systems experience contact with a noisy 
environment. Much of today's entanglement theory comes from attempts to preserve 
the precious correlations, by determining how entanglement can be created, manip- 
ulated and quantified. Bi-partite entanglement is well understood and, now, acces- 
sible from an experimental viewpoint also. The field of multi-partite entanglement 
characterization is currently one of ever expanding horizons, strongly motivated by 
recent observations of three- to eight-party entangled states [36] . 



2.2 Fundamental Properties 

The main property of an entangled state is that its component subsystems are no 
longer independent. Formally, this means an entangled n-party state p cannot be 
written as a mixture of separable states, i.e. 

pA,B...N Y^^^p^^ (2.1) 
i 

with 

P^ = pt®pf<^■■■p^. (2.2) 

It follows that all tensor product (or separable) states contain no entanglement, while 
all non-separable states are entangled, and therefore allow some tasks to be achieved 
better than by local operations and classical communication (LOCC) alone [HZlESj- 
Separable states can evolve into entangled states, however this requires non-local 
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operations such as two-qubit gates or interaction with ancillary systems. Indeed, 
the total entanglement of an n-party state cannot increase under LOCC alone [39] . 
This goes back to entropy, as there exist theorems stating that a pure state ip can 
only be transformed into another (p if the subsystems of tp are less ordered than those 
of [381 Ho]. Therefore, as the entanglement of a state is directly proportional to the 
entropy of its subsystems, a locally ordered, separable state cannot be transformed 
into a locally disordered entangled state. LOCC can however change the distribution 
of entanglement within a many-particle system; this property enables protocols such 
as quantum teleportation and quantum dense coding, which allow Alice to transmit 
a message to Bob by performing local measurements on her half of a previously 
shared entangled state [HI |l2]. It has also been shown there exists a special class 
of entangled states that increase the scope of LOCC if employed as ancillae. These 
are called bound entangled states, sometimes referred to as 'activators', owing to 
their ability to facilitate entanglement-boosting protocols such as distillation. An 
interesting result is that any bipartite pure state can be transformed into another 
by LOCC alone, provided one has access to bound-entangled states [HZl SSI Sll US] • 

The invariance of entanglement under LOCC can be related to a property known 
as entanglement monogamy According to this property, there is a limit on the 
amount of information a particle can share with others. If a particle is maximally 
entangled with another, its entanglement with any other system must be zero. This 
is owed to the fact that a maximally entangled state contains all the information 
available on that state; there is none left to share with other parties. If, on the other 
hand, two particles are only partly entangled, correlations with other parties are 
allowed. The bound on the amount of entanglement between any pair of particles of 
a multiparty state has recently been formalized in terms of an entanglement measure 
known as the tangle, which will be introduced in the following section. 



2.3 Measures 

Analyzing the entanglement of a quantum state can mean several things. On the one 
hand, one could be content with establishing whether the state is in fact entangled or 
separable. For bipartite systems whose dimensions d do not exceed 6, this is easily 
decided via a criterion known as the Positive Partial Transpose (PPT) |m 1^ . 
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The criterion states that if the eigenvalue spectrum of a state pab becomes negative 
under partial transposition with respect to either system, the state is entangled. This 
does not mean a state with positive partial transpose is necessarily separable; the 
aforementioned activator states, for example, are entangled states with PPT. The 
situation becomes more complicated for multipartite states; in this case, the PPT 
criterion can tell us whether any two parties are entangled, but not whether genuine 
multipartite entanglement exists. It is, therefore, often simpler to quantify the 
entanglement of a state via ad hoc functions known as entanglement measures. Not 
all measures have the same operational meaning, therefore different entanglement 
measures can afford different state ordering. Calculating a spectrum of entanglement 
measures for a given state can then provide a useful means of characterization. 

Entanglement measures E satisfy the following fundamental axioms: 

1. E{p) = if p is separable; 

2. E{p) does not increase over LOCC, and is therefore an entanglement mono- 
tone. 

As regards pure states, the best known and most commonly used entanglement 
measure is the Von Neumann Entropy, defined as 



This is the quantum equivalent of the Shannon entropy, which measures the un- 
certainty associated with a classical probability distribution [19]. In the context of 
pure-state entanglement, S (p) represents something of a 'reference', as all entangle- 
ment measures E{p) must reduce to S{p) if p is pure. 

In more general terms, any function that adheres to the two conditions listed above is 
a valid entanglement measure. Therefore, the set of measures is in principle infinite 
|50j . However, for reasons of mathematical convenience, it is often useful to choose 
functions which satisfy more specific conditions [38], for example: 

1. Normalization: E{p) = 1 if p is a Bell state, or related to a Bell state by a 
local unitary transformation; 



S (P) = -Tr^ {B) [pA (B) logs PA {B)] ■ 
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2. Convexity: For a mixed state p = YliPiPi^ where the Pi are the weights of the 
different pi components, an entanglement measure is convex if 



Some entanglement measures have an operational definition, others are constructed 
on the basis of the above axioms, and not all have a physical interpretation. The brief 
summary that follows is intended as a largely qualitative overview rather than an 
exhaustive and mathematically detailed discussion, as I have focussed on measures 
that are either relevant to the present work or have an intuitive physical interpreta- 
tion. I refer the reader to [3H1 EI] and references therein for a more comprehensive 
and fully mathematical treatment. 

The Distillable Entanglement Entanglement distillation (or concentration) is 
an entanglement purification protocol whereby Alice and Bob share n copies 
of an entangled state p, perform LOCC and obtain k < n copies of a Bell 
state [22] ■ The distillable entanglement is defined as the optimal ratio k/n 
yielded by this process in the limit of large n. It is interesting to note that all 
entangled bipartite states can be distilled, however this may require the help 
of an activator state. Activator states, on the other hand, cannot be distilled; 
their entanglement is inaccessible, hence the term 'bound entanglement'. 

The Entanglement Cost represents the number of singlet states n one has to 
share to distill an arbitrary state, such that errors become infinitesimal in the 
limit of large n [52]. It is therefore defined in relation to a process which is 
the opposite of distillation. 

The Entanglement of Formation was introduced in the context of entangle- 
ment purification, with the principal objective of quantifying the entanglement 
of mixed states, i.e. states which are a mixture of entangled pure states [39] . 
The entanglement of formation is then defined as the minimal average entan- 
glement over all pure state decompositions. In a two-qubit system, this can 
be calculated exactly from a second quantity, the concurrence [53], defined as 




(2.4) 



C{p) = max{0, Ai — A2 — A3 — A4}. 



(2.5) 
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Taking p to be the density matrix of the two qubits, the Aj are the eigen- 
values in decreasing order of the matrix R = a/ \fpP\fPi which has p = 
{o-y<S>(^y)p*{(^v<S>(^v) and 




ay 



where p* is the element-by-element complex conjugate of p. 

The Logarithmic Negativity A bipartite state is entangled if it has a negative 
partial transpose. The magnitude of the negative eigenvalue can be related to 
the degree of entanglement using the logarithmic negativity [5^ [55] . defined 
as 

E;v(p) =log2||p^^|| (2.6) 



where ||^|| = TrvJOX, and Tb indicates the partial transpose with respect 
to qubit B. 

The Robustness of Entanglement This quantifies how much noise it is possible 
to mix in with a given entangled state before it becomes separable 



The Squashed Entanglement is related to the intrinsic information, which mea- 
sures the correlations between two discrete random variables whose value is 
conditioned by a third ^7\. It is proportional to the minimal quantum intrinsic 
information over all states Pabc satisfying Tr^ Pabc = Pab- 

E,, = inf{,,,,|i [S{AC) + SiBC) - S{ABC) - 3(0)] , (2.7) 

where S {■) is the Von- Neumann entropy. The squashed entanglement is re- 
lated to the distillable entanglement and the entanglement of formation, as 
Ed < Esq < Ep. 



Strictly, the above quantities relate to bipartite entanglement only, though multi- 
partite generalizations can sometimes be derived. An example is the 'global entan- 
glement' [58], which is simply the sum of the pairwise concurrences of a many qubit 
state. The squashed entanglement can also be extended to many particle states 
by replacing the bipartite mutual information with its equivalent for multipartite 
systems [59] . 
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A commonly used measure is the tangle t [16], defined for tliree qubits as 

t{A, B,C) = C'^^^Q — C^^^ — C'^ Q, (2.8) 

wfiere C is tlie concurrence. It lias been shown that C]^ bc ~ PAi where pa 
is the reduced density matrix of the first subsystem. This is owed to the fact that 
the pair BC can be viewed as a 'unit'; therefore, the four- dimensional subspace 
of BC effectively becomes two-dimensional. The tangle is then the entanglement 
between system A and the unit BC which is not accounted for by the separate 
entanglements of A with B and C individually. This is sometimes referred to as the 
'residual entanglement'. In its original context, the tangle was used to define the 
monogamy inequality for three qubits 

C^, — ^1, B + ^A, c- (2-9) 
Recently, this has been generalized to m-qubit systems [60]. One then has 

Ci% + --- + Ci2,™<Cj(2,..^), (2.10) 

hence the entanglement of one qubit with all the others individually cannot exceed 
the entanglement of that qubit with the others taken as a block. The quantity 
^1 (2 -m) defined in analogy to via the so called convex roof construction 

[6T] . Taking qubits {2,m} to represent a 2'^™~^^-dimensional system B, one finds 

^Ib = ^T^kp^^^^y^p^si^'^B p^J, (2.11) 

X 

where the infimum runs over all pure state decompositions of pi^ = YlxP^Pip^ 
yielded by the pure state ifjx- 

Another interesting measure, showing some formal similarities to the bipartite en- 
tanglement of formation, has been defined in the context of spin chains. This is the 
localisable entanglement [^2], and represents the maximum bipartite entanglement 
one can obtain by choosing two subsystems of a multipartite entangled state and 
performing LOCC to optimize their correlation. For three qubit states, the localis- 
able entanglement is equivalent to the so-called entanglement of assistance, i.e. the 
amount of entanglement Alice and Bob can establish between their qubits if Charlie 
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measures his and communicates the outcome [63] . 



For reasons of relevance and brevity I shall stop here, although the above is by no 
means an exhaustive list. For a comprehensive discussion, I once more refer the 
reader to [38l [5T] and references therein. 



Entanglement measures provide a theoretical characterization of entanglement dy- 
namics, but cannot be measured directly because they do not represent observables. 
To verify the entanglement of a system in an experimental setting, it is therefore 
useful to define a set of hermitian operators termed entanglement witnesses [IHlEl], 
whose expectation values represent average values of measurable quantities, i.e. ob- 
servable properties of a physical systems. A valid entanglement witness must satisfy 



if pab is separable, hence the expectation value of a witness must be positive on all 
product states. A witness therefore provides an experimentally accessible means of 
establishing whether a state is entangled or not. 

The estimation of the entanglement properties of a system from expectation value 
measurements bears an interesting relation to the Jaynes maximum entropy principle 
[65] . It is known that the state of a system is completely characterized by the 
ensemble of the expectation values of all commuting observables relating to that 
state. Suppose, however, that one disposed of an incomplete data set which did 
not allow a unique characterization; the Jaynes principle affirms that, of all possible 
states of the system, the most probable one is that which maximizes entropy. Let us 
now imagine we hold an incomplete data set relating to some quantum state; what 
could we infer about its entanglement properties? According to the Jaynes principle, 
if the data set is compatible with a range of entanglements e G [e_,e_|_], one should 
always assume e = e_|_. This is due to the correspondence between entanglement and 
entropy discussed above. However, it has been shown that such a maximization is not 
always 'truthful'; in other words, by applying the Jaynes principle to quantum states 
one runs the risk of overestimating the amount of entanglement at one's disposal 



2.4 Witnesses 



{W) = Tr [Wpab] > 



(2.12) 
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[66] . It has therefore been suggested that in a quantum setting one should apply 
minimization of entanglement under experimental constraints. In other words, if 
an incomplete data set is compatible with the range of entanglements e G [e_ , e+j , 
one should always assume the true entanglement of the system is represented by the 
lower bound. 

Witness operators can be ordered by the entanglement they detect. A witness is 
finer than another if it detects more entanglement, and is termed optimal if no finer 
witness exists, i.e. if it detects all the entanglement of a given state. Witnesses are 
usually represented by a combination of observables, rather than a single one. As the 
measurement of each component of the witness typically requires a different setting 
of the laboratory apparatus, much effort is dedicated to finding optimal witnesses 
that can be measured with a minimal number of device settings. This problem is 
sometimes referred to as 'finding the optimal decomposition' [671168]. Constructing 
a witness capable of detecting the entanglement of an arbitrary state is generally 
an arduous task. However, the problem simplifies greatly if one has some a priori 
knowledge of the structure of the target state. In this case, it is possible to define 
witnesses that are measurable with as few as three device settings |67j . 

In practise, the spectrum of measurable witnesses for a given system is restricted by 
experimental constraints. However, witnesses find application as theoretical tools 
also, because for any valid witness W, condition fl2.12|] provides a separability cri- 
terion. Conditions such as f l2.12p are known as linear separability criteria, because 
only the first power of the expectation value of the witness appears. Non-linear 
criteria feature higher powers of (W). Typical examples involve the uncertainties 
associated with measurements of local variables, expressed in terms of the variance, 
which is a quadratic function of the expectation value [69l 170] . 

2.5 Entanglement Generation 

A vast proportion of all research in quantum computation and communication relates 
to entanglement generation and control, with a view to find experimentally feasible 
ways of producing robust entanglement on demand. Nearby qubits are best entan- 
gled via two-qubit gates; however, gating operations become inefficient on distant 
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qubits, whose direct coupling is usually extremely weak. This problem represents 
one of the main obstacles to upscaling current devices, and has therefore been looked 
into by several authors. Many existing proposals for creating long distance entangle- 
ment rely on single- or two-photon interference effects [7T1 [721 [73l [TH [751 [76] . In these 
types of schemes, distant qubits are entangled via the detection of photons result- 
ing from their radiative decajQ- Recently, however, the focus seems to have shifted 
somewhat to the alternative scenario of employing entanglement mediators. These 
are ancillary qubits, made to interact successively with two non-communicating 
qubits A and B; the entanglement between A and B arises when measurements on 
these ancillae are made. Various ways to implement this idea have been proposed 
[771 [751 [791 [80] . In [77], for example, the authors describe an experiment to en- 
tangle two optical cavities by repeated interaction with a mediator. The success 
of the scheme is shown to depend on the interaction times and strengths, and on 
the measurement of the mediator. Yuasa and Nakazato [75] build on this scheme 
to demonstrate that maximally entangled states can be created by three successive 
interactions of identical mediators with the cavities, provided one assumes rotating 
wave couplings of a specific strength, acting for specific times. In later work, the 
same authors find that repeated resonant scattering of a mediator also generates 
maximal entanglement, particularly if the momentum of the mediator can be fixed 
to a desired value [7^. Costa et al. have analyzed the more complex scenario of 
an injected electron being spatially scattered by two spin-| magnetic impurities in 
a one-dimensional, non-magnetic, metallic chain. In this case, it was found that the 
electron could almost perfectly entangle the impurities, conditional to a favourable 
(and reasonably probable) outcome of a measurement on the electron's spin. 

Very recently, some authors have been investigating the issue of entanglement on a 
macroscopic scale [5T 1 [521 [53] . with a view to underline that: 

... entanglement is not an artificial mathematical property but [..] can be 
extracted and therefore used for quantum information processing in the 
same way as heat can be extracted and used for work in thermodynamics 



The situation described in [HI] is that of two probe particles simultaneously scatter- 



^The entanglement arises because it is impossible to trace the provenance of the individual 
photons, hence the qubits are left in a superposition of ground and excited state. 
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ing from two entangled spins in a solid. It is shown that the amount of entanglement 
extractable from the solid depends mainly on the interaction Hamiltonian and the 
joint state of the spins. 

Another interesting proposal is that described by Cunha and Vedral [H2] • This builds 
on previous work by the same authors, in which the Pauli principle is invoked as 
a means to create entanglement between fermions in a degenerate Fermi gas. Be- 
cause fermions only have spatial and spin degrees of freedom, it is argued that if 
two fermions were simultaneously detected at the same position, their spin degrees 
of freedom would necessarily be entangled. The authors show this entanglement 
persists provided the fermions remain within a certain distance of each other. The 
entanglement distance is found to be inversely proportional to the Fermi wavevector 
kp, which imposes constraints on the type of fermion one could use for an experi- 
mental verification of this scheme. The setup described involves a vessel containing 
ultracold neutrons, with a small hole. The authors suggest that provided certain 
precautions are observed in preparing the Fermi gas, the neutrons emitted from this 
hole within a given time interval of each other will be entangled. The scheme is 
simple and efficient, however an experimental realization would be extremely chal- 
lenging, owing to the difficulty of preparing a Fermi gas of neutrons in a laboratory 
environment. 

The Pauli principle is used as a means to create entanglement also in [S3]- Here, 
the authors consider two electrons with opposite momenta and spin, trapped in a 
potential well inside an optical cavity. Two independent photons with orthogonal 
polarizations are fired into the cavity, creating two excited electrons with different 
spin and momenta. The electrons are made to interact with some external energy 
absorbing reservoir, which "drags" them into the same zero-momentum state, caus- 
ing their spins to become entangled. The entangled electrons then decay according 
to well-defined selection rules, emitting spin-polarized entangled photons. Interest- 
ingly, the dissipative interaction with the environment is essential to the success of 
the scheme, because the photons become entangled only if the electrons' momenta 
are equal. This is ensured by the dissipative interaction with the environment, hence 
decoherence is quite surprisingly viewed as a 'pro' rather than a 'con'. 
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2.6 Experimental Achievements 

The study of entanglement was in origin a largely theoretical enterprise. The ad- 
vancement of experimental prowess in the past few years now means the scientific 
community is poised to mount a two-pronged attack on the subject: theory tells us 
what we should see, experiment tells us what we do see. Attempting to reconcile 
the two produces a positive feedback, whence an understanding of the relationship 
between the microscopic and macroscopic world is born. 

Recent experimental studies of entanglement generation have pursued two main 
objectives. On the one hand, great effort has been dedicated to the realization of 
many-particle entangled states. These are essential to measurement-based quantum 
computation with graph or cluster states [HI], which provides an alternative to the 
more widespread gate model. Achievements in this area include the realization of 
four- to eight-ion W-states [36], four- ion [85], five-photon [HS] and three-atom [H7j 
GHZ states . 

Much research is also directed at sourcing novel systems to entangle. Indeed, entan- 
gled particles are not only at the heart of quantum computation and communica- 
tion, but also find application in more fundamental areas, such as investigating the 
non-locality and non-contextuality of quantum mechanics. Photons, ions, atoms, 
quantum dots [HE] and superconducting qubits [HSj are all demonstrably feasible 
possibilities, but neutrons, for example, have never been entangled [101 [HI - an is- 
sue this thesis will attempt to address. An area of rapidly growing interest is that of 
single-particle entanglement [02], which involves entangling distinct degrees of free- 
dom of the same particle. Systems of this type clearly cannot exhibit non-locality, 
though non-contextual behaviour has been observed for both single neutrons and 
single photons [90l [93]. One of the most impressive achievements to date involves 
the creation of a three-way entangled GHZ-tjpe state of a single neutron, the three 
'qubits' being the spatial and spin parts of the neutron's wavefunction, and its energy 
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2.7 Conclusions 



The study of entanglement occupies a prime position in the field of quantum in- 
formation processing. From a practical viewpoint, it is clear that generating and 
preserving entanglement in a controllable fashion are paramount to building quan- 
tum computers that rival existing classical machines. In addition, entanglement is 
a probe of nature, given its potential to reveal whether the world around us has 
objective reality or whether, to paraphrase slightly, reality is in the eye of the be- 
holder. Finally, entanglement is a probe of ' quant umness', because it is destroyed by 
contact with the macroscopic. I believe the practical and fundamental aspects are 
closely interwoven: to use a resource to its full potential, one must first understand 
its properties. However, if cracking an RSA-encrypted message and understanding 
'life, the universe and everything' truly are two sides of the same coin, I'd place 
strong odds on achieving the former before the latter. 



Chapter 



3 



A Novel 



Proposal for Neutron 



Entanglement 



In this chapter, I put forward a proposal for a scattering experiment aimed at achiev- 
ing entanglement between two distinct neutrons. This has never been accomplished to 
date; indeed, neutron entanglement has only been observed between distinct degrees 
of freedom of single particles. The realization of entangled states of two neutrons is 
therefore a highly desirable objective, and would open up the possibility of testing the 
non-local nature of quantum mechanics with an as yet untapped resource. After a 
general overview of neutron scattering as an experimental technique, I will describe 
the protocol, and discuss the assumptions and approximations underlying it. In the 
concluding part of the chapter, I outline the mathematical formalism employed to 
model the experiment, and comment on its relationship to the traditional approach 
to scattering problems, which will be treated in detail in the next chapter. 

3.1 Neutrons and Scattering 

§ A great deal of our current understanding of the structure and dynamical proper- 
ties of matter comes from neutron scattering. This technique, originally developed 
to further the study of condensed matter physics and crystallography, now finds 

-'^The majority of this material was obtained from [94] , 
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application in the most diverse areas of science, including biology, nanotechnology 
and planetary science. 

Neutrons are cliargeless particles which, together with protons, form atomic nuclei. 
The standard model of particle physics classifies them as baryons, the subclass of 
hadrons composed of three quarks - two down and one up in this specific case. Free 
neutrons are unstable, and have a typical lifetime of the order of 15 minutes, after 
which they decay into a proton, an electron and a neutrino. Scattering events are 
usually concluded well within this time frame. Owing to the absence of charge, 
neutrons are largely insensitive to stray electric fields and surface defects, and do 
not themselves create fields capable of disturbing the environment of the scatterer. 
Consequently, neutrons are excellent non-destructive probes both of surface and 
bulk properties of matter, and can be used to study samples immersed in complex 
environments |95j . 

Typical scattering experiments are aimed at the characterization of a sample. For 
reasons that go back to wave-particle duality, this is achieved by analyzing the mo- 
mentum distribution of the outgoing neutrons. Indeed, the neutron scatters as a 
wave but is detected as a particle. The De Broglie wavelength of a neutron is de- 
termined by its energy; so-called thermal neutrons have wavelengths comparable to 
interatomic spacings, and can therefore penetrate bulk materials. The energy of a 
thermal neutron is of order O (10~^) eV, which is similar to the energy scales of 
elementary acoustic or magnetic excitations (phonons and magnons). Momentum 
changes of scattered neutrons can then be used to map dispersion curves for lattice 
vibrations and spin fluctuations. Elastic scattering, on the other hand, provides 
information on the structural properties of the sample, and, on a more fundamental 
level, on the optical potential that couples nucleons to nuclei and to each other. 
Finally, if the resolution of the experiment is sufficient to distinguish elastic and in- 
elastic processes, scattering can also identify whether or not the sample is entangled 

m- 
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Figure 3.1: A visual description of neutron scattering as an experimental technique. Image repro- 
duced with kind permission of Professor Sunil K. Sinha, University of California San Diego and 
Los Alamos National Laboratory. 



The interaction between the neutron and the sample is determined mainly by two 
mechanisms: the strong interaction, which couples the neutrons to the nuclei, and 
the magnetic dipole interaction, which couples the neutrons to the sample's magnetic 
moment [97] . The latter mechanism is owed to the neutron's non-zero magnetic 
moment, given by 

s 

fj^ = gNiJ'N^, (3.1) 

where Qn = —1.913 is the neutron gyromagnetic ratio, s is the neutron spin operator 
and /iAT is the nuclear magneton, whose magnitude is approximately 1000 times 
smaller than the Bohr magneton. The potential experienced by a neutron during 
a scattering event is then a combination of a short-range interaction, arising from 
the nuclei, and a long-range potential, arising from the sample's unpaired electrons. 
Both interactions are typically rather weak, hence it is common in mathematical 
models of experimental data to treat scattering as a first order process [HH] ■ 
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Neutron scattering is unfortunately a costly technique, which must be carried out in 
apposite facilities. These come under two categories: fission research reactors, and 
spallation sources. Fission reactors, also known as continuous sources, produce a 
steady stream of neutrons via a controlled chain reaction. When heavy elements such 
as uranium or plutonium absorb a slow neutron, the nucleus splits into two lighter 
components, releasing two or three fast neutrons in the process. These are referred 
to as prompt neutrons, because they are the first products of the reaction. They 
are subsequently followed by delayed neutrons, produced by the decay of the fission 
fragments. Each of these neutrons can trigger a new fission event, giving rise to a 
chain reaction. This is kept in check by removing excess neutrons, which are slowed 
by passage through a moderator and absorbed by control rods. The moderator is a 
large volume of material such as water or solid graphite, that absorbs weakly but 
scatters strongly. The rods, on the other hand, are typically boron-rich, as boron 
has a high neutron capture cross section. 

It is sometimes desirable to produce pulses of neutrons rather than a continuous 
stream. This offers a means to distinguish the processes contributing to the scat- 
tering event, which do not typically occur on the same timescalej^. Fission reactors 
can be made to operate as pulsed sources with the help of a chopper. Essentially, 
this is a shutter positioned downstream from the source, which is opened and closed 
periodically so that only neutrons with well-defined velocities pass through. Pulsed 
beams can also be produced by spallation, which involves bombarding a heavy nu- 
cleus with a high-energy proton beam, usually produced in a particle accelerator. As 
a result of the collision, several neutrons are knocked out of the atom, while others 
boil off as the nucleus heats up. Spallation typically produces 20 or 30 neutrons per 
collision, but requires more energy than it delivers, therefore chain reactions cannot 
occur. 

Once the neutrons have been produced, they are typically directed away from the 
source by a neutron guide. This is a channel with high reflectivity at small incidence 
angles, often curved to reduce the line of sight from the source to the sample, hence 
avoiding contamination of the scattered signal by stray high-energy neutrons. Upon 
emerging from the guide, the neutrons often undergo some kind of velocity selection, 
after which their fate is determined by the purpose of the experiment. 

^Inelastic and fast-neutron processes occur shortly after the pulse, whereas neutron capture 
processes follow some time after. 
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Finally, once scattering has taken place, the outgoing neutrons are detected. This is 
not a trivial task, because most currently available detection techniques are designed 
for charged particles. High-efficiency detection can be achieved using hydrogen-rich 
scintillators, because hydrogen is a strong scatterer. Otherwise, detection occurs via 
a two-step process, typically involving a reaction that produces charged by-products 

PIES]. 

The field of neutron scattering is incredibly vast; the above is intended to provide 
only the most basic overview. Where necessary, I will go into more detail with 
regard to instrumentation and experimental techniques. The formal description of 
neutron scattering will be the subject of a separate discussion. 

3.2 Notation 

Having reached this point, it is convenient to define some 'shorthand' notation, 
which will be employed throughout this thesis. Unless otherwise specified, quantum 
states will be expressed in accordance with the Bloch sphere convention, which la- 
bels the spin-up and down eigenstates of cr^ with the kets |0) and |1), respectively. 
The wavefunction of a many-body system fully polarized in the positive (negative) 
z direction will be represented as |0) (|1)). Depending on the magnitude and orien- 
tation of an external magnetic field along z, either |0) or |1) may correspond to the 
system's ground state. Excitations from the ground state are obtained by rotating 
individual spins through an angle vr. A particle whose spin is anti- aligned with the 
ground state configuration will be referred to as a 'spin fiip'. A fully aligned state 
with a single spin fiip at site j will be represented as a fully aligned state with 
spin fiips at sites j and / will be represented as \jt) , and so on. 

For reasons of mathematical simplicity, the majority of the analysis carried out in 
this thesis is concerned with pure states. Often, these will be described in terms of 
a density operator, conventionally labeled p. For a state vector this takes the 
form 

p=|^)(^|. (3.2) 

If p represents a many-particle system, the density operator of a subsystem is ob- 
tained by performing a partial trace. For a system of two qubits A and B, one then 
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has 

Pa (b) = Tr^ (s) [pab] , (3.3) 
where pab represents the composite density matrix. 

Variables and operators will usually be assigned a set of qualifying indices; how- 
ever, these will be omitted whenever the context allows, to make for more compact 
notation. 

3.3 Approximations 

To provide a more complete picture of the scope and interpretation of this work, it 
is useful to head the description of the protocol with an explanation of the approx- 
imations that underlie it. These have been chosen according to the twin criteria of 
analytical computability and experimental feasibility, and much care has been taken 
to strike a balance between the two, as the ultimate goal is to design a feasible ex- 
periment. 

As regards the operating conditions, it will be assumed the sample is at zero tem- 
perature and completely isolated from surrounding noise. This removes the risk of 
thermally induced spin flips or other forms of decoherence. The sample and the 
neutrons are immersed in a static, uniform magnetic field aligned with the z direc- 
tion. This is chosen as the quantization axis. The effect of the field on the neutrons 
will be neglected, as the neutron magnetic moment is O{10^) times smaller than the 
electron magnetic moment. The neutrons themselves are produced by an ultra-cold 
source capable of imparting an arbitrary, user-defined polarization. Neutron veloci- 
ties are chosen not to exceed 10 ms~^, with corresponding energies of order 10~^ eV. 
Sources of this type have been in operation for some years now, an example being the 
PF2 source at the Institute Laue-Langevin in Grenoble (see Appendix [B] and |100] ). 
It is assumed the intensity of the source is low, so that very few neutrons interact 
with the sample at any chosen moment. This should not be a problem, as scattering 
facilities go to great lengths to avoid producing dilute beams. High-intensity sources 
could also be adapted to serve our purpose, provided they were operated in pulsed 
mode - indeed, when a continuous beam is chopped, a large proportion of the flux 
is lost. In addition, neutrons are weakly interacting particles with an extremely low 
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scattering cross section; therefore, the probabihty of many neutrons scattering at 
once is small. 

The sample from which the neutrons scatter is modeled as a magnetically ordered 
crystal, in which the unpaired electrons are localized near the equilibrium positions 
of the ions in the lattice. This is known as the Heitler-London model, and provides a 
good description of ferromagnetic insulators such as EuO, LaMn03.i2 and few other 
materials |101l 11021 1103] . Magnetic scattering is provided by the long range dipolar 
force between the neutron and the unpaired electrons, while nuclear scattering is 
provided by the short-range strong force. I assume LS coupling, i.e. that orbital 
and spin angular momenta of the individual ions can be summed to yield total 
orbital and spin angular momenta, defined by the quantum numbers L and S. I 
work in the limit of L = 0; this could describe a lattice of s-shell electrons, and also 
a more general situation in which the resultant orbital angular momentum is either 
zero, or quenched by the internal electric field of the crystal. Magnetic scattering 
of the neutrons will be teated as a purely electronic effect; the nuclear spins will 
be completely neglected. This is reasonable, as the strength of magnetic scattering 
from nuclear dipole moments is approximately two orders of magnitude weaker than 
that from electronic dipole moments |104] . 

The electrons in the sample are coupled by a translationally invariant potential 
which conserves the total spin quantum number Sz- A reasonable example of such 
a potential is Coulomb nearest-neighbour exchange, given by the Heisenberg Hamil- 
tonian 



where the notation (ij) indicates the sum over all nearest-neighbouring pairs, and 
Jij represents the coupling strength between nearest neighbours. For simplicity, I 
will take this to be constant and positive, so as to describe a ferromagnetic sample. 
The full Hamiltonian of the system then reads 



where 5^ > is the strength of the external field and is the number of spins in 




(3.4) 



N 




(3.5) 



the sample. Each spin is represented by the operator a 



H^^y^o-^). where the 
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0"^ are the Pauli matrices for spin j. 

The sample itself is a macroscopic, periodic system experiencing negligible boundary 
effects. Each neutron scatters from the sample as a wavepacket of spatial width 
6x and momentum width 6p, whose relationship is determined by the uncertainty 
principle 

SJp > |. (3.6) 

It will be assumed 6^^ is comparable to the size of the sample, or, equivalently, that 
the neutron interacts with all the sample spins at once (figure I3l2l) . The validity of 
this assumption will be discussed in section 15.51 

k 




Figure 3.2: A one-dimensional snapshot of the scattering geometry. The sample sees the neutron 
as a wavepacket, the peak of which travels with momentum kg. 



The scattering events are described as sequential, time-dependent interactions with 
the potential occurring such that no direct interaction between the neutrons takes 
place at any time. Only neutrons scattered into a specific Bragg peak will be de- 
tected. Upon emerging from the source, the neutrons are uncorrelated and polarized 
in the same (arbitrary) direction, specified by the coefficients a and which can 
take on any real or complex value satisfying |ap + = 1. To begin with, I will 
consider a simplified model in which only two neutrons scatter from the sample. The 
initial state of the system as a whole is then a product of the states of the neutrons 
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and the sample. It is convenient to represent this as 

l^o) = («|0) + (a|0) + \r), (3.7) 

where |V'"')m = (q^|0) + is the state of neutron m, and s stands for 'sample'. 

Hereafter, I will refer to the system of neutrons plus sample simply as 'the system'. 
With these details in place, I now describe the entanglement scheme. 

3.4 The Protocol 

The protocol can be thought of in terms of four phases: 

1. The first period of free evolution of the sample, before the arrival of the first 
neutron pulse; 

2. The interaction between the first neutron and the sample; 

3. The second period of free evolution of the sample, before the arrival of the 
second neutron pulse; 

4. The interaction between the second neutron and the sample. 

From a qualitative viewpoint, the first scattering event has the effect of depositing 
the 'spin information' carried by the first neutron in the sample. As we will see in 
chapter El this signature can remain intact during the second period of free evolution. 
The sample therefore behaves as an 'entanglement safety-deposit box', which retains 
the information deposited by the first neutron until the second neutron comes to 
claim it. However, the entanglement between the neutrons is conditional on the 
second scattering event being able to 'undo' -in part or in full- the transformation 
on the sample spin precipitated by the the arrival of the first neutron (see section 

EM- 

The protocol begins by initializing the sample and the neutrons to a known state, 
chosen such that \iPq) contains at most two spin flips. The motivation for this re- 
striction is mathematical rather than physical, and stems from the hope of attaining 
an analytically solvable model by reducing the size of the working Hilbert space (see 
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section H75|l . Two possible combinations of |^/;") and {tp'^) are considered here, and I 
label these \iPq) and I^Pq). In state I'iPq), the neutrons are polarized parallel to the 
positive i:-axis and the sample is prepared in a single spin-flip state. In the 
neutrons have arbitrary polarization while the sample is prepared in a fully aligned 
state. Hence 



1 ^ 



lO = (a|0)+/5|l)), ia\0) + m)i |0). 
From equations fl3.7p and fl3.8l) . for \iPq) we have 



(3.8) 
(3.9) 



a 

\r)i 

W) 



1, 

0, 

r)2 = io), 
1 ^ 



(3.10) 
(3.11) 
(3.12) 

(3.13) 



The state contains a single spin flip, which is localized on the sample. This is 
the only realistic means to create a single spin flip state, as the alternative would 
involve polarizing the neutrons in opposite directions while taking = |0). While 
the latter is certainly feasible, the former is decidedly not. The form of l^/'*) is 
familiar from linear spin wave theory |105j : indeed, equation (13.131) represents the 
zero-momentum single magnon eigenstate of the Hamiltonian (13. 5p . In zero field, 
this is degenerate with the ground state, but a finite field breaks the degeneracy. 

Now, comparing equations (13. 7p and (13.91) shows that, for \iPq) 



a\0) 
|0). 



(3.14) 
(3.15) 



In its most general form, this state contains zero, one and two-spin-fiip components. 
Both spin flips are localized on the neutrons, while the sample is in a fully polarized 
state. In zero field, this is one of the degenerate ground states of the system. In a 
finite field, it's energy is determined by the relative magnitude of J and B^. 
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In the time interval between state initiahzation and the arrival of the first pulse, the 
sample undergoes free evolution under the effect of fl3.5p . When the pulse reaches 
the sample, all neutrons but one pass through without scattering. The neutron in 
question remains coupled to the sample for a time r, and the interaction generates 



sample undergoes a second period of free evolution, which ends with the arrival of 
a second pulse. The second scattering event is identical to the first, and has the 
effect of redistributing the entanglement in the system. As a result, part of the 
entanglement between the first neutron and the sample is swapped to the second 
neutron. 

The next step involves establishing whether the two scattered neutrons have acquired 
a correlation as a result of the scattering process. This can be achieved by measuring 
a witness operator W, as defined in chapter O Let us call the density operator of 
the scattered neutrons pn', if this state is entangled, the eigenspectrum of its partial 
transpose should contain a negative eigenvalue A_, corresponding to eigenvector 
|e~). The witness operator is then defined by 



such that Tr [VT^'^pn] = Tr [W^p^-*] < 0. It will be shown in chapter [5] that, when 
the protocol has the desired outcome, this witness can be measured with as few as 
three device settings. 

The protocol is deemed successful if the entanglement measured by the witness 
exceeds some threshold value £. In general, this depends on two factors: (i) creating 
at least S units of entanglement between the first neutron and the sample during the 
first scattering event, and (ii) swapping the state of the second neutron and the state 
of the sample during the second scattering event. Point (i) follows from theorems 
on entanglement monogamy [l6], according to which the total entanglement of the 
first neutron with the rest of the system must remain constant throughout the 
protocol. Therefore, if the two neutrons are found to share X ebits of entanglement 
after the second scattering event, the first neutron must have shared that same 
amount of entanglement with the sample prior to this event. Point (ii) implies the 
second scattering event must change the distribution of entanglement within the 

■^The mechanism whereby this occurs will be explained in more detail in chapter [5l 




When the neutron departs, the 



W = \e-){e 



(3.16) 
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system. In an ideal scenario, the second neutron and the sample would experience a 
complete role reversal; in practise only a partial swap takes place, because it is not 
possible to tune neutron interaction times separately (see section EH]). Therefore, 
the entanglement between the neutrons is typically less than X, owing both to the 
residual correlation between the first neutron and the sample, and to the creation 
of correlations between the second neutron and the sample. 

3.5 The Hamiltonian 

During each phase of the protocol described in section 13.41 the evolution of the 
system as a whole is governed by a specific Hamiltonian. In phases 1 and 3, this is 
just the free Hamiltonian of the sample, given in equation (13.51) . Phases 2 and 4, on 
the other hand, require a Hamiltonian that describes both the intrinsic dynamics of 
the sample and its interaction with the incoming neutron. If the initial state of the 
system is fully separable, the Hamiltonians of phases 2 and 4 must not commute, else 
no correlation can be established between the scattered neutrons [lU6j . Therefore, 
the sample needs to 'know' which neutron arrived first; in other words, the sample 
must have some memory. 

Let us label the Hamiltonians associated with phases 2 and 4 with Tii and 0,2, to 
indicate the neutron involved in that particular phase. Then 



where Vi and V2 represent the coupling between the neutron and the sample, and 
0£){r) is a heaviside function defined to be unity inside the sample and zero elsewhere. 
This accounts for the finite size of the sample. The dominant terms in the potential 
are the short-range, strong interaction with the nuclei, and the long-range magnetic 
dipole interaction with the unpaired electrons. The potential experienced by neutron 
m can then be written as 



ni = {Ho + Vi)9D (r), 
n2 = {Ho + V2) 60 (r) , 



(3.17) 
(3.18) 
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^6{r- r^) - firn ■ B(r) 
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The first term represents the point-hke, spin-independent coupling of the neutron 
to the nucleus at position rj, often referred to as the Fermi pseudopotential, and the 
second term is the spin-dependent interaction of the neutron's dipole moment 
with the magnetic field created by the sample electrons. The quantity bj is known 
as the scattering length. This measures the strength of nuclear scattering, and must 
be determined empirically for each nuclear species. For the remaining terms, one 
has 

Mm = fi'TV/iJV^, (3.20) 

n 

S-^EVx(^). (3.21) 

where is the spin of neutron m, Sj is the spin of the electron at position r^, and 
Rj = r — Tj is the displacement of the neutron from this electron [97] . 

The state of the incoming neutron can be represented as a plane wave with momen- 
tum k carrying spin (Tm 

\r,) = \ka^)^^e^-\r), (3.22) 

where {i/)"') is the spin state of the neutron, as defined in the previous section, and 
L is an arbitrary normalization length. As we only consider neutrons scattered into 
a specific Bragg peak, it is convenient to integrate the potential over the spatial 
coordinates of the neutron. This gives a spin- dependent spatial matrix element of 
the potential, V^^. Note this does not in itself represent the action of the potential 
on the spin state of the system. The true form of the spin potential warrants a 
separate discussion, and will be derived in the following chapter. 



Neglecting for the moment the subscript m, for a cubic sample of side D centred at 
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the origin one finds [971 1107] 



1 

L3 



all space 

2 N 



m 



N N 

b.e'"^'"' + A ^ s • [Q X (s, X q)] e'^" 



27rn' 



m 



2 





T»r Sa 


Q D 






2 



Af TV 

^ 6j + A ^ s ■ [Q X (^Sj X Q 



(3.23) 



where Q is the direction of the neutron scattering wavevector Q 
the hyperbolic sine function 



kf, Sa I 



IS 



Sa 



Q D 



sine 



smc 



QyD 



smc 



(3.24) 



and A is related to the magnetic coupling strength, whose precise value will be 
discussed in the following chapter. The last step of equation (13.231) follows from the 
Bragg condition, which states that peaks in the intensity of the scattered beam are 
obtained when the momentum transfer is a reciprocal lattice vector. 

If the sample is composed of a single nuclear species, the short-range term in the 
scattering potential can be approximated by a constant. Therefore, V* depends only 
on the direction of the momentum transfer. Expanding the cross product in equation 
f l3.23p shows the component in the direction of Q must disappear. Assuming for 
simplicity that Q is parallel to the z-direction, one finds 



L3 



Vo + A 




(3.25) 



where J2f=i ^} is the ^-component of the total spin of the sample, and Vq 



Ylf=i bj- Note this is simply a spatial matrix element of the total interaction, 
not the spin potential acting on the system at finite momentum transfer. 

Special care must be taken if the momentum transfer is zero, because the spin- 
dependent part of V'^ is undefined. To account for the rapid variation of Q in the 
vicinity of Q = 0, one can replace the double cross product in equation f l3.23p with 
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its average value over a small sphere S of radius e <^ r^^, centred on Q = 



JgQx (^Sj X Q j dQ 



(Q X (s, X Q)) = j^dQ • ^^-^^^ 



By carrying out the substitutions 



Q X (^Sj X Qj = Sj -\sj\ cos6'Q, (3.27) 
dQ = dQ sine d9d(j), (3.28) 

with IQp = 1, Q G [0, e], ^ G [0, vr], and G [0, 27r], it can be shown that 

(Q X (s, X q)) = - ^\s^\z = ^ Sj. (3.29) 

The matrix element of the potential at Q = can therefore be approximated by the 
expression 



L3 



2A / ^ jv N 

j=i j=i j=i 



(3.30) 



3 

As we will see (cfr. section . this gives a true spin potential V^'''' of the form 

2A 



^ N N ^ \ 

sl. E 4' + ^™ E + E ^.0 ' (3-31) 
^ i=i j=i i=i / 



which represents an exchange interaction between the neutron and the sample. Phys- 
ically, Q = describes forward scattering, i.e. scattering in the direction of the 
incoming momentum. 

Equation (13.311) represents an interaction capable of creating spin excitations in the 
sample. This process costs energy, which in zero external field is detracted from the 
kinetic energy of the neutron. It may then seem counterintuitive that scattering at 
zero momentum transfer - that is, elastic scattering - creates spin flips. To resolve 
the controversy, we recall that V^"^ is an average over a range of Q around Q = 0. 
The averaging procedure accounts for the fact that a Bragg peak is not infinitely 
sharp, an effect owed to the finite size of the sample. Averaging over a range of 
momenta centred at a certain Q is therefore equivalent to detecting all neutrons 
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scattered into a Bragg peak of finite width about Q. The width of the peak is 
of order D~^, as expressed by the range of the sine function of equation (13.231) . 
The variation of the Hamiltonian over this range is neghgible, hence the effective 
Hamiltonian is well-approximated by its value at the centre of the peak. 

3.6 Quantum State Evolution 

The final point I would like to make regards the method I use to derive the scattered 
state of the system. As the initial state \iIjq) is pure, the scattered state can be 
obtained by straightforward time-evolution using the canonical operator f/(7i, t) = 
exp [—iHt]. This operator must refiect the four stages of the protocol, hence UiH, t) 
can be expressed product of four terms 

U{n, t) = U{7i2, r)U{Ho, r})U{Hu r)U{Ho, r/), (3.32) 

describing phases 4, 3, 2 and 1, respectively. None of the terms commute (save, 
trivially, the second and fourth), so U(T-C,t) truly represents a specific sequence of 
events. It is furthermore evident that, in zero field, the first and second periods of 
free evolution only contribute a phase to the scattered state. During phase 1, this 
is true in finite fields also, because both I^Pq) and I'iPq) are eigenstates of Hq. The 
scattered state \ipf'^) is then, modulo a phase factor 

\^l:f'^) = U{n2,T)U{Ho,r'f)U{nuT)\ij^'''). (3.33) 

From here, the state p„ of the neutrons can be obtained by calculating the system 
density matrix p = \^f''^){4'f'^\ and tracing out the sample. 

Ultimately, the outcome of the protocol is determined by the neutron polarization 
{a, (3) and the interaction time r. It is therefore essential to relate these parameters 
to the experimental settings of our equipment. The physical interpretation of a 
and (3 is immediate, as they are simply the direct numerical expression of the real 
property that is the neutron polarization. For r, on the other hand, the situation is 
a little more hazy, because the interaction time between the neutron and the sample 
is not in itself a well-defined quantity. In fact, scattering problems are not typically 
solved in terms of state vectors. Conventionally, one works with cross sections, 
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which represent the probabihty that a certain scattering event will result in neutrons 
being detected in a specific region of space, and with a certain polarization. The 
wavefunctions \i/ of the scattered neutrons satisfy the Lippmann-Schwinger equation 

m 

^ = $ + Gy^, (3.34) 

where $ is a solution to the time-independent Schrodinger equation, V is the scat- 
tering potential and G is known as the retarded Green function. To obtain these 
solutions, one acts on the initial state of the system with an operator known as the 
S'-matrix, given by 

S = I-2mTS{E' - E), (3.35) 
where T is the transition matrix, related to the scattering potential asT = V {1 — GV) 

The S'-matrix is universally acknowledged as providing the correct description for the 
observable outcome of a scattering experiment. It follows that any results obtained 
with a quantum-mechanical treatment should measure up to the predictions of the 
S'-matrix method. Comparing the outcomes of the two techniques may therefore 
provide an indication of the relationship between r and real, physical quantities one 
might be able to measure and tune. I will explore this possibility in the following 
chapter. 



3.7 Conclusions 

In this chapter, I have described an original protocol aimed at creating entangle- 
ment between two distinct neutrons. The success of the protocol would provide the 
quantum information processing community with a novel source of entanglement, 
as well as a playground for the test of quantum effects with neutrons. 

The present scheme bears some superficial similarities with a proposal for entangle- 
ment extraction from a solid, advanced some years ago by De Chiara et al. [81]. As 
the name suggests, the protocol involves sending two uncorrelated probe neutrons 
to interact separately with two entangled spins in a bulk solid. It is shown that, if 
the neutrons and the sample spins are coupled by an interaction similar to (13.311) . 
for some values of the interaction time it is possible to transfer the entanglement 
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of the sample spins to the probes. In addition, in the hmit of multiple collisions, 
the entanglement between the probes approaches unity. The authors proceed to 
describe a more realistic situation, in which each probe neutron, being a wavepacket 
of finite width, interacts with a subset of spins in the sample, which is prepared in a 
state analogous to fl3.13p . In this case, the probes are entangled only if the separate 
spin subsets are entangled; in other words, the success of the protocol relies on the 
presence of entanglement between different parts of the sample. 

Three main points distinguish this protocol from the scheme I have suggested. First, 
the interaction of the neutrons with the sample is sequential, meaning the respective 
interaction Hamiltonians do not commute. This is not the case in [HI], because 
the probes interact simultaneously with spatially separated regions of the sample. 
Second, I assume the incoming neutron couples to all the sample spins, rather than 
a subset. This creates a three-body entangled state of the neutrons and the sample, 
which for certain values of the interaction time is close to being separable with 
respect to the neutrons/sample bi-partition. Finally, and most important of all, no 
prior entanglement is necessary for the protocol to succeed: the scattered neutrons 
can become entangled even if the sample is prepared in a fully separable state. This 
allows for greater flexibility in state initialization, and in general renders the protocol 
a more versatile and realistic alternative to that of |81j . 



Chapter 4 



The Formal Theory of Neutron 
Scattering 

In this chapter, I illustrate the relationship between time- dependent and time-independent 
treatments of scattering as applied to the entangling protocol discussed previously. 
The time-independent counterpart of the scattered state \'ipf'^) will be calculated for 
both Q = and Q 7^ scattering geometries. This form will be compared to an 
expansion of equation i\3. 33\) in terms of the acting potentials to yield the physical 
significance of the time parameter r in this context. A review of the traditional, 
time-independent scattering formalism developed by Heisenberg, often referred to as 
the S -matrix method, can be found in Appendix\M and follows the work of Gell-Mann 
and Goldberger I109f . 

4.1 A Preliminary Study 

To first order in the scattering potential, the S'-matrix representing scattering be- 
tween initial and final eigenstates of the free Hamiltonian of the system is 

S^^^ = I-2mV6{Ei- Ej) , (4.1) 

where Ej and Ei are the initial and final eigenenergies of the free Hamiltonian, 
respectively. In the case of forward scattering, the derivation of scattered states 
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according to the ^-matrix formalism reduces to solving the familiar problem of a 
particle of energy E traversing a potential barrier of height V < E. Let us then 
consider a one-dimensional Gaussian wavepacket of width w, centred on position 
Xq and momentum /cq, traversing a region of uniform magnetic field of strength 
oriented along z, as shown in figure HTH 




Figure 4.1: Transmission and reflection of a wavepacket of incident energy E from a barrier of 
height B < E. 



At time t, the state of the wavepacket at any point along the x axis is given by 

cr k 

where the Ek represent the energies of each /c-component, corresponding to space 
and spin eigenf unctions 4'k,a- The (j)k,a are easily (if somewhat laboriously) derived 
by expressing the wavefunctions outside and inside the barriers in terms of incident, 
transmitted and reflected waves, and matching wavefunctions and derivatives at the 
boundaries. One has 

<Pi^={e^''' + R^e~-^'^)\a), (4.3) 
<j>l^={T,e''^-^ + Ke-''^-^)\a), (4.4) 
0t = ne''1^), (4.5) 



Chapter 4. The Formal Theory of Neutron Scattering 



47 



with 



R(7 



(4.6) 
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(/c^ + ql) sin g^L + 2iqf^k ' 
2A; (A; + q„) 



(4.7) 



e2ig.L (A;2 _ g2) _ (^2 + ^2) ' 



K 



(4.8) 



e2i,.L (^2 _ ^2) _ (^2 + ^2) ' 



(4.9) 



g2*9.L (^2 _ ^2) _ (^2 + ^2) ' 



where L is the extent of the field region, and R'^ are the spin-dependent refiec- 
tion coefficients in regions 1 and 2, T„ and are the spin-dependent transmission 
coefficients in regions 2 and 3, and q^ is the spin- dependent wavevector inside the 
barrier. In natural units, this is given by 



for particles of spin up and down, respectively. Values of k can be fixed by imposing 
periodic boundary conditions over some length C ^ L. For a set of integers m, one 
then has 



The group velocity of each km component is a linear function of km] therefore, the 
wavepacket spreads as it travels, as shown in figure 14. 1[ Let us assume that, at 
some time t = 0, the wavepacket is located at a point Xj„ to the left-hand side 
of the barrier. After a time tg, we measure the position of its peak and find it to 
be at Xout, on the right-hand side of the barrier, tg then defines the time required 
for the wavepacket to travel a distance \xin — Xout\- Provided Xj„ is sufficiently far 
from the barrier, the spin component of the wavepacket is initially oriented along 
the positive x axis. By the time we reach Xout, the spin has been rotated to some 
new orientation, which can be estimated by calculating the expectation values of 
the Pauli operators a^, <Jy and from the density operator of the system at time tg 



q^ = {e-2B,y, 
q^ = {e + 2B,y\ 




(4.10) 




(4.12) 



(o-q) = Tr [aaP (ts)], 



(4.13) 
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with a G {x,y,z}. The square modulus of the wavefunction at therefore de- 
fines the probabihty of finding the wavepacket at position Xout with a certain spin 
orientation {9,(f>), defined by I4.13[ 

Now, let us put aside the spatial part of the wavepacket's evolution, and focus 
on spin. Changes in the coordinate of the spin presumably took place while 
the wavepacket was passing through the field. We can estimate the approximate 
duration of this phase from the wavepacket's average velocity inside the barrier 



Is it then valid to identify Ib with the time taken by the spin to rotate from {it/2, 0) to 
{6, 0), or is there some other timescale associated with spin rotation? In other words, 
do the 'space' and 'spin' clocks tick at the same rate? To answer this question, I 
calculate how much a single spin initially oriented along the positive x axis precesses 
about z under the effect of the field B^z. If this value is compatible with the 
(p angles given by the {<Ja), there is good reason to assume the 'space' and 'spin' 
clocks measure time in the same way. 

Figure 14.21 shows a few examples at fixed C, L, and xq. Here, the three axes 
represent the eigenstates of ax-, Uy and cr^ in the Bloch sphere, and the red vector 
indicates the spin orientation of the wavepacket at time t = 0. The green vector 
shows the expected precession around z over a time t^? and the orange segment 
illustrates the range of possible spin orientations over the width of the peak. Note 
that the states of the system at time ts are not always pure, hence the magnitude of 
the vectors forming the orange segment can be less than unity. It is clear from this 
figure that for the chosen values of the input parameters one finds good agreement 
between the 'space' and 'spin' clocks of the system. 

What remains now is to verify whether this holds in a more general context also. 
Therefore, I return to the entangling protocol, and to a more formal approach, by 
deriving the scattered state of the system according to the ^-matrix formalism. 
For simplicity, I will assume elastic scattering at both zero and finite momentum 
transfer. As discussed at the end of section 13.51 under these circumstances spin flips 
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can still be occasioned by neutrons scattered slightly away from the centre of the 
Bragg peak. 



2 




5 10 15 20 25 

Is 



(b) 




Figure 4.2: A comparison between the 'space' and 'spin' clocks, for the situation shown in figure 
14.11 In all cases £ = 40, L = 10, xq = —8, ko = ^j^. It is assumed the barrier originates at 
position a; = 0. (a) The im-normalized probability distribution j'l'p as a function of travel time tg 
at position Xout = 16 for w ~ 2 and ~ 0.05. To the right, the distribution of measured spin 
states over the range Atg = [19, 25], compared with the initial spin state (red line) and the expected 
final spin state (green line), (b) The un- normalized probability distribution j^Pp as a function of 
travel time is at position Xout = 19 for w = 3 and = 0.02. To the right, the distribution of 
measured spin states over the range Ats = [21,31], compared with the initial spin state (red line) 
and the expected final spin state (green line). All quantities are expressed in natural units. 



4.2 Time at Zero Momentum Transfer 



I will derive the scattered states of the system using only the first-order expansion 
of the S'-matrix, as expressed by equation (14. ip . At zero momentum transfer, the 
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S'-matrix representing the first scattering event is given by 

S = l-2mVi6{E' - E) , (4.16) 

where Vi is the scattering potential of equation (13.191) . and E is the energy of the 
initial state of the system. For incoming neutrons with momentum ko, given the 
dispersion relation E = this can be written as 

S = l-2mV,^-^^^, (4.17) 
k 

where ko and k are the magnitudes of the neutron's initial and final momenta, 
respectively. The general state of a neutron with momentum k and spin am is given 
by equation (13.221) . Neglecting for the moment the subscript m, the ort honor mality 
of the neutron states obeys the relation 

(k'a'lka) = (^^^ 6ik' - k)6^^,, (4.18) 

where 6{k' — fc) is a three-dimensional delta function in momentum and [^)^ is 
the number of states per unit volume of momentum space. The orthonormality 
condition of equation (14.181) also applies to the states of the system as a whole, 
provided one now takes the ket \a) to label the overall spin state of neutron plus 
sample. Indeed, as the spatial part of the sample's wavefunction is not changed by 
the scattering event, the sample is entirely described by its spin state. 

Following section 13.51 we focus on neutrons scattered into the Bragg peak at Q = 
and calculate the spin state produced by the action of the S matrix by integrating 
the scattered state over the spatial coordinates of the neutron. Recall, however, that 
the peak has finite width, due to the finite size of the sample. Our intention is to 
collect all neutrons scattered into the peak, not just those precisely obeying kj = k/-. 
As measurements will be made after both neutrons have scattered, it is in our 
interest to maintain a certain amount of uncertainty in these measurements, because 
the spatial and spinor parts of the wavefunction of scattered neutrons are usually 
entangled [llOj . Therefore, an exact measurement of momentum would destroy any 
indeterminacy in the spin state of the neutron, proving counterproductive to the 
success of the entangling protocol. To account for this effect, we take the outgoing 
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spin state of the system to be an average over all outgoing spatial states q consistent 
with scattering into the Q = peak 



m) 



L 
2^ 



{q\S\ka) dq, 



(4.19) 



which becomes, with the help of (14.171) 



l*?) = k>-S/'{q|V,M^<i^dq^ 



(4.20) 



Note the function (g|Vi|A;o") differs from the matrix element of equation (13.231) be- 
cause it is evaluated at a generic value of Q within the Bragg peak, rather than at 
its centre. 

Let us assume the scattering length bj is a constant for all nuclei j, and set 



m 



N N 

b-e^Q-r, + A 5^ s ■ [Q X (s, X q)] e^Q-^ = / (Q) 5 (Q - G) , (4.21) 
j=i j=i 



where G is a reciprocal lattice vector. The three-dimensional delta function 5 (Q — G) 
accounts for the fact that the right hand side of equation (14.211) is small unless the 
momentum transfer coincides with a reciprocal lattice vector. From equations (I3.17P 
and (I3.23p . the matrix element (g|Vi|A;cr) then becomes 



(g|Vi|M = ^ 



/(Q)5(Q-G)*Sa 



Q D 



JjfiG) Sa 



with Sa [■] given by equation (13.241) . and 



(Q-G) D 



(4.22) 



/(G) 



(4.23) 



To within a constant, the function / (G) represents the spatial matrix element of 
the potential evaluated at the centre of the Bragg peak. At G = 0, the form of this 
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matrix element is known from equation (13.301) . With hindsight, let us then set 
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(4.24) 



Substituting equations (I4.22p and (14.241) into (14.201) . one finds 



Q D 



S{q - ko) 



dq. 



(4.25) 



At fixed kj, the integral of equation (I4.25P can be written in terms of the momentum 
transfer as 



Sa 



Q D 



ko-Q 



(4.26) 



where 



Qo = q- ko, 

q = ko - Q. 



To define a convenient and visually descriptive set of integration coordinates, I make 
use of a construction known as the Ewald sphere, shown in two dimensions in figure 
14.31 This is a sphere of radius /cq, centred at the origin of reciprocal space. If the 
scattering process is elastic, the tip of the final wavevector kj must lie somewhere 
on the surface of the sphere, as ko = kf. The momentum transfer Q is the vector 
joining ko and kj. The integration region dQ can then be approximated by a small 
cube of volume dQadQpdQ^ centred on the tip of the vector kj, where the directions 
{a, (3} and 7 lie parallel and perpendicular to the surface of the sphere, respectively 
(fig. 14.31) . In these coordinates, the integral of equation (I4.26P at fixed ko becomes 



Sa 



Sa 



Sa 



ko — Q 



dQadQ/sdQj. 



(4.27) 



The integration with respect to is immediately resolved by the delta function 



Sa 



SjQo) 
ko - 



dQ^ 



1 

ko' 



(4.28) 



The remaining terms must be integrated over the range of the sine function, which 
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Figure 4.3: The Ewald construction in two dimensions (left), and the set of integration co- 
ordinates a, P and 7 (right). The left-hand diagram is adapted from a figure taken from 
http://perso.fundp.ac.be/ jwouters/DRX/ewaldl.GIF 



Let us now consider the scattering event in a time-dependent frame. If the neutron 
is scattered by a spin- dependent potential in a time r, the spin state of the system 
at time r is described by 

= e-*^"|a) ^ (I-zl/r) (4.31) 

with h = 1. Comparing equations (14.30p and fl4.3ip shows that first-order agreement 
between the time-independent and time-dependent pictures is achieved when 

V = V''% (4.32) 

(4,33) 

The form of V^^'^ is known from equation fl4.24p : the quantity ^ is simply the free- 
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flight time of tlie neutron througli the sample. These results lead us to conclude 
that: 

1. The potential that determines the evolution of the spin state of the system is, 
from equations ([SISD, dSHD and (OD 

N 

(ij) i=i 

2A / ^ ^ ^ \ 

V j=i i=i i=i / 

2. To first order in the scattering potential, the interaction time of a forward- 
scattered neutron coincides with its free-flight time through the sample. This 
can be tuned by simply adjusting the neutron momentum, thus providing a 
realistic means of setting the interaction time to the desired value. 



4.3 Time at Finite Momentum Transfer 

The derivation of the scattered state at finite momentum transfer Q = Q^z is 
virtually identical to that outlined in the previous section. What changes is the 
form of the scattering matrix, now given by 

5 = -2vrzV,fcM. (4.35) 

rv 

The disappearance of the first term reflects the absence of forward scattering; it is 
now assumed all neutrons scatter away from the incident direction. The scattered 
spin state, averaged over the width of the Bragg peak, is then 

I^J) = -^ /(q|Vi|M^^i-M,q. (4.36) 
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Let us repeat the integration of equation fl4.22p . To within a constant, the form of 
/ (G) is known from equation (13.251) . It is therefore convenient to set 



(N ^ \ 

j=i j=i J 



(4.37) 



Substituting equations fl4.22p and (14.371) into (14.361) . one finds 



47^2 



V'^y / Sa 



Q D 



6{q - ko 



dq. 



(4.38) 



This can be resolved as before, with the help of equations (14.261) . (14.271) . (14.281) and 
Km . to give 



N ^ \ D 

i=l j=l J 



a). (4.39) 



Owing to the lack of a non-imaginary component proportional to the identity, this 
cannot be immediately traced back to an expansion of the time evolution opera- 
tor to any order in the scattering potential. One is therefore left with a series of 
possibilities, among which: 



1. The neutron interaction time is not a well-defined quantity for this process; 

2. The correct interpretation of the interaction time emerges from a higher order 
expansion of the 5* matrix; 

3. The present model is too simplistic; the presence of an apparently complex- 
valued time parameter may then result from failing to take into account addi- 
tional processes, such as decoherence. 

Related to this final point, the significance of imaginary time in quantum mechanics, 
and particularly in scattering theory, is a topic of intense debate, which has yet 
to find a unique resolution. It has been suggested that the complex part of the 
time parameter arises from changes in the scattering probability with potential. 
Specifically, Im[r] is thought to describe the variation in the number of particles 
scattered per unit time, per unit change of the scattering potential [111] . The 
derivation of this result is based on a spin-independent potential, therefore it may 
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not be of strict relevance to the present discussion. However, it is interesting that 
quahtatively similar findings should crop up also in the context of tunneling [1121 

[IS]. 

We can conclude that scattering at finite momentum transfer cannot immediately 
be described in a physically intuitive time- dependent frame. Consequently, it is 
difficult to see how parameters such as the neutron interaction time may be tuned 
in an experimental setting. It may be that the entangling protocol of the previous 
chapter translates into a feasible scheme at zero momentum transfer only. Therefore, 
I will restrict subsequent analysis to realizations of the protocol at Q = 0. 



4.4 The Magnetic Coupling Strength 

We conclude from sections 14.21 and [4.31 that in the limit of weak scattering (T ^ V), 
one can ascribe a well-defined physical meaning to the time-evolution of the spin 
part of the system wavefunction only in the case of forward scattering. The Hamil- 
tonian representing this process resembles an exchange interaction, whose strength 
we define as A^;. To complete our understanding of the scattering process, I will 
now discuss the value of this parameter. 

The Hamiltonian of fl4.34p is derived by taking a spatial matrix element of the 
potential of equation (13.191) . Let us isolate the magnetic part of this potential, 
writing 

Vb = -lira ■ B(r) On (r) . (4.40) 

We then have 

{k\VB\k) = --^ / e-'^^r^ [m™ ■ B(r) Od (r)] e^^-Mr = 

^ J all space 

= ■ j e^Q-'-B(r) Od (r) dr. (4.41) 
At zero momentum transfer, the exponential term reduces to unity, leaving 

{k\VB\k) = -ji- [ B(r) do (r) dr. (4.42) 

J all space 
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The integral on the right-hand side of this equation is simply the average magnetic 
field of the sample multiplied by its volume. One can then write 

{k\VB\k) = -^{fim-B). (4.43) 

Comparing equations fl3.23p . (13.311) . (14.341) and (14.431) . it becomes clear that for 
Q = 

2A ^ ^ 



j=i j=i j=i 



Expanding the left-hand side of this equation in terms of the average magnetization 
M gives 

- /x^ ■ B = [Sm ■ M) . (4.45) 

Then, recalling the magnetization of the sample is defined as the magnetic moment 
per unit volume, we obtain 



-lXm B = -— \sm-^sA. (4.46) 



N 



It follows from equations (14.441) and (14.461) that 



The action of the spin Hamiltonian V^""^ on the combined spin state of the interacting 

4 



neutron and the sample produces terms proportional to The coupling energy A^; 



= ~8 ^liV ' ^^-^^^ 



can therefore be defined as 

3 gNfJ'NjM)9efJ'B 

where the sample volume has been replaced with the number of spins multi- 
plied by the volume of the unit cell a^. We finally conclude that the Hamiltonian 
describing the evolution of the system spin state during forward scattering of a 
neutron is given by 

3A 4 ^ ^ 



^"^ = ra7^ + ^4.49) 

j=i j=i j=i 
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with 



A = 3 , (4.50) 



where the value of A can be interpreted as the strength of the magnetic interaction 
of a neutron with an electron at a distance of one lattice constant. 



4.5 Symmetries and Operator Representation 

The full Hilbert space of a system comprising sample spins and two neutrons 
has dimensions d = 2^+^. Fortunately, the need to deal with exponentially growing 
matrices is removed by exploiting the symmetries of equation f l4.34p . and the periodic 
nature of the sample. Let us define the operators Sz and Sz, which represent the 
2;-components of the total spin of the sample and of the system, respectively. It is 
easily shown that 

[Ho, §,] = [Ho, 5J = [V^ §J = 0, (4.51) 

therefore both the free Hamiltonian of the sample and the interaction potential have 
excitation-number conserving symmetries. Bearing in mind our chosen initial states 
contain at most two spin flips, this allows us to define a reduced computational basis, 
comprising all possible zero-, single- and double-spin-fiip states. The dimensions of 
the reduced basis can be calculated by summing the binomial coefficients (^n+2)Co, 
{N+2)Ci and (Ar+2)C'2, to obtain 

d=^ + ^ + 4. (4.52) 

This grows quadratically rather than exponentially with A^, but still blows up as 
becomes large. 



A further reduction is possible if one abandons the canonical basis. The single- and 
double- spin-flip eigenstates of a periodic sample can be expressed in terms of Bloch 
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vectors \s 



I) and \sl) 





J2 e''^-(^^+r')|j/). 



N 



(4.53) 



(4.54) 



iV(A^- 1) 



As the initial states of the system are translationally invariant, and neither Hq nor 
V"* contain any spatial dependence, coupling to states \sl) and \sl) is null unless 
k = 0. The Hilbert space of the sample is, therefore, fully specified by the vectors 
|0), \sl) and \ sq). Combining these with the four canonical two-neutron basis states, 
it is possible to define a complete 8-dimensional basis for the system, comprising the 
following eigenvectors 



where |000) = |0)2|0)i|0), and so on. Unless otherwise stated, from this point on- 
wards equations fl4.55p - fl4.62p will be my computational basis. For ease of reference, 
I call this the Block basis, and drop all indices to simplify notation. Henceforth, the 
values inside bras and kets should be interpreted as labeling \second neutron, first 
neutron, sample), respectively. 

The dimensions of the Bloch basis are fixed, manageable, and independent of the 
size of the sample. As a result, many calculations can be done analytically and 
operator representation is greatly simplified. The initial states of the system, for 



1) 
2) 
3) 
4) 
5) 
6) 
7) 
8) 



000), 

001, ), 

002, ), 
010), 

Oils), 

100), 

101s), 

110), 



(4.55) 
(4.56) 
(4.57) 
(4.58) 
(4.59) 
(4.60) 
(4.61) 
(4.62) 



Chapter 4. The Formal Theory of Neutron Scattering 



60 



instance, can now be written as 



|2), 

a^\l) + a(] (|4) + |6)) + /32|8). 



(4.63) 
(4.64) 



Also, the zero field free Hamiltonian becomes a multiple of the identity 



Ho = -JNI, 



(4.65) 



owing to the fact that the sample basis states |0), |sq) and |sq) are degenerate 
eigenstates of i^o with eigenvalue —JN. Equation (14.651) underlines a very important 
point: the internal coupling of the spins in the sample plays no meaningful role in 
the dynamics of the system. Hence, both the numerical value of J and its sign 
are irrelevant. This warrants two observations. First, provided one could prepare 
the same initial states, ferromagnetic and antiferromagnetic samples would produce 
exactly the same outcome. Second, in zero field the evolution of the system between 
scattering events is trivial. Therefore, if one did succeed in producing a correlation 
between the neutrons, this correlation would not depend on the time elapsed between 
the departure of the first neutron and the arrival of the second. Qualitatively, this 
means the 'information' deposited in the sample by the first neutron can survive 
indefinitely; in other words, the sample functions as a perfect quantum memory 
between scattering events. Note this is a specific feature of the present model, and 
should not be interpreted as a result of general validity. 

4.6 Conclusions 

The scattering of a neutron from a sample is usually described as a time-independent, 
probabilistic process. In the case of forward scattering, it is possible to recast the 
problem in a temporal frame, thus showing that the interaction time of a neutron 
with a scatterer is, in the limit of weak scattering, the time taken by the neutron to 
traverse the sample. At finite momentum transfer, the interpretation of r remains 
unclear. It emerges from a first-order calculation that r may be a complex-valued, 
or even completely imaginary quantity. This possibility has been the focus of myr- 
iad studies, but no unique resolution has been found to date. Indeed, a unique 
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resolution may not exist, as the details of the scattering geometry or the scattering 
potential may be a factor. Some answers may lie in a higher order treatment of 
the problem; alternatively, it may be necessary to look into the dynamical theory of 
neutron scattering, which accounts for the possibility of the incident neutron being 
re-scattered within the sample before emerging [110] . These options are not explored 
in the present thesis, but may become the topic of future work. 



Chapter 5 



Two-Neutron Entanglement via 
Scattering 



This chapter considers two alternative realizations of the neutron entanglement 
scheme outlined in chapter In order of discussion, these are: (i) scattering at 
zero momentum transfer from a sample in a single magnon state; (ii) scattering 
at zero momentum transfer from a sample containing no spin excitations. In each 
case, the performance of the system in zero and finite field is analyzed. After de- 
scribing the quantum state of the system at each stage of the protocol, I will discuss 
the behaviour of the entanglement between the neutrons as a function of the input 
parameters. It will emerge that, for certain combinations of these parameters, the 
concurrence attains high values. Furthermore, the concurrence proves independent 
of the time of free evolution rj between scattering events. To investigate this prop- 
erty further, I will compare the time- evolution of the concurrence to that of the 
logarithmic negativity, which is not independent of rj. I then discuss parallels with 
the optical entanglement scheme proposed by Haroche et al. (TJ^, and evaluate the 
robustness of the protocol and its experimental feasibility in light of current neu- 
tron scattering facilities. Finally, I propose a simple model to describe how sample 
anisotropics might affect the performance of the protocol. 



62 
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5.1 Scattering from a S ingle- Magnon State 

Let us first consider forward scattering from a single-magnon state, taking the input 
wavefunction of the system to be \iPq) = |001) = |2). One sees from equation 
f l3.33p that the first non-trivial step of the protocol is the first scattering event. It 
was shown in sections 14.21 and 14.41 of the previous chapter that the Hamiltonian 
describing this process is 



It is clear from this equation that both Vq and the exchange coupling between the 
sample spins simply map the state of the system onto itself. In the interaction 
picture, Hi can then be written as 



N 




(5.1) 



(5.2) 



with 



N 




(5.3) 



(5.4) 
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Absorbing the factor of 3/8 into the value of A, the non-zero matrix elements of Ti'i 
are 



{2\n[ 
{3\n[ 



1) 

4) 
2) 

3) 

5) 

4) 

5) 



{6\n[\6) = NB, + \, 
{8\n[\8) = NB, - A, 



(7|H;|7) = A(l--)+i?,(iV-2), 



A(l--)+5,(iV-4), 



A( ^-1) +5.(iV-2), 



2A 



N 



{5\n[\3) = 2X. 



N 



1 - 



N 



(5.5) 
(5.6) 

(5.7) 
(5.8) 
(5.9) 
(5.10) 

(5.11) 



The time evolution of state I'iPq) is dominated by only two of the eight eigenstates 
of this potential. We write these as |i) and |2): 



|i) = c|2) + rf|4), 
|2) = ci|2)-c|4), 



(5.12) 
(5.13) 



with 



'fi?-25. ( 1-4^ A + fl+ ^ 



N 



N 



^ = ^'%N,\,B,)=\ 



N 



B. 



v/4A^TiV 

2A 

v/4A2 + A^<^2 ■ 



A2, (5.14) 
(5.15) 
(5.16) 
(5.17) 



The scattered states of the system can then be derived by expanding the 'effective' 
time-evolution operator U {H[,t) in terms of its eigenstates \i) and eigenvalues Ei, 
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as 



uin[,T) = J2e-'''n^){^\- (5.18) 

i 

To within a global phase, the first scattered state is therefore 

1^^) = e-'^^[{c^ + d^e-^'^^)\2) + cd{l - e-^''^^)\A)], (5.19) 



where A = + B,{N - 1) - 0. 

During the second period of free evolution, each component of \ip2) simply acquires 
a phase proportional to the corresponding eigenvalue of Hq , thus yielding the state 

1^^) = e-'^r^-iB.Nr'^^^2iB.^r'^^^2 ^ ^2g-2i0r ^) 1 3^ + _ e~2i<^-) |4)] . (5.20) 

The final stage of the protocol is the second scattering event. As before, this is 
governed by an effective Hamiltonian Tig*^^ = Hq + , which is formally similar 
to Ti']^ but for the placement of signs and off-diagonal matrix elements. The three 
eigenstates of Ti^ contributing to the evolution of 1-03^) are 

k) = |4), (5.21) 
|5) = c|2) + d|6), (5.22) 
|6) = rf|2)-c|6), (5.23) 

hence the final state of the system becomes 

= e-^(2Ar+B.A^r;)[g2^B.r;^^2 ^ ^2g-2i<^r ^2 1 3^ + cde-'"^(l - e-2''^")|4) 

+ crfe''^^"/(l-e-2^'^0(c' + (iV2^'^0|6)], (5.24) 

with 

y = \[l + ^+B, + cp. (5.25) 

For brevity, let us label the coefficients of l^jf) by F, A and 6, respectively. The 
reduced density matrix of the neutrons in the canonical basis is calculated by per- 
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forming a partial trace over the system, which gives 



+ |enio)(io| 



= |rnoo)(oo| + |Anoi)(oi| 

A0*|O1)(1O| + A*0|1O)(O1|. 



(5.26) 



The concurrence of the neutrons is then calculated from the eigenvalues of the 
spin-flipped density matrix described in section 12.31 If p = p„, the only non-zero 
eigenvalue of this matrix is A = 4|Ap|0p. Substituting the values of A and 0, one 
obtains 



8V2A2sin2 (j)T 



2A2 



4A2 



sm 0T 



(5.27) 



Equation fl5.27p immediately reveals several important characteristics. First, the 
concurrence does not depend on the time between scattering events r'j. This is es- 
sential for the purpose of experimental implementation, because rj- cannot be tuned. 
In zero field, this invariance follows from the fact that the sample Hamiltonian Hq 
is a multiple of the identity. In finite field the interpretation is more complex, as we 
will see. 

Second, the concurrence is an oscillating function of time, whose periodicity is de- 
termined by the cosine function cos(20r): 



T^(iV,A,i?,) = ^ 



TT 



(5.28) 



The quantity 2(f) is the energy splitting of the eigenstates of Ti.'^ corresponding to 
the spin-flip being shared between the interacting neutron and the sample. There- 
fore, during each scattering event the system oscillates between states |2) and |4) 
(1 6)), in which the spin flip is localized on the sample or the first (second) neutron, 
respectively. This is a promising beginning, as it suggests there will be times at 
which some entanglement is present in the system. To verify how this entanglement 
is distributed, however, we must focus in detail on the behaviour of the concurrence 
as a function of time, spin density and magnetic field. 
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5.1.1 Zero-Field Evolution 



If the applied field is zero, equation fl5.27p becomes 



8iVsin^[A(l + i)r]y/iV2 + l + 2iVcos[2A(l + i) ^ 
CiN, A, r) = (N + ir • ^^-^^^ 

At fixed and A, the only free parameter in this expression is the interaction time. 
Maximizing with respect to r and neglecting un-physical solutions yields 

r- = . (5.30) 

hence in zero field the optimal interaction time is half the oscillation period. The 
peak concurrence Cp is then found to be 

Cp{N,X,T ) = ^^^^p • (5.31) 

It is evident from equation (15.311) that Cp falls as increases, converging to a value 
8/A^ as A^ ^ oo (fig. 15. ip . This phenomenon can be explained in terms of the 
concurrence of the first neutron and the sample prior to the second scattering event, 
which is 



4ViV| sin[A (1 + i) t]\JN^ + 1 + 2A^cos[2A (l + r] 
CM A, r) = )L__ . (5.32) 

The maximum value of Ci drops with A^ as A^^^. The behaviour of Cp therefore 
stems from the fact that the total amount of 'available' entanglement falls as the 
number of spins in the sample becomes large. 
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Figure 5.1: Input state iV'ff), Q = 0, = 0, X = I. (a) The evolution ol the coneurrence as 
a lunetion of interaction time r for iV = 4. = 10 and = 20 (blue, red and green curves, 
respectively), (b) The peak concurrence as a function of A^. All quantities are expressed in natural 
units. 



5.1.2 Finite-Field Evolution 

We conclude from the previous section that if no field is applied there is a limit to 
the amount of information the sample can store; but is this true in finite field also? 
And if not, is there an optimal field value? To answer this question, let us consider 
the first scattered state. The success of the protocol relies on the system being as 
entangled as possible before the second scattering event. From equation (15.191) . the 
state \ip2) is maximally entangled if 

2c2(l-c2)(l-cos20r) = ^. (5.33) 

If this is true, the derivative of equation (15.331) with respect to c must vanish. 
Therefore 

4c2(l - 2c^) (1 - cos 2(pr) = 0. (5.34) 

The only reasonable solution to this equation is c = i"^? hence d = also. 
Therefore, the concurrence of the neutrons is maximized when the eigenstates of 
the interaction Hamiltonian in the single-excitation sector correspond to the spin- 
flip being equally shared between the interacting neutron and the sample. From 
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equations fl5.14l) - fl5.17l) . one finds 



) 



(5.35) 



It follows from equation fl4.5Up that, in the large limit, the optimal field represents 
the interaction of a neutron with an electron at a distance of one lattice constant. 
Substituting the value B* into Ti'i, the difference between the diagonal elements 
{2\T-l[\2) and (4|7i'^|4) becomes zero, which ensures the states |i) and \z) take on 
the form of equal superpositions. As we will see, 5* represents the optimal field 
value for input state \iPq) also, because in both cases the first scattered state can 
be written in terms of eigenstates formally equivalent to |i) and 

Substituting = B* into equation fl5.27p . the concurrence takes on a particularly 
simple form 



For all values of A^, the peak concurrence between the neutrons is now fixed at 0.77 
(fig. 15. 2p . This can be expressed in units of entanglement by making use of the 
relationship between the concurrence C and the entanglement of formation S: 




(5.36) 



This is now easily maximized with respect to r to yield 




(5.37) 



8 




(5.38) 



h{x) 



a:log2 X — {1 — x) log2(l — x). 



(5.39) 
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Figure 5.2: Input state Ii/jq), Q = 0, A = 1. The evolution of the concurrence with interaction 
time at optimal field i?* for = 4, = 10 and iV = 20 (thick blue, dashed red and dot-dashed 
green curves, respectively). All quantities are expressed in natural units. 

Hence, for a concurrence Cp = 0.77 the neutrons share approximately 0.68 ebits of 
entanglement, and are found to be in state 

= i/xnoo)(ooi + iz.noi)(oii + ienio)(ioi 

+ i^r|01)(10| + z/*{|10)(01|, (5.40) 

with = |, |z/|2 = |, |^|2 = |, = ^ 0.385, and Arg [u^] = g(iV), where 
the function q{N) is shown in figure 15.31 
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Figure 5.3: Input state iV'cf), Q = 0, A = 1. The argument of the coefficient of state as a 
function of the number of spins in the sample, N. 
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Let us now consider the evolution of the concurrence for a generic value of B^. By 
maximizing equation f l5.27p with respect to r, one finds the optimal interaction time 
function of applied field to be 



If lies outside a certain range, tb is equivalent to half the oscillation period T^; 
hence, the concurrence oscillates sinusoidally with time [fig. 5.4(a)| . Inside this 



range, however, tb and no longer coincide and the concurrence develops a 

double-peaked structure, showing a higher frequency oscillation modulated by an 
envelope with period T^. The limits of this range are found by calculating when the 
equality tb = ceases to apply. This yields the interval Bz G [B^-, Bz+] shown 



in figure 5.4(b) , with 

B« = ^(l-^F + yi). (5.43) 

The peak concurrence outside and inside this interval then takes on the form 



8\/2A2|(^- 



2A£ 

N 



C,iN,X,Bz) = —-^ , (5.44) 

QiN,X,Bz) = (5.45) 

Equations (15.311) . (15.441) and (I5.45P then show the peak concurrence is always im- 
proved by applying a field, provided this does not exceed an upper limit of B^ = 
2A(1-^) (fig. EH). 

To summarize: if the system is initialized to a single magnon state, with the spin-flip 
localized on the sample, the concurrence of the scattered neutrons is an oscillating 
function of time, which ranges from zero to a maximum of 0.77. In zero applied 
field, the oscillation is sinusoidal, and its amplitude falls as the spin density of the 
sample increases. The same result applies if the field lies outside a certain interval 
[Bz-, Bz+]. Conversely, if the field lies within this interval, it is always possible to 
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find an interaction time for wliicli the concurrence attains its optimal value. 




2 4 6 8 10 0.0 0.5 1.0 1.5 2.0 

r B, 

(C) 



1.0 P — ' — ' — ' — ' — • — ' — ' — ' — ' — • — ' — ' — ' — ' — • — ' — ' — ' — ' — ^ 

0.8- 

0.6 

0.4 - 
0.2 - 

o.o l - 

0.0 0.5 1.0 1.5 2.0 

Figure 5.4: Input state \iIjq), Q = 0, N = 1Q, X~1. (a) The evolution of the concurrence in 
time for i?^ = (sohd blue line) and = I (dashed green line), showing the development of 
a double-peaked structure when B^^ < B, < Bz+. (b) The discrepancy between the optimal 
interaction times, and the range of the interval [B^-, Bz+]. (c) The peak concurrence outside this 
interval (blue and green dot-dashed lines) and inside this interval (orange dashed line) as a function 
of the applied field. The solid red line represents the peak concurrence at Bz = 0, and shows the 
existence of the limiting field B^ = 2A (l — ;^). All quantities are expressed in natural units. 



5.1.3 Concurrence in the Thermodynamic Limit 

The expressions derived so far to describe the behaviour of the concurrence in dif- 
ferent parameter regimes are valid at all A^. In a macroscopic sample, may be of 
order 10^'^. Consequently, it is useful to note the functional form of these expressions, 
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and other relevant quantities, as becomes large. We find 

(j)^ B,~\ (5.46) 
ip^2{B,-X) (5.47) 

T^{N,X,B,)^^^ (5.49) 

C(iV,A,0,r)-^^^ (5.50) 

B: ^ A (5.51) 
B,± A (5.52) 

The only form to remain unchanged is that of C{N, A, BI,t)] hence equation fl5.36p 
describes the time-evolution of the concurrence in the thermodynamic limit also. 



5.2 Scattering from a fully Polarized Sample 

The derivation of the system wavefunctions for initial state is rather elaborate, 
because it involves the zero-, one- and two-spin-flip sectors of the Hilbert space. Let 
us write \iPq) explicitly in terms of the basis states: 

1^^) = a^ll) + a/3|4) + a(3\Q) + f3^\8). (5.54) 

There are now six eigenstates of T-C[ that have non-zero overlap with {iPq): 





= 


(5.55) 




= |6), 


(5.56) 


|3) 


= c|7)+c/|8). 


(5.57) 


|4) 


= c\2)+d\4), 


(5.58) 


|5) 


= d\7)-c\8), 


(5.59) 


|6) 


= d\2)-c\A), 


(5.60) 
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where 0, c and d are given by equations fl5.14p . fl5.16p and fl5.17p . respectively. The 
state of the system after the first neutron has scattered is 

+ a(3{d^e"^ + c^e'y^) |4) + a/5|6) 

+ p^cdie'''^ - e'y^)\7) + /^^(d^e-- + c^e^^^lS)], (5.61) 

with 



X = B, + X[l + ^] + (5.63) 



y = B, + A(l + i-)-0. (5.64) 



A = X + B,N, (5.62) 

The second period of free evolution simply introduces extra phases between the 
components of 1"?^^): 

+ e-''^^"/a/3(d2e»^ + c^e'^^K) + e-^'^'^^'f a(3\Q) 

+ p^cdie'""^ -e'y^)\7) + e-^'^'^'ff3\d^e'''^ + c^e'y^)\8)]. (5.65) 

Hence, the final state of the system is found to be 

+ cdapie"^ - e*^")(c2e"" + d^e'^^ + e-^'^^^f )\2) 
+ cdfgp^e'''^ - e'y^){e''^ - e™0|3) 

+ e-''^^"/a/5(d'e^"" + c2e'^0|4) 

+ e-^'^^^'fcdp\d^e'^^ + c2e'^")(e"" - e'^")|5) 

+ apc^d^e'""^ - e'y^f + a(3e-^'^^^'f {d^e'''^ + c^e^^") |6) 

+ cd(3\e'''^ - e'y^){g^e''^ + fe''"^)\7) 

+ e-''^^"/732(ci2e-- + c2e^2'-)2|8)], (5.66) 
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with 



3 \ , / 1 ^ ^ 



7 = ^^^{X,N,B,) = \lB^,-2B^i^l--jX+i^l + -j X\ (5.67) 
^ = ^'^{X,N,B,) = x(^^ -1^ + B, + j, (5.68) 

i:}N 

f = (5 59) 

V8A2(iV-l) + iV2i92' 

2Xy/2{N-l) 



(5.70) 

w = 3£, + A ( 1 + 4 ) - 7, (5.71) 



35, + A ( 1 + 4 ) +7- (5.72) 



v/8A2(iV- 1) + Ar2^2 

The form of 1^/;^^) is unfortunately too complex to yield a manageable analytical 
expression for the concurrence, so much of the analysis must be carried out using 
numerical simulation. Fortunately, the results of the previous section provide a 
useful starting point and somewhat reduce the magnitude of the task. 



5.2.1 Zero-Field Evolution 

The main features of zero-field time evolution of the concurrence are illustrated in 
figure [531 We note that: 

1. The concurrence does not depend on the duration of the period of free evolution 
between scattering events. 

2. The concurrence is an oscillating function of interaction time. 

3. The period of the oscillation is defined by the energy splitting of the single- 
excitation eigenstates that correspond to the spin fiip being shared between 
the interacting neutron and the sample. 

4. The amplitude of the oscillation decays as the density of the sample increases. 
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Therefore, the behaviour of the system mirrors that observed for input state \iPq), 
except the peak concurrence for input state \iPq) is considerably lower for all and 
seems to decay more slowly as the number of spins in the sample increases. 
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Figure 5.5: Input state IV'cf), Q = 0, = 0, a = /3 = A = 1. (a) The evolution of the 
concurrence as a function of interaction time for = 4 and iV = 10 (blue and green curves, 
respectively), (b) The concurrence as a function of the time between scattering events for N = 10. 
(c) The peak concurrence as a function of (blue dots). To ease comparison, the same curve for 
\iPq) is also shown (green dots). All quantities arc expressed in natural units. 



The only outstanding question is then the relationship between the concurrence and 
the neutron polarization. Figure 15.61 shows the peak concurrence as a function of a 
for different values of A^, and clearly illustrates that, for each A^, the concurrence 
is non-zero for a specific range of a. The upper boundary of this range is fixed at 

= 1; the lower boundary a__ increases with between approximately 0.4 and 
converging to at large A^. The value of a that maximizes the peak concurrence 
varies with A^, and tends to a ~ ^ in the large A^ limit. The peak concurrence 
decreases as a tends to unity, finally returning to zero at a = 1 (figs. 15.61 15.71 and 
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Figure 5.6: Input state \iPq), Q = 0, Br = 0, A = 1. The peak concurrence as a function of a for 
(a) iV = 4. (b) N=:10, and (c) N=30. All quantities are expressed in natural units. 
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Figure 5.7: Input state \iPq), Q = Q, B, ^ 0, N ~ 8, X ~ 1. (a) The peak concurrence as 
a function of a, showing a-{N 8) ~ 0.6. (b) The concurrence as a function of interaction 
time for a =^ 0.6. (c) The concurrence as a function of interaction time for a ~ 0.62. Note the 
non-analytical behaviour of the curve at a ~ q:_. All quantities are expressed in natural units. 
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Figure 5.8: Input state |V'(f), Q = 0, = 0, N ~ 9, X = 1. (a) The peak eoncurrence as a 
function of a, showing a-{N = 9) « 0.62. (b) The concurrence as a function of interaction time 
for a = Q_. Note the concurrence is no longer discontinuous at this value of a. (c) The concurrence 
as a function of interaction time for = (solid blue curve), a = 0.65 (dashed red curve) and 
a = ^ (dot-dashed green curve). The optimal neutron polarization is the latter. All quantities 
are expressed in natural units. 



The absence of entanglement when a = 1 is evidently due to the fact that when the 
neutron spin is parallel to the positive 2;-axis the input state is an eigenstate of the 
interaction Hamiltonian. To explain why similar behaviour is observed at a = 0, it 
is useful to examine the state of the system at the different stages of the protocol. 
The first scattered state can be written as \ip2) — r|7) + A|8), and for the protocol 
to succeed this state must be entangled. The concurrence of the first neutron and 
the sample can be expressed in closed form as 



VN\ sin [A (1 + ^)t]\ JN^ + 1 + 2N cos [2A (l + r] 
C. = 2\rA\=4 , (5.73) 

which is once again a regularly oscillating function of r with period T^. For > 5, 
Ci peaks at time t = with a maximum value A\/N {N — 1){N + 1)~^, while for 
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N < 5 the peak concurrence reaches unity. We can verify this is true by showing 
there exist times at which the square moduh of the coefficients of state \ip2) ^ire 
equal to ^, imposing for instance 



(5.74) 



Substituting the values of c and d and solving for r then yields 



1 N 



T 



AA^+ 1 



sm 



N+l \ 



(5.75) 



which has real solutions for < 5. Hence, the concurrence of the first neutron and 
the sample does not reach unity for N > 5, but can assume non-zero values for all 
N. 

Evidently, then, it is not the first stage of the protocol that is at fault. Let us 
therefore consider the state of the system after the second scattering event; this 
takes the form \tpf) = 0|3) + T|5) + H|7) + Q\8). The concurrence of the neutrons 
can be expressed analytically as C = 2 (|TH| — |0f2|), which substituting the zero- 
field values of the coefficients takes the rather complicated form 



C 



SGjN) sin^ [A (1 + jj) t] 
(iV + l)4 



NF{N) - ^/2N{N -1)G{N) 



(5.76) 



with 



G{N) 




'a^2 _2Ar + 5 + 4(A^- l)cos 



2A(l + -,r 



(5.77) 
(5.78) 



The pre-factor to equation (15.761) is always positive. In the limit of large A^, the 
function enclosed in square brackets tends to A^^ (l — v^) , which is always negative. 
Therefore, the concurrence between the neutrons must always be null. 

In summary, if the applied field is zero the concurrence of the neutrons is deter- 
mined by the spin density of the sample and the neutron polarization. The peak 
concurrence is maximal when the neutrons are polarized at approximately 60° to 
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the quantization axis, but this maximum never exceeds 0.16 and falls roughly as 
5.2.2 Finite Field Evolution 

Let us now consider the effect of applying a finite field. Bearing in mind the results 
of sections 15.1.21 and 15.2.11 there are several features we might expect to see: 

1 . Independence of the concurrence from rj- ; 

2. Discontinuities and null concurrence when a < 

3. Increase or decrease of the peak concurrence according to the applied field and 
the neutron polarization; 

4. Maximal concurrence when = B*; 

The first point holds for all and B^, hence the concurrence is once again inde- 
pendent of the time between scattering events (see fig. 15.91) . It is also undoubtedly 
true that if a = 1 the neutrons must remain unentangled, because the initial state 
of the system is an eigenstate of the interaction potential. In addition, we know 
the final point must apply, owing to to structure of the potential. The remaining 
items are easily verified from figures 15.101 and 15.111 which illustrate the evolution of 
the concurrence in time for different field strengths and polarizations. We observe 
that, for fixed a, the peak concurrence Cp improves when a field is applied, provided 
the strength of the field does not exceed a rough upper limit of 3A. Surprisingly, 
however, the concurrence in finite field is now non-zero for all a 7^ 1. 
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Figure 5.9: Input state Q = 0, iV = 10, A = 1. The evolution of the concurrence as a 

function of r and for different field strengths and neutron polarizations, (a) Bz ~ 0, a = 
(b) = 0.3, a = ^, (c) B, = 1, a = ^, (d) = 0, a = ^, (c) B, = 0.3, a = ^, (f) B, = 1, 
a = All quantities are expressed in natural units. 
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Figure 5.10: Input state |'0(f ), Q ~ 0, a ~ 0.7, iV = 10, A = 1. The evolution of the concurrence 
with interaction time at different Bz, showing the existence of a limiting field of approximately 3A. 
(a) Bz = 0, (b) Bz ~ 0.3, (c) Bz = 1, (d) Bz = 2, (e) Bz = 3.5. All quantities are expressed in 
natural units. 
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Figure 5.11: Input state |?/'(f )i Q = 0, A = 1, = j^, N ^ 10. The evolution of the concurrence 
with interaction time at different a for (a) a = 0, (b) a = 0.2, (c) a = (d) a = (e) 
a = and (f) a ~ 0.95. Note the concurrence is finite for all a ^ 1. All quantities are expressed 
in natural units. 



In fact, Cp is maximal when a = 0. To explain this behaviour, it is instructive to 
analyze the state of the system at each stage of the protocol. From equation fl5.66l) . 
neglecting the state of the second neutron, the first scattered state l?/;^) can be 
written as 

l^f ) = a|00) + Pcd (e"" - e'y^) |01) + /3 {d^e"^ + c^e'^^) |10). (5.79) 
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Evidently, for a given c and d the entanglement of this state is directly proportional 
to (3, and is maximal if (3 = 1. Therefore, to optimize the first stage of the protocol 
the neutron spin must be anti-aligned with the 2;-axis. Let us then assume (3 = 1 
and consider the state of the system after the second scattering event; this has the 
form \ipf) = 6'|3) + T'|5) +S'|7) + f2'|8). As before, the concurrence of the neutrons 
is C = 2 (|T'S'| — |6'f2'|), hence substituting the values of the coefficients at finite 
Bz yields 



C 



(4A2 + iV^2)3[8A2(iV 
Nip'^ sin^ (pT 



A2 (iV- 1) 



Ar2 



^/2eN{N^^ I (/j'i? sin 0r sin 7r | 

(5.80) 



A2 



-32 (iV- I)i92__sin2 7r 



with 



16A^ 4 8A2 

^ + cos2</.r. 



(5.81) 



Both terms in the square brackets are always positive. The concurrence will therefore 
be finite if there exist values of Bz and r for which the magnitude of the second 
term is greater than that of the first. Let us set a lower limit for the concurrence 
of C'min = 10~^. For any arbitrarily small value B^ < -B*, C exceeds Cmin at 
certain values of r satisfying BzT > O.IA. Not all values of r described by this 
condition yield a finite concurrence, because C is not a smooth function of time (fig. 
I5.12p . Equation fl5.8Up peaks when the difference between the bracketed terms is a 
maximum, which for B^ = B* first occurs at time 



TlN 



T 



2XV2N 



(5.82) 



Hence, as for initial state iV'o^), setting r = r* the peak concurrence converges to 
0.77 in the limit of large N (fig. Km . 
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Figure 5.12: Input state Q = 0, a = 0, = 10. The development of a non-zero concurrence 

as B,T exceeds the hmiting value O.IA. (a) B,_ = lO^^A, A = 1, (b) = O.OIA, A = 2, (c) 
Bz — 0.003A, A = 3. The insets show that, on a sufficiently small timescale, the 'spikes' pictured 
have finite width. All quantities are expressed in natural units. 
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Figure 5.13: Input state Q ^ 0, a = 0, B^ ~ B*, \ = 1. The peak concurrence as a function 

of N (red dashed line), which at large N converges to 0.77 (blue solid line). All quantities are 
expressed in natural units. 
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We conclude that, from a purely mathematical viewpoint, the performance of the 
system is optimal when the neutrons are polarized in the negative ^-direction. How- 
ever, in practical terms this may not be the most advantageous choice, because when 
the field is weak the concurrence at short times is poor. Comparing the evolution 
of the concurrence for different field strengths at a = and a = ^ illustrates that, 
for BzT < 2A, the optimal spin orientation is in fact the latter (fig. I5.14p . 

(a) 

,1 I - 0.4 - I - 
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Figure 5.14: Input state Q = 0, N = 10, X ^ 1. The evolution of the concurrence at short 

times for a = (blue curve) and a = (green curve) at (a) Bz = 0.01, (b) Bz = 0.1, and (c) 
Bz — 0.2. All quantities are expressed in natural units. 



Figure 5.14(c) shows evidence of a rather interesting beating effect. To explore this 
further, let us examine the concurrence for generic a and B^; this is illustrated 
in figure ISTTSl If the field lies within an interval [Bz-,Bz+] ~ [fi^A (l + j^)], no 
evident systematic behaviour is observed. Conversely, if Bz is outside this interval 
there emerges a pattern which resembles a high frequency oscillation modulated by 
a low frequency envelope. The period of the rapidly oscillating component is the 
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familiar quantity T^. The period of the envelope function is defined by 

T, = (5.83) 
10-71 

where 7 = 7^^(A^, A, -B^). The behaviour of the system is therefore determined 
by the beating of four eigenstates, of which two correspond to the spin-flip being 
shared between the first neutron and the sample, and two correspond to the spin-fiip 
being shared between the second neutron and the sample. The optimal interaction 
time can then be estimated by tracking the evolution of the peaks of the envelope 
function. 
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Figure 5.15: Input state \i/jq), Q — 0, a — N = 10, A = 1. The emergence of a beating effect 
in the evolution of the concurrence as a function of interaction time. In figures (a) to (f), the 
applied field sweeps through the range < B^- to B^ > Bz+, with (a) Bz ~ 0.1, (b) B^ = 0.3, 
(c) Bz = 0.5, (d) Bz = I, (e) Bz = 2, (f) Bz = 3. All quantities are expressed in natural imits. 



In summary, when a field is applied the peak concurrence always increases, provided 
the field does not exceed an upper limit. This limit depends on the neutron polar- 
ization, but is generally of order 4A. For a given field, Cp is a maximum when the 
neutron spin is anti-aligned with the quantization axis. If the field lies outside the 
interval [Bz^,Bz+] this maximum occurs at approximately r = However, if 

the field is weak then is large, hence at short times the system performs better 
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if both a and (3 are non-zero. If the field hes within the interval [Bz-,Bz+] the 
behaviour of the concurrence is more difficult to predict, but setting Bz = B* and 
a = the concurrence takes on a well-defined analytical form and peaks at pre- 
dictable times. For N > 5 the peak value exceeds 0.7, and converges to 0.77 as 
becomes large. The performance of the system at zero momentum transfer for 
input state therefore matches in efficiency that observed for input state I'iPq), 
although its behaviour is overall more complex and less predictable. 



5.3 Concurrence and Negativity 



One of the most interesting properties of the neutron concurrence is the invariance 
with respect to the second period of free evolution. In zero field, this arises from the 
form of the free Hamiltonian Hq. In finite field, however, Hq is no longer a multiple of 
the identity, nor does it commute with Tii or 7-^2; therefore the explanation is not as 
evident. To gain insight on the origin of this behaviour, it is useful to investigate its 
dependence on the structure of the interaction potential. Specifically, let us focus on 
the second scattering event. One finds that, provided the interaction Hamiltonian is 
hermitian and displays an excitation number conserving symmetry, the concurrence 
is always independent of the second period of free evolution. It would therefore 
seem that the sample behaves as an 'entanglement safety-deposit box', which holds 
the spin information deposited by the first scattering event until a second neutron 
bearing the right 'key' comes along to claim it. In other words, the entanglement 
between the neutrons is conditional on the second scattering event being able to 
'undo' -in part or in full- the transformation on the sample spin precipitated by the 
the arrival of the first neutron. 

It is interesting to ask whether the invariance with rj- extends to all entanglement 
measures. Figure 15.161 illustrates the evolution of the logarithmic negativity during 
the second period of free evolution for input state and reveals a definite vari- 

ation with t'j-. The discrepancy between the concurrence and the negativity may 
arise from the difference in the physical significance of the two quantities: given an 
arbitrarily large number of maximally entangled states, the concurrence quantifies 
the entanglement cost of generating an arbitrary state p under local operations and 
classical communication, whereas the negativity quantifies the entanglement cost of 
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generating an arbitrary state p under positive-partial-transpose-preserving opera- 
tions. The dependence of the negativity on rj- then becomes a useful indicator of 
the inequivalence of the two classes of operations. 
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Figure 5.16: Input state Q ^ 0, = 0.35, = 10, q = A = 1. Figure (a) shows the 

concurrence of the neutrons remains constant during the second period of free evolution. Figure 
(b), conversely, shows the negativity of the neutron state varies with rj. Figure (c) shows the 
spread of negativity values for states with a given concurrence, sampled at random values of r. 
The red and green points are maximum and minimum negativity values, respectively. The solid 
and dashed lines are the bounds for the maximum and minimum negativity of a state with a given 
concurrence |115| . All quantities are expressed in natural units. 

5.4 Parallels with Optical Entanglement Genera- 
tion 

The entanglement of massive particles via successive interactions with a macroscopic 
mediator is by no means a novel idea. Almost ten years ago, a group led by Haroche 
of the Ecole Normale Superieure (ENS) reported the generation of EPR pairs of 
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atoms via the exchange of a single photon in a high- (5 cavity according to the 
following scheme |114] . A single atom is prepared in an excited state |e) and is 
made to interact with a cavity in its ground state. Assuming the cavity and the 
atom exchange photons at a rate Q, after a period of time Ti = ^ the combined 
state of the atom and the cavity becomes 

\^l;,) = ±{\eO) + \gl)), (5.84) 

where the right-most number refers to the photon number inside the cavity. After a 
delay time T2, a second atom, this time prepared in its ground state \g), is made to 
interact with the cavity for a period T3, which is twice the duration of ri. If the first 
atom has emitted a photon, the second atom becomes excited; otherwise it remains 
in l^f). As a result, the composite state of the system is 

|^2) = ^(|^7eO) + |e(70)), (5.85) 

which is a maximally entangled state of the atoms. 

It is shown in figure 15.171 that if the ENS protocol were exactly mapped to the 
scattering scheme I have described, it would be possible to generate fully entangled 
neutron states. Indeed, if we restrict ourselves to the single excitation subspace by 
working with input state I'iPq), near unit concurrence is achieved when = 2ti 
[figure 5.17(a)| . However, such a mapping requires the possibility of separately 



tuning the interaction times of consecutive neutrons, which is not currently possible. 
Therefore, to compare the ENS protocol and the scattering scheme on a more even 
footing, it is useful to examine the evolution of state |eeO), which is contextually 
equivalent to state |8) = |110) of the Bloch basis, for equal interaction times ri = t^. 
In this case, the concurrence of the atoms is given by 

Ch = 2 ^cos^ Ar| sin V2Arsin Ar| - sin^ Ar| cos V2Ar cos Ar| j . (5.86) 

This expression peaks at a value 0.77 for r ~ and, as we will see, greatly 
resembles the evolution of the concurrence in optimal field for input state |8), when 
scattering takes place from an anisotropic sample [see equation (15.1181) and section 
15. 7j . The ENS scheme and the scattering protocol are therefore comparable in all 
but the timescale. 
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Figure 5.17: The time-evolution of the neutron concurrence assuming independently tuneable 
interaction times, for TV = 10, = B* and A = 1. Figure (a) relates to input state IV'o^), for 
which optimal concurrence is first achieved at T3 ~ 2ti. Figure (b) relates to input state |^/'(f), 
with a = and (3 = \. All quantities are expressed in natural units. 



5.5 Experimental Feasibility 



From the results presented so far in this chapter, it would seem that the scattering 
protocol of chapter [3] offers a valid means of creating entanglement between distinct 
neutrons. The question to ask is therefore: can it be done? Let us review the 
quantitative requirements of the protocol: 



1. The spin lattice relaxation time (Ti) must be longer than the first period of 
free evolution, which roughly coincides with the time taken by the neutrons 
to reach the sample; 

2. The phase coherence time (T2) must be longer than the period of free evolution 
between scattering events; 

3. The neutron coherence volume must be comparable to the size of the sample; 

4. The optimal field B* and the optimal interaction time r* must be attainable 
values; 

5. The protocol must be robust against experimental uncertainties; 
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6. The final state of tlie neutrons must yield a measurable witness. 

To address these points, I take as a reference the technical specifications of the PF2 
source of UCNs at the Institut Laue-Langevin (ILL) in Grenoble (see Appendix |B] 
and [TUU] l 

For UCNs with velocity v = 7 ms~^ and a flight path of 10^^ — 1 m, the time taken 
to reach the sample is of order 10~^ — Is. Ti times consistent with these figures 
are achievable in materials such as phosphorus-doped silicon or NtQiCeo |116l 1117] . 
Furthermore, for a neutron flux F = 10^ m~^s~^ and a sample area of, say, 10~^ m^, 
the time between scatterings might be of order 10~^ s. T2 times similar to this are 
also attainable in these systems at low temperature |116j . 

Next, we require the neutron coherence volume to be comparable to the size of the 
sample. This condition is expressed by the uncertainty relation of equation (13.61) . 
which in terms of the neutron velocity becomes 

Av = — (5.87) 

For a sample of length AL = 10 cm, Av is of order 10~^ ms~^. For neutrons traveling 
at 7 ms~^, therefore 

— ^ 10~l (5.88) 

V 

Currently, this is somewhat unrealistic. Resolutions of order 10~^ have been achieved 
in the field of neutron spin-echo spectroscopy |118] , but this yields coherence lengths 
of only AL ^ 10^^ — 10^^ m. The situation could be improved if the neutron velocity 
were reduced, which might be achieved by prolonging or improving the efficiency 
of the cooling stage described in Appendix [HI or by allowing the neutrons to slow 
down under the effect of gravity (see below). 



For the remainder of this discussion, I will consider scattering from a sample in a 
single magnon state. To estimate working values of B* and r*, it is useful to re- write 
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equations (15.351) and fl5.37p in SI units and take the large limit: 



dNfJ'NfJ'O 



(5.89) 



T* = 



hVn 

4A 



^QNl^'NlJ'OgelJ'B 



(5.90) 



Hence 



Bl ^ 10 




-32 (^y\ 



(5.91) 



T 



10 



(5.92) 



Assuming a lattice constant Oq = 10~ m and a macroscopic number of scatterers 
]\f = 10^^ gives an optimal field -B* = 10^^ T and a sample volume = 10~''m^. 
For these parameters, the neutron velocity is v ~ 10~^ms~^, corresponding to an 
interaction time r* ~ 100s. Currently, such low values of v are not attainable. This 
is not a fundamental constraint, but a technical limitation. In principle, UCN could 
be slowed to this energy regime by causing them to decelerate in a gravitational po- 
tential |119j : a UCN with v = 7ms~^ would stop at a height of approximately 2.5m, 
having converted all its kinetic energy to gravitational potential energy. Unfortu- 
nately, there are several problems with this approach. First, the neutron absorption 
cross-section is inversely proportional to the neutron velocity. Therefore, the slower 
the neutrons, the greater the chance they will be absorbed by the scatterer or indeed 
by any other component of the experimental setup. Second, even if it were possible 
to produce such a slow neutron beam, one would expect its intensity to be extremely 
low. Therefore, it is difficult to see how one could produce enough neutrons in the 
required energy range to carry out a statistically significant measurement of the 
witness (see discussion below on how to measure the witness). 

It was stated in chapter O that the protocol would be deemed successful if the 
neutrons were found to share an amount S of entanglement, corresponding to a 
concurrence Cs. I now set = |, which implies we still consider the protocol 
successful if the concurrence deviates from its optimal value by 0.1. For simplicity, 
the possibility of simultaneous errors in the field and in the timing will be excluded. 
In other words, I assume we operate either at optimal field, or at optimal time. 
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At optimal field, the time evolution of the concurrence is described by equation 
f l5.36p . An acceptable timing error can then be estimated by measuring the half- 
width of this function in the vicinity of r*, at C = C^-. For the set of parameters 
defined above, one finds r* ~ 150 s and Ar ^ ±70 s [see fig. 5.18(a)| , hence 



At 

T* 



Av 

V 



0.5. 



(5.93) 



Note this value does not take into account the restrictions on the neutron coherence 
volume. 



Let us now fix r = r*, and consider the effect of errors in field calibration. The 
concurrence as a function of applied field is given by equation fl5.27p . As before, the 
allowed spread in can be estimated from the half-width of the curve at C = C^: 
around Bz = B*. From figure 5.18(b) 



ABz 
AB^ 
B* 



10 



10 



-13 



-11 



(5.94) 
(5.95) 



To recap, for a completely successful realization of the protocol we require 

^ « 10-. 

V 

AB. 



Bt 



11 



(5.96) 
(5.97) 



The second condition may be attainable, but the first is currently not |120lll211ll22] . 
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Figure 5.18: Input state Q = 0. (a) The spread in concurrence around t = t* at B, ^ B*. 

The origin of the ordinate axis is set at Cf = 0.67 to better illustrate the tolerance interval At. 
(b) The spread in concurrence around B^ ~ B* at r = t*. 



The final point regards the measurement of the witness. We saw in section 15.1.21 
that a successful realization of the protocol yields the scattered state of equa- 
tion (15.401) . As expected, the partial transpose of is negative. The eigenvector 
corresponding to the negative eigenvalue has the form 

|e-) = a|00) +6|11), (5.98) 

where both a and h are approximately equal to It has been shown that witness 
operators constructed from states of this form can be decomposed as follows [67] : 

W = a^\z+z^){z-^z-^\+b'^\z-z-){z-z-\ (5.99) 

+ ab (^\x~^x~^) {x~^x~^\ + \x~ x~) {x~ x~ \ — \y~^y~) {y~^y~\ — \y~y~^) {y~y^\) , 

where \x^), \y^) and 1^^) are the spin-up and down eigenstates of the Pauli matri- 
ces ax, (Ty and a^, respectively. Such a witness could be measured with as few as 
three device settings, provided one could detect all outgoing neutrons and measure 
each component of their spin. A Stern-Gerlach experiment in an arbitrary direction 
would achieve this objective. Bearing in mind the stochasticity of existing neutron 
sources, (W) could then be evaluated on pairs of neutrons detected in rapid suc- 
cession following a gap long enough to reset the sample state. This is conceptually 
possible, though very challenging from a technical viewpoint. 
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5.6 Scattering from an Anisotropic Medium in a 
Single Magnon State 



We have assumed so far that the sample is an isotropic medium. This is reflected in 
the form of the interaction Hamiltonian of equation (15.11) . which is proportional to an 
electronic gf-factor with equal components Qe along all three axes. In general terms, 
the (^-factor is a tensor with diagonal elements Qx, Qy and Qz, whose magnitude can 
vary in space according to the chemical nature of the environment, or the presence of 
external fields. This property is of great relevance in the field of spintronics, which 
exploits g-tensoT anisotropies to manipulate electron spins using electrical signals 
[123] . It is then interesting to explore how the scattering scheme might perform 
if the sample were anisotropic. To describe this situation accurately, one would 
have to account for the possible breakdown of the LS model. However, as a very 
basic approximation, we can assume LS coupling still holds, and simply choose a 
(yf-tensor with unequal components along the three spatial directions. Let us choose 
an extreme example: Qx = 9y = Qe^ cind Qz = 0. In this limit, one might imagine a 
spin potential of the form 



As before, this can be expressed in terms of an 'effective' interaction potential similar 
to equation (15.21) . provided one takes the correct form for Vm- Quantitatively, V^*""" 
differs from V'^"'^ because it lacks a diagonal component; hence the diagonal elements 
of H'„^ are simply those of Hq. As a result, the equations describing the states of 
the system in the different stages of the protocol are identical to (15.191) . (15.201) and 




(5.100) 
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(I5.24p . provided one sets 



4A2 

= ^^y(N,\,B,) = ^B^, + —, (5.101) 
A = B,{N-1)-(I), (5.102) 

(5.103) 

(5.104) 
(5.105) 




The concurrence between the neutrons is then 



CiN, A, B^, r) = ^-^^^^\jBl + ^ cos^ 0r. (5.106) 

There are striking similarities between this and equation (15.271) . First, fl5.106p is an 
oscillating function of the interaction time, and independent of the time between 
scattering events t'j. The period of the oscillation has the form of equation fl5.28p . 
with the value of (f) given by equation fl5.10ip . Second, (15.1060 reduces to the form of 
(I5.36P if the field is set to zero, suggesting the concurrence is maximal when there is 
no applied field. Indeed, as discussed in section 15.1.21 at = the single excitation 
eigenstates of Ti^ are reduced to equally weighted superpositions. 

In zero field, the evolution of the concurrence in time shows the familiar double- 
peaked structure observed in figure 15.21 This structure disappears when the applied 
field exceeds a certain threshold, to be replaced with a sinusoidal oscillation [fig. 
5.19(a)| . To determine the threshold field we repeat the analysis of section [5.1.21 



maximizing equation (15.1060 with respect to r and imposing equality with half the 
oscillation period [fig. 5.19(b)| . One finds 



1 _i 

"^B = TTT COS 



iV52 + A^ 
3A2 



(5.107) 



which reduces to the form of (I5.37P when Bz = 0. From equations (I5.10ip . (I5.107P 
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and fl5.28p . the threshold field Bt is then 



B, = (5.108) 



Substituting into equation fl5.106p . the peak concurrence above and below Bz = Bt 
can therefore be expressed by the following: 

C+(iV,A,i?.) = (5.109) 
C_(7V,A,i?,) = (5.110) 



This is illustrated in figure 5.19(c) The trends described by fl5.109p and fl5.110p 
are qualitatively similar to those observed for an isotropic sample. Apart from the 
obvious equivalence of Ci and C_ at optimal field, the behaviour of C+ strongly 
resembles that of Co for Bz > Bz+, and both tend to the same limit as i?^ — oo 
[fig. [519(d)] . 



To summarize: if the system is initialized to a single magnon state, with the spin-fiip 
localized on the sample, the concurrence of the scattered neutrons is an oscillating 
function of time, which ranges from zero to a maximum of 0.77. Above a threshold 
field Bt, the oscillation is sinusoidal, and its amplitude falls as with N. Below this 
threshold, the concurrence can always be boosted to the peak value. The only salient 
difference between scattering from an isotropic or anisotropic sample therefore lies 
in the value of the optimal field, which is finite in the former case, but zero in the 
latter. 
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Figure 5.19: Input state \^q), Q = 0, N=10, X=1. (a) The evolution of the concurrence as 
a function of interaction time for Bz ~ 0, ~ 0.1 and Bz — 1 (solid blue, dashed red and dot- 
dashed green curves, respectively), (b) The discrepancy between the optimal interaction times as 
a function of applied field, showing the existence of a threshold field Bt. (c) The peak concurrence 
as a function of Bz, for Bz < Bt and Bz > Bt (red dashed and blue solid lines, respectively), (d) 
The convergence of C+ (blue curve) and Co (red curve) in the limit of large N. All quantities are 
expressed in natural units. 



5.7 Scattering from an Anisotropic Medium in a 
fully Polarized State 

The derivation of the system states for initial state is identical to that outlined 
in the previous section, hence to avoid unnecessary repetition I will simply quote 
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the results of the calculation. The first scattered state can be written as 

l^f) = e-'^^[a\-'y'\l) + apcd{l-e-^"^^)\2) 

+ (3^cd{l - e-^''^^)\7) + (3\d^ + c^e-^''^^)\8)], (5.111) 

where = 0^^(A, A^, B^), A = Bz{N — 1) — 0, and c, d and y are defined by equations 
(15.1031) . (15.1041) and (I5.105p . respectively. The second period of free evolution simply 
adds phases: 

\^^) = e-'^^e-'''^^^'f[e-'y^a^\l) + «/?crfe''^^"/(l - e^2'^")|2) 
+ apid^ + c2e-2^<^") |4) + e-'y^af]\6) 

+ (3^e^'''^^fcd{l - e~^'^^)\7) + (3\d^ + c^e-^'^^)\8)]. (5.112) 
The final state of the system is therefore 

+ cdape^'^^^'f{l - e-2'^")(c2 + d^e''"''^^ + e-^'^^^'fe-'y^)\2) 

+ cdfgp^e^'^^^'fil - e-2^^<^-)(e-™" - e-'"")|3) 

+ e-'^"a/?(d2 + cV^")|4) 

+ cd(3\d^ + c^e-^''f'^){l - e-^'^^)\5) 

+ cdf3^e^''''^'f{l - e"2^<^0(/'e-'"" + ^'e"™")!?) 
+ (3\d^ + c^e-^'^^f\8)l 
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with 




B'. + ^{N-1), (5.113) 

(5.114) 

(5.115) 

(5.116) 
(5.117) 

The resuhs of section show the behaviour of the system in the presence of sample 
anisotropy is similar to that observed for an isotropic sample, offset by a field B*; this 
is true here also. The evolution of the concurrence in time at zero field can therefore 
be summarized in the following observations: (i) the concurrence is independent of 
rj.; (ii) for a 7^ 0, the concurrence is not a periodic function of the interaction time; 
(iii) the peak concurrence is finite for all a 7^ 1, but maximal when a = 0, in which 
case it takes the form 



C = 4 



4Ar / 2Ar . 2X^/2{N - 1)t . 2Ar 2\^2{N - 1)t 



sm — = cos —;= sm — sm — = cos 



N \ VN N ViV N 

(5.118) 



This peaks above 0.76 for all N, converging to a value 0.77 as — * 00. These 
characteristics are illustrated in figures 15.201 and 15.211 
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Figure 5.21: Input state |?/'(f ), Q ^ 0, Bz = 0, N = 10, X = 1. The evolution of the concurrence 
with interaction time at different a for (a) a 0, (b) a = 0.2, (c) a = (d) a = (e) 
a = and (f) a — 0.95. Note the concurrence is finite for aU a ^ 1. All quantities are expressed 
in natural units. 



The behaviour of the concurrence in finite field is more complex. For all a, and 
Bz the concurrence is an oscillating function of time. At Bz < A, the oscillation 
is irregular and lacks an easily discernible periodicity, but the peak concurrence 
remains roughly constant with Bz- Conversely, for Bz ^ X, there emerges a beating 
effect similar to that observed in figure 15.151 for Bz > Bz+ and Bz < Bz^ (see figure 
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I5.22p . The timescales involved are once again (rapid oscillation) and (envelope 
function), with = (p^^ and 7 = 'y^^. In this range of Bz, the concurrence drops 
rapidly with increasing field (figure I5.23P . Bearing in mind the form of c and d, 
this is due to the fact that, as the field is raised, the single-excitation eigenstates of 
the potential deviate increasingly from being equal superpositions. Indeed, by the 
same argument the concurrence at fixed increases with increasing A^. For all B^, 
the peak concurrence is strongly related to the neutron polarization, and attains a 
maximum when a = 0. The concurrence in this case takes the form 



In summary, with the exception of a few quantitative details, the performance of 
the protocol is not significantly affected by a strong sample anisotropy. As observed 
previously, if the sample is in a fully polarized state the behaviour of the system 
becomes more complex and less predictable, but no less effective. 



C 




(5.119) 
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Figure 5.22: Input state I'tpo), Q = 0, a = = 10, A = 1. The evolution of the concurrence 

with interaction time at different B,_ for (a) B,_ = 0, (b) B, = 0.2, (c) B, = 0.3, (d) B^ = 0.5, (e) 
Bz = I, and (1) B^ = 2. Note that, for < A, the peak concurrence remains roughly constant 
with Bz- All quantities are expressed in natural units. 
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Figure 5.23: Input state Itpo), Q = 0, Br = 0.5, iV = 10, A = 1. The evolution of the concurrence 
with interaction time at different for (a) = 0, (b) = 0.2, (c) B^ = 0.3, (d) B^ = 0.5, (e) 
Bz = 1, and (f) B^ = 2. Note that, for B^ < A, the peak concurrence remains roughly constant 
with Bz ■ All quantities arc expressed in natural units. 



5.8 Conclusions 



Two distinct neutrons can be entangled with respect to their spin degree of freedom 
if they are made to scatter sequentially from a macroscopic sample having the prop- 
erties of a ferromagnetic insulator. I have considered forward scattering from an 
isotropic or strongly anisotropic sample. In both cases, the system must be initial- 
ized to a well-defined state, chosen such that the sample contains at most one spin 
excitation. More freedom is permitted in the choice of the neutron states, however 
there does exist an optimal polarization, which is determined by the initial state 
of the sample and the magnitude of the external magnetic field. It must also be 
possible to tune the time each neutron spends interacting with the sample. This 
can be done by regulating the neutron velocity. Provided these conditions can be 
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met, the neutrons can come to share as much as 0.68 ebits of entanglement. 

In hght of present-day achievements across the range of disciphnes employing the 
tools and methods I have discussed, a realization of the entanglement scheme in its 
present form does not seem feasible. However, future progress in neutron cooling 
and manipulation may change this state of affairs. 

Leaving aside for a moment the problem of producing neutrons with sufficiently large 
coherence volumes, the most difficult part of the protocol is the entanglement detec- 
tion stage. Evaluating the witness of equation f l5.99p involves measuring expectation 
values of the polarizations of the scattered neutrons along all three axes. Currently 
operating neutron spectrometers are capable of full XYZ polarization analysis |124] , 
but with two limitations. First, for a given scattering event, it is only possible to 
measure the spin of the outgoing neutrons along one axis - either X, Y or Z. Second, 
only neutrons with spin in one direction along the chosen axis are detected. For 
example, suppose one had a beam of scattered neutrons containing a mixture of 
particles with spin up or down along the z-direction, and wished to measure the 
^-component of the spin of each neutron. One could either set the analyzer to detect 
neutrons with spin up or to detect neutrons with spin down. If the setting were 
chosen as spin up, none of the neutrons with spin down would be detected. Instead, 
the relative occurrence of spin up and down in the beam would be inferred by relat- 
ing the measured flux to the total expected scattered flux. Now, for the purposes 
of measuring the witness of equation (15.991) . this is not a problem until one comes 
to the term in \y^y~){y~^y~\ + he, which can only be measured reliably if both the 
scattered neutrons are detected. It may then be necessary to employ more than one 
measurement strategy to evaluate the witness. One might, for example, measure 
the ffist four components using a spectrometer such as D7, currently in operation at 
the ILL |124j , and measure the final two by routing the scattered neutrons through 
a rotated Stern-Gerlach apparatus. This would certainly be extremely challenging, 
but is not a priori impossible. 



Chapter 6 



Many- Neutron Scattering 



In this chapter, I extend the protocol described in chapterlM to study the concurrence 
of an arbitrary pair of neutrons scattered from the sample. Once again, I consider 
sequential scattering at zero momentum transfer from a sample in one of two possi- 
ble states - the polarized, and single-magnon states of equations (13.131) and (13.151) . / 
find that, under optimal circumstances, pairs of neutrons can share a finite concur- 
rence even if they are not scattered in succession. It is therefore possible to create 
entangled states of many neutrons; in fact, for certain choices of input parameters, 
the amount of entanglement of any neutron with any other becomes a slowly-varying 
function of the neutron indices. Consequently, the probability of creating an entan- 
gled neutron pair is high, provided the sample can be periodically reset to its initial 
state. In addition, measurement of the witness is no longer subject to the constraint 
of detecting the first two scattered neutrons. 



6.1 Introduction 

I showed in the previous chapter that, conditional on some hmitations, the first two 
neutrons to scatter from a sample prepared in a specific eigenstate of its internal 
Hamiltonian can be measurably entangled. For an optimal realization of the proto- 
col, this entanglement is detected by the witness operator of equation (15.991) . whose 
expectation value is calculated on the states of the two outgoing neutrons. As the 
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neutron scattering cross-section is small, if the beam is sufficiently dilute and the 
sample state is periodically reset, it is reasonable to assume the ffist two neutrons 
detected in rapid succession following a reset pulse are indeed the ffist two neutrons 
to scatter from the sample. 

But what if this were not the case? Equation (15.991) represents a particularly simple 
decomposition of a generic witness operator capable of detecting the entanglement of 
quantum states having some overlap with the state p|01)-|-g|10). This decomposition 
is valid provided the eigenvector corresponding to the negative eigenvalue of the 
partially-transposed density matrix of the neutrons has the form a|00) -|- fe|ll). The 
state p"^ of equation f l5.40p falls into this category, however the state pm^n of two 
arbitrary scattered neutrons m and n may not. This means the witness in its 
present form is no longer guaranteed to detect the entanglement of Pm,n- In addition, 
detection of the wrong neutron pair could lead to an inaccurate assessment of the 
performance of the protocol. Suppose, for example, that pi,2 is entangled, but p2,3 is 
not. By accidentally detecting the second and third neutrons, rather than the ffist 
and the second, we might conclude the protocol has failed, whereas in fact it has 
succeeded. 

It may then be useful to know exactly how much rests on the ability to detect the 
first two neutrons. To this end, I study the concurrence Cm,n of an arbitrary neutron 
pair m and n, and calculate the average concurrence over all detected neutron pairs 
{m, n}. 

6.2 Scattering from a Sample in a Single Magnon 
State 

We saw in section ISTTl that when the system is initialized to \iPq), quantum state 
evolution is restricted to the single excitation sector of the Hilbert space. This 
allows the concurrence of the scattered neutrons, and the behaviour of the system 
in general, to be described in terms of analytical expressions, making clear how 
the performance of the protocol is affected by changing the input parameters. The 
same argument applies to the situation we now wish to describe. By deriving the 
state of the system after the first few scattering events, with the help of equation 
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fl3.32p . there quickly emerges a pattern in the coefficients of the scattered state. This 
pattern allows us to extrapolate the form of the wavefunction, hence of the reduced 
density matrix of a chosen neutron pair, after an arbitrary number of scattering 
events. I will spare the reader a detailed account of this derivation, as it is completely 
straightforward but very lengthy. Suffice it to say that, omitting global phase factors, 
one finds the ffist four scattered states of the system as a whole to be 

|^i) = P|01) + Q|10); (6.1) 

1^2) = p2|ooi) + ge"|010) + i?|100); (6.2) 

1^3) = p3|oooi) + ge^^lOOlO) + Pe^lOlOO) + QP^IlOOO); (6.3) 
\ip4) = P^lOOOOl) + ge^'lOOOlO) + Pe^^lOOlOO) 

+ gPV|01000) + gPV|10000); (6.4) 

with 

s = A (^1 + + P, + (6.5) 

P = c^ + d?e-^'^\ (6.6) 
Q = cd{l- , (6.7) 

R = c^d (1 - e-^''^^) . (6.8) 

The basis I have used is an extension of the Bloch basis, where the right-most qubit 
represents the state of the sample and all qubits to the left of this label the neutrons, 
disposed in order of arrival from right to left. It then becomes clear that scattering 
a fifth neutron from the sample will produce the state 

= p5|000001) + ge^'lOOOOlO) + Pe^'lOOOlOO) 

+ gp^e^^ looiooo) + gpV|oioooo) + gp^|iooooo). (6.9) 

Therefore, the general state \iPn) can be written as 

N 

1^^) = p^|l) + e(^-2)^P|4) + Q^'^'e(^"'^^|2^), (6.10) 

i=l, 17^2 
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where the number inside the ket represents the binary value of the state in a single- 
spin-fiip basis. 

In theory, knowing the form of the {ipi) provides a means to calculate the reduced 
density matrix pm.n of any neutron pair, hence the concurrence Cm.n of that pair. 
At first glance, this appears to be quite a task, however things are much simpler 
than they seem. One finds from numerical simulations that the concurrence of any 
neutron and any other can be related to the concurrence of the first neutron and 
those scattered after. Therefore, one can obtain from the set of density matrices pi^„ 
all the information necessary to provide a full characterization of the concurrence 
dynamics of any neutron pair. 

It can be shown that, for all n, the structure of pi^n in the canonical basis is of the 
form 



Pl,n = 



/ 


p 











\ 







Q 


r 












r* 


s 







\ 














/ 



where the coefficients p, q, r and s are obtained from the by tracing over 
all neutrons save the first and the n^^. From the eigenvalues of the spin-fiipped 
counterpart of the concurrence Ci^„ is 

= 2v/^, (6.11) 

with qs = |P"^i|^ \Q\^- Therefore 

Ci,n{N,X,B,,r)=2\P^-'\\Q\\ (6.12) 

as P and Q are in general functions of A^, A, and r. Finally, from numerical 
studies, one finds for the concurrence of an arbitrary neutron pair {m, n} 

(6.13) 

This equation reveals two interesting properties of our system. First, at fixed m, 
the concurrence Cm,n will always fall as \m — n\ increases. Therefore, neutrons 
separated by long time gaps can never be more entangled than neutrons detected in 
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rapid succession. Second, the concurrence of a given pair of neutrons not scattered 
in succession is optimal if one of those neutrons is the first. 

Let us now consider how Cm,n is affected by the neutron interaction time and the 
strength of the apphed field. 



6.2.1 Zero Field Evolution 



In zero field, equation (16.121) becomes 
Ci,„ {N, A, 0, r) = 8N{N + 1)"""^ sin^ 

iV^ + 1 + 2 cos 



2ll + -lAr 



From fl6.13p . the concurrence of an arbitrary neutron pair is then 
Cm,n (N, A, 0, r) = 8N{N + l)"™"" sin^ 

iV^ + 1 + 2 cos 



n-l 
2 



. (6.14) 



2ll + -IAr 



m+n — 2 
2 



(6.15) 



Generalizing the remarks of section 15.1.11 we now see that the concurrence of any 
two neutrons is roughly proportional to A^"^. For all m and n, Cm,n is an oscillating 
function of time, with periodicity T^{N, X,0) [see equation (15.281) ]. The mode of 
these oscillations is determined by the relative values of and m + n. Let us define 
the quantity ( = m + n. If 

{N + lf 



(6.16) 



2N ' 

Cm,n is sinusoidal, and peaks at time tq = Otherwise, the concurrence develops 

a double-peaked structure, whose maxima shift further apart as ( increases [figures 
6.1(a) and 6.1(b)| . The position of the first peak is determined by the relation 



N 



2A {N + r 



sec 



NC 



l + N^-N{C-2)_ 



(6.17) 
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where „ at fixed m is shown in figure 6.1(c) , The peak concurrence of an arbitrary 
neutron pair is, therefore 



if equation fl6.16p is satisfied, or 

QP- 



(C-2 



(6.18) 



W (A^ 



(6.19) 



{N + 1)^ ' 

if it is not. The transition between the two regimes can be observed by comparing 



the form of the purple and orange curves of figure 6.1(a), Interestingly, is 



independent of A^, suggesting that as ( becomes large, the maximum achievable 
entanglement between any two neutrons depends only on the number of scattering 
events which have taken place. The peak concurrence at fixed m for = 10 is 



shown in figure 6.1(d) 



In general, it is not possible to maximize the concurrence of every neutron pair. 
Equations (16.181) and fl6.19p describe the results one might obtain if it were possible 
to tune the interaction time r of individual neutrons. Of course, this cannot be done; 
as discussed in the previous chapter, r is determined by the neutron momentum 
and the size of the sample, and remains constant throughout the realization of the 



protocol. From figures 6.1(a) and 6.1(b) it is clear that an optimal choice of r 



depends on m and n, that is, on the neutron pair we feel most certain of detecting. 
Setting T = To, for example, will maximize the concurrence of {m,n} pairs {1,2}, 
{1, 3}, {1,4}, {1, 5}, {2, 3}, and {2, 4}, but is far from ideal for {1, 15}, {6, 10}, etc. 
This makes us heavily reliant on detecting either the first or the second neutron, 
together with a neutron scattered very soon after. On the other hand, going back 



to figure 6.1(a), by moving towards the outer edge of the curves the concurrence 
may not be so high, but Cm,n falls much less sharply with m and n (see figure . 
Given the form of a similar result can also be obtained by increasing N. 

We conclude from figure \U7I\ that, provided m and n remain small, the chances of 
detecting an entangled neutron pair are reasonably fair. This entanglement may be 
low, but by repeating the protocol a sufficient number of times one could produce 
enough copies of the scattered state to perform distillation. However, this would 
rely on detecting the same neutron pair with every run of the experiment and, as 
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we know, such conditions are difficult to meet. 
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Figure 6.1: (a) The concurrence Ci.„ for n = {2, 3, 4, 5, 6, 7, 15} and N = 10. Increasing values of 
n are indicated by the blue, red, green, purple, orange, yellow and magenta curves, respectively. 
Note C = 6.05, hence the relative structure of the purple and orange curves, (b) As figure (a) 
for N = 20. Note in this case C, = 11.03, hence the relative structure of the yellow and magenta 
curves, (c) The optimal interaction time „ at = 10 and large for m = 1 (blue dots), to = 2 
(red dots), m = 3 (green dots), to = 4 (purple dots) and to = 5 (orange dots), (d) The peak 
concurrence as a function of n at A = 10 for m = 1 (blue dots), to = 2 (red dots), to — 3 (green 
dots) and to = 4 (purple dots). All quantities arc expressed in natural units, with A = 1. 
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Figure 6.2: (a) The concurrence Ci^„ for iV = 10 at r = 7^/2, t = 0.8 and r = 0.6 (blue, red and 
green curves, respectively), (b) Ci_„ for iV = 20 at r = 10/2, r = 1 and t = 0.8 (blue, red and 
green curves, respectively). All quantities are expressed in natural units, with A = 1. 



6.2.2 Finite Field Evolution 

In finite field also, the evolution of the concurrence of any neutron pair can be traced 
back to Ci^n- Let us then consider the behaviour of Ci,„ as a function of the applied 
field and the neutron interaction time, as shown in figure I6.3[ A slice through these 
curves at any point along the axis illustrates the concurrence remains a periodic 
function of r irrespective of the field (figure WM . does, however, determine the 
timescale of the oscillation, which is given as always by equation (15.281) . Comparing 
the shape of the curves for different values of n, one notes two main features. First, 
the peak concurrence falls with n. Second, the essential attributes of the curves do 
not change a great deal as n is increased. In all cases, the concurrence at short 
times is characterized by a concave, approximately conical dip, which expands into 
a 'valley' as n becomes large. By inspection, maximum values of Ci^„ lie along the 
'mouth' of this valley, which is mirror-symmetric about Bz = B^. We know this field 
maximizes the concurrence of the first two neutrons; therefore, from the form of Ci^„ 
and the trends shown in figure we conclude that 5* represents the optimal field 
value for all m and n. Then, with the help of equation (I6.13P and setting B^ = B^, 
Cm,n takes on the particularly simple form 



Cm.n — 2 sin 



2Ar 



cos'"+"-' I —= I 



2Ar\ 



(6.20) 
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which peaks at time 



T* = 




m + n — A 



) 



(6.21) 



4A 



cos 



m + n 



As expected, both of these expressions reduce to fl5.36p and fl5.37p when m = 1 and 
n = 2. Finally, from (16.201) and (I6.2ip . the maximum concurrence of any neutron 
pair is given by 



which, for ease of comparison, I have written in terms oi C, = m + n. Indeed, a 
quick glance to section 16.2.11 shows equations (I6.22p and (I6.18P are identical. This 
highlights an unexpected property of the system: applying a magnetic field can only 
improve the maximum concurrence of a given neutron pair with respect to its zero 
field value if the number of scattered neutrons is, roughly speaking, smaller than 
the number of spins in the sample. Above this, we run up against a fundamental 
limit, whereby the maximum concurrence is independent of the external parameters. 
Therefore, represents an upper bound to the amount of entanglement any two 
neutrons can acquire. It is important to note that equations (16.211) and (I6.22p do 
not in general describe the performance of the protocol in finite field; rather, they 
define how the concurrence of a given neutron pair would scale if neutron interaction 
times were individually tuneable. 
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(a) (b) 




Figure 6.4: A slice through Ci.„ at randomly chosen values of (different for each plot) for initial 
state \ipo)- AH quantities are expressed in natural units, with A = 1. (a) rt = 2; (b) n = 3; (c) 
n = 4; (d) n = 5; (e) n = 6; (f) n~ 7. 



Let us then consider the concurrence of the pair {m, n} at an arbitrary value of B^. 
This has the form 

CmAN,X,B,,T) = 32iV^ sm' 0r / ^^^^^ ^ ^^^^ ^ ^^^^^^ ^^^^^.^^ 

(6.23) 

where and if are given by equations (I5.14p and (15.151) . respectively. This is the 
generalization of equation (I5.27P for all m and n. We know from section 15.1.21 that 
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within the interval G [B^- , 5z+] the concurrence Ci^2 develops a double-peaked 
structure; if m, n 7^ [1? 2], the limits of this interval are given by 



B-. 



B. 



A 1 





/2(C-2) 


N 


V N 




/2(C-2) 




V N 



(6.24) 
(6.25) 



Within this interval, Cm,n peaks at time 



1 



COS 



16A^ + iVV - 4A^AV' (C - 2) 



4iVAV^C 

Otherwise, maxima occur at T0/2. Therefore, one finds 



(6.26) 



a = 4i 



(C-2) 



C-2 



(6.27) 
(6.28) 



For a chosen neutron pair, these relations define the maximum achievable concur- 
rence for any value of N and B^. As in the previous section, these limits can be 
phrased in terms of the number of scattered neutrons, rather than the applied field. 
Defining 



c 



B 



1 1 

2iV ^ 2A 



2B, 



N{B,-\y 

A 



one has 



if C < Cb, and 



otherwise. 



CP — C 



nv — n. 



(6.29) 

(6.30) 
(6.31) 



Finally, to assess the performance of the protocol in a finite field, let us consider 
how the concurrence scales with m and n. If the field is weak, (b is proportional to 
A^. Therefore, the optimal interaction time at fixed N will depend on the number 
of scattered neutrons relative to the number of spins in the sample. If this ratio 



Chapter 6. Many-Neutron Scattering 



122 



is small, the optimal time is the same (i.e. for all m and n. Otherwise, as 

discussed in section 16.2.11 one must choose a time which yields appreciable values of 
Cm,n for the broadest possible interval \m — n\. As the field is raised above -Bz-, this 
optimization strategy becomes compulsory to avoid the concurrence plummeting 
after the first few scattering events. When Bz exceeds -Bz+, Cb becomes large and 
optimal measurement times come once more to coincide. However, the concurrence 
of any neutron pair is extremely low no matter how values of r are chosen. These 
results are shown in figure 16.51 

In conclusion, if the sample is prepared in a single-magnon state the performance 
of the protocol can be optimized by applying a magnetic field within the range 
[Bz-i Bz+]. In general the concurrence of an arbitrary neutron pair decreases with 
the total number of scattering events. However, by choosing the neutron velocity 
with care, one can ensure this decline is moderately slow. Therefore, the likelihood of 
detecting an entangled neutron pair is reasonably strong, though it remains desirable 
to detect neutrons scattered in quick succession. Peak values of the concurrence do 
not typically exceed 0.5 unless one of the neutrons detected is the first. Consequently, 
the yield of highly entangled pairs could be improved by periodically re-setting the 
sample, i.e re-initializing the system to the original input state. 
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Figure 6.5: The concurrence Ci^„ in different field regimes for A^=20 and n = {1, 2, 3, 4, 5, 6, 7, 15}. 
Increasing values of n are indicated by the blue, red, green, purple, orange, yellow and magenta 
curves, respectively. All quantities are expressed in natural units, with A = 1. (a) Ci_„ at = 0.3. 

(b) Ci,n at Bz = 0.3 and r = (blue curve), r = 0.9 (red curve), and t = 0.7 (green curve). 

(c) Ciji at Bz = 0.9. (d) Ci^„ at B^ = 0.9 and r = Tg'^ (blue curve), t = 1.7 (red curve), and 
T = 0.8 (green curve), (e) Ci_„ at B^ = 3. (f) Ci,„ at = 3 and r = (blue curve), r = 0.6 
(red curve), and r = 0.4 (green curve). 
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6.3 Scattering from a fully Polarized Sample 

Scattering from a fully polarized sample, with the system initialized to state I^Pq) 
[see equation fl3.9l) ]. cannot be described analytically, as the system wavefunction 
becomes increasingly convoluted with each step of the protocol. This occurs because 
the system can now accommodate as many spin flips as scattered neutrons, therefore 
the size of the working Hilbert space grows linearly in the total number of incoming 
neutrons. Consequently, the concurrence must be characterized using numerical 
methods. Owing to the increased computational demands as the number of scattered 
neutrons becomes large, I will limit my study to short interaction times bounded 
by r = 20 (expressed in natural units, and for A = 1). As this interval suffices to 
observe the salient features of the concurrence of the first two neutrons, and the 
interaction of each neutron with the sample is the same, it is not unreasonable to 
assume that any interesting features of the function Cm,n might emerge within this 
time frame. 

At fixed A^, we expect the concurrence of a given neutron pair to depend on three 
parameters 

1. The interaction time; 

2. The applied field; 

3. The neutron polarization. 

Let us then observe how changing one or more of these parameters can affect the 
performance of the protocol. 

6.3.1 Zero-Field Evolution 

It is already known from section [5.2.1l that when the applied field is zero the concur- 
rence of the first two scattered neutrons is a decaying function of A^. To study the 
concurrence of an arbitrary neutron pair, it is therefore convenient to work at fixed 
A^, chosen such that the peak value of Ci^2 is non-zero within the allowed range of 
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neutron polarizations (see figure [5^ . and is large enough to accommo- 

date several spin fiips. I work with = 10, for which the peak concurrence of the 
first neutron pair falls just below 0.1 (see figure . 

It is useful to begin with a brief recap of the properties of Ci,2- These are summa- 
rized in two points: (i) oscillatory behaviour as a function of interaction time; (ii) 
maximum amplitude for a = Let us then observe the concurrence of different 
neutron pairs at a = ^ ; this is shown in figure 16.61 It is immediately evident that 
the concurrence maintains its characteristic periodicity irrespective of the indices 
m and n, however these indices do affect the mode of oscillation. At fixed n, the 
concurrence falls with m and, with the exception of m = 1, the peak value of Cmn 
falls with increasing \m — n\. More interesting features develop as a It now 

emerges that the concurrence of distant {\m — n\ > 2) neutrons can exceed that 
of neutrons scattered in succession (figures 16.7116. Ill) . The most prominent example 
occurs at m = 1 and = for which the concurrence of the first and second neu- 
trons represents a minimum value. This may seem promising, however the neutron 
concurrence remains low no matter which neutron pair we detect. Furthermore, 
even for an optimal choice of input parameters, the resulting state of the scattered 
neutrons is not compatible with the witness decomposition of equation (15.991) . 
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Figure 6.6: The concurrence of different neutron pairs for incoming neuron polarization a = 
/? = ^ and -B^ = 0. Increasing values of n are indicated by the blue, red, green, purple, orange 
and yellow curves, respectively, (a) The concurrence of neutrons m and m + 1 for to < 4. (b) The 
concurrence Ci„ for 1 < n < 5. (c) The concurrence C2n for n = {1,3,4,5}. (d) The concurrence 
Csn for n = {1,2,4,5,8}. (e) The concurrence C4„ for n = {1,2,3,5,7,10}. All quantities are 
expressed in natural units, with A = 1. 
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Figure 6.7: The concurrence Cmn at fixed r and -B^ = 0, for m G [1,4] and increasing n. Values 
of T are inset, with the exception of figure (a), where r = Tfj,/2. All quantities are expressed in 
natural units, with A = 1. Results are for incoming neuron polarization Q = -^,/3=i. (a)m = l. 

(b) TO = 2. (C) 771 = 3. (d) 771 = 4. 
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Figure 6.8: The concurrence of different neutron pairs for incoming neuron polarization a ~ 
/? ^ and = 0. Incr easing values of n are indicated by the blue, red, green, purple, orange 
and yellow curves, respectively, (a) The concurrence of neutrons m and m + 1 for m < 5. (b) 
Cin for 1 < n < 5. (c) C2„ for n = {1,3,4,5}. (d) Csn for n = {1,2,4,5,7}. (e) dn for 
n = {1, 2, 3, 5, 8, 10}. All quantities are expressed in natural units, with A = 1. 



Chapter 6. Many-Neutron Scattering 



129 



0.04 . 
0.02 

0.00 



(a) 



(b) 




(c) 




(d) 




Figure 6.9: The concurrence Cmn at fixed t and = 0, for m € [1,4] and increasing n. Results 
arc for incoming neuron polarization a = /3 = (a) m = 1 at t = Tif,/2. (b) m = 2 at r = 1 
and T = (blue and red dots, respectively), (c) m = 3 at t = 0.65, t = ?0/2, and t ~ 2 (blue, 
red and green dots, respectively), (d) m = 4 at r = 0.7. All quantities are expressed in natural 
units, with A = 1. 
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Figure 6.10: The concurrence of different neutron pairs for incoming neuron polarization a = 
P = ^ and Bz = 0. Increasing values of n are indicated by the blue, red, green, purple, orange 
and yellow curves, respectively, (a) The concurrence of neutrons m and m + 1 for m < 4. (b) 
Cin for I < n < 5. (c) C2„ for n = {1,3,4,5}. (d) Csn for n = {1,2,4,5,7}. (e) dn for 
n = {1, 2, 3, 5, 8, 10}. All quantities are expressed in natural units, with A = 1. 
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Figure 6.11: The eoncurrence Cmn at fixed r and i?^ = 0, for m e [1, 4] and increasing n. Results 
are for incoming neuron polarization a — /3 = (a) m = 1 at r = (b) m = 2 at r = '7^/2. 

(c) TO = 3 at T = 7^/2, T = 0.75 and t = 1.825 (blue, red and green dots, respectively), (d) 
TO = 4 at r = 0.5, r — 0.8 and t = 2.15 (blue, red and green dots, respectively). All quantities are 
expressed in natm'al units, with A = 1. 



6.3.2 Finite-Field Evolution 

When a finite field is applied, the behaviour of the concurrence is highly unpre- 
dictable. Comparing the concurrence of an arbitrary pair {m, n} with Ci 2 at 
Bz < Bz+1 it is no longer a given that the peak value of Cm,n will increase if 
Bz is raised [see section 15.2.21 and figure 6.6(b)| . As expected, the details of the time 



evolution of Cjn,n are determined by the strength of the field and the neutron po- 
larization; however, exactly how these parameters affect the concurrence is a priori 
unclear. To make our task more manageable, let us use the strength of the external 
field to define three regimes: 



1. Weak field: < |; 



Chapter 6. Many-Neutron Scattering 



132 



2. Strong field: B, > 3A; 

3. Optimal field: = 5*, 

and study the features of the concurrence in each case for different neutron polar- 
izations. 

In the weak field regime, the concurrence of the first and second neutrons shows a 
well-defined periodicity as a function of the interaction time (cfr. figure ISTTSl) . As m 
and n are increased, this is no longer true unless a = 0. In general, the concurrence 
falls with increasing \m—n\, though when m > 2 instances of Cm,n > Cm,m±i become 
more common. 

A strong dependence on the neutron polarization is also observed and, in analogy 
with the results of section 15.2.21 the performance of the system markedly improves 
as a (figs. I6.12H6.151) . Indeed, the a = case merits a separate discussion, as 
several interesting effects are observed. First, as noted above, the concurrence of 
any neutron pair shows traces of the periodic behaviour observed for Ci^2- Second, 
for any m and n the concurrence evolves through a series of peaks. At fixed m, 
peaks in Cm,n for different n are aligned [see figure 6.15(a)| , though as m increases 



the position of the strongest peak can change. In correspondence of this peak, the 
concurrence at fixed m falls as (n + 1)"^. The scattered state of the neutrons at 
this time is always compatible with the witness decomposition of equation (15.991) . 
however the concurrence of any non-successive neutron pair never exceeds 0.6. This 



is shown in figure 6.15(d) 
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Figure 6.12: The concurrence Cmn for = A/2 and a = at (a) m = 1 and (b) m = 2. 
Increasing values of n are indicated by the blue, red, green and purple curves, respectively. All 
quantities are expressed in natural units, with A = 1. 




Figure 6.13: The concurrence Cmn for Bz ~ A/5 and a = at (a) m — 1, (b) to = 2, (c) to = 3 
and (d) to = 4. Increasing values of n are indicated by the blue, red, green, purple and orange 
curves, respectively. All quantities are expressed in natural units, with A = 1. 




Figure 6.14: The concurrence Cmn for = A/3 and a = at (a) m = 1, (b) to = 2 and (c) 
TO = 3. Increasing values of n are indicated by the blue, red, green, purple and orange curves, 
respectively. All quantities are expressed in natural units, with A = 1. 
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Figure 6.15: The concurrence for = A/3 and a = 0. Increasing values of n are indicated by 
the blue, red, green and purple curves, respectively. All quantities are expressed in natural units, 
with A = 1. (a) Ci^n- (b) C2,„. (c) C3,„. (d) The peak concurrence for m = {1,2,3} (blue, red 
and green dots, respectively) at r = 17 (blue dots) and t = 13.6 (red and green dots). Dashed 
lines indicate the function F(n) = , °, , + 6, evaluated at a = 2.49 and b = —0.75 (blue curve), 
a — 1.42 and b = —0.37 (red curve) and a = 1.11 and b = —0.29 (green curve). 



If the field is strong, the concurrence is generally low - maximum values do not 
exceed 0.2. Excluding for the moment the case a = 0, for an arbitrary neutron 
polarization Cm,n retains some traces of periodicity provided ??7 < 2. Otherwise, the 
concurrence is a disordered function of the interaction time. As in weak field, peak 
values of Cm,n increase as a ^ 0, but fall with \m — n\. However, is now extremely 
difficult to find instances of r and n for which the concurrence at fixed m increases 
with respect to Cm,m±i- As opposed to Ci^2, if m > 1 the concurrence of a given 
neutron pair can improve with respect to its zero-field value, even if Bz exceeds the 
threshold Bz+. However, the concept of a limiting field still applies, as it is not 
possible to increase the peak concurrence indefinitely by making the field stronger 
(see figure 16.161) . 
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Figure 6.16: The concurrence of different neutron pairs for incoming neuron polarization a = 
/? = ^. Increasing values of n are indicated by the blue, red, green, purple and orange curves, 
respectively, (a) The concurrence of neutrons m and m + 1 for m < 5 and ~ 3A. (b) As (a) 
for ^ 4A. (c) Can for n = {1,3,4,5} and B^ = 3A. (d) As (c) for B^ ^ 4A. (e) Cg^ for 
n — {1,2,4,5,7} and B^ — 3A. (f) As (e) for B^ — 4A. All quantities are expressed in natural 
units, with A = 1. 



As always, unusual features are observed at a = 0. In this case, Cm,n retains the 
same periodicity as Ci^2 for all m and n (see figure I6.17p . In analogy to its weak- 
field behaviour, the concurrence evolves through a series of peaks which, at fixed m, 
are aligned for different n. If r is chosen to coincide with a peak, the concurrence 
at fixed m falls with a very gentle gradient as n~^. This gradient does not change 
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significantly as m is increased, therefore the concurrence of neutron m with any other 
is generally a slowly- varying function of n. In other words, provided |m — ?i| remains 
sufficiently small, the concurrence of any neutron pair is roughly constant (see figure 



6.17(d)). Consequently, for a group of neutrons scattered in quick succession, the 



entanglement with respect to any bipartition is also approximately constant, and 
can be detected using the witness of equation 05.991) . 
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Figure 6.17: The concurrence for Bz = 3A and a ~ 0. Increasing values of n are indicated by 
the blue, red, green and purple curves, respectively. All quantities are expressed in natural units, 
with A = 1. (a) Ci,„. (b) C2,„. (c) C3,„. (d) The peak concurrence for m = {1,2,3,7} (blue, red, 
green and purple dots, respectively) at r = 6.4. Dashed lines indicate the function F{n) — an + b, 
evaluated at a = —0.006 and b = 0.164 (blue curve), a = —0.005 and 6 = 0.157 (red curve), 
a = —0.005 and b = 0.152 (green curve), and a = —0.004 and b = 0.127 (purple curve). 



Finally, let us look to the properties of Cm,n at 5^ = = A (l — -^). On the 
basis of our discussion so far, it is reasonable to assume the protocol will be most 
successful at a = 0, therefore I will focus on this limit. It was shown in section 
15.1.21 that the concurrence of the first two scattered neutrons is a maximum when 
the field is set to i?*. At general m and n, this is no longer the case, because the 
concurrence of neutron n with any neutron scattered before it will depend on the 
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sample state it encounters. For the case of Ci^2, represents an optimum because it 
allows us to 'deposit' entanglement in the system by creating a maximally entangled 
state of the first neutron and the sample. For all practical purposes, at this stage of 
the protocol the sample behaves like a two-level system; it is therefore conceivable 
that the second scattering event might swap - at least in part - the state of the 
second neutron and the state of the system. After all, the entanglement has few 
other places to go. As the sample acquires more spin-flips, several effects come into 
play. First, the entanglement we deposit acquires new channels over which to re- 
distribute itself. Second, the sample ceases to behave like a qubit, and becomes a 
bona fide multilevel system. At that point, even if neutron n were to encounter a 
fully entangled state of neutron m and the sample, the protocol could only succeed if 
most of the sample wavefunction were concentrated in the single-excitation subspace. 
As the number of events grows, the chances of this happening become very slim. 
Therefore, in the context of this more general version of the protocol, -B* simply 
defines an 'intermediate' field regime. As one might expect, the concurrence in this 
regime shows a mixture of strong- and weak-field features. Among the former, we 
note the existence of interaction times for which the concurrence of any neutron 
pair is linearly dependent on \m — n\ [see figure 6.18(b)| . Among the latter, we find 
instances in which the concurrence of a non-consecutive pair is greater than that of 
two neutrons scattered in succession [see figure |6.18(c)| . 

In summary, if the sample is prepared in a fully polarized state, the concurrence of 
neutrons m and n falls with the number of scattering events between them. In some 
cases, it is possible to predict how Cm,n scales with |m — n|, however the evolution 
of the concurrence is generally unpredictable, unless either the applied field is zero, 
or the neutrons are polarized in the negative ^-direction. In the former case, the 
concurrence of any neutron pair is low; in the latter, the performance of the protocol 
depends on the strength of the applied field. Best results are obtained in the weak- 
field regime, where the concurrence at fixed m scales roughly as n~^. Maximum 
values of the concurrence of any non-consecutive pair do not exceed 0.5, hence in 
this respect a realization with input state is equivalent to a realization with 
input state What we lose is the ability to predict with certainty when the 

concurrence of a given pair will peak, although from the behaviour of Cm,n with 
r (cfr. figure I6.15P one could formulate a rough guess. Unfortunately, given the 
non-analytical form of Cm,n in weak field, the margin of error is rather narrow. 
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Figure 6.18: (a) The concurrence Ci^„ for Bz — 0.9A and a = 0. Increasing values of n are 
indicated by the blue, red, green and purple curves, respectively. All quantities are expressed in 
natural units, with A = 1. (b) The concurrence Cm,n at t = 14.41 for m = 1 (blue dots), m = 2 
(red dots), to = 3 (green dots), and m = 4 (purple dots). Note the linear decrease with n at fixed 
m. (c) Ci,„ at T = 13.1. Note Ci,4 > Ci,5 > Ci,2. 

6.4 Average yield of Entangled Pairs 

I have shown so far that the entanghng protocol of chapter [S] can be effective for a 
range of neutron indices m and n, on the proviso that |m— ?7,| is not large compared to 
the number of spins in the sample, A^. To assess the true performance of the protocol 
in a realistic situation, any reference to the indices of the scattered neutrons must 
now be removed. The question to ask becomes: will a pair of neutrons chosen at 
random from those which have scattered be entangled? 

The results of sections 16.21 and 16.31 indicate that the entanglement of two neutrons 
m and n decreases with \m — n\. Let us therefore assume we apply a resetting pulse 
to the sample whenever the number of scattered neutrons approaches A^, such that 
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|m — n| < A^. I introduce the quantities t and C to represent the number of scattered 
neutrons on which measurements are performed, and the concurrence of the average 
state of any pair, respectively. For ease of comparison, I will compute the average 



concurrence relative to the situations presented in figures 6.2(a), 6.2(b), 6.5(b) 



6.5(d), 6.5(e) [6.151 16.171 and 16.181 These all correspond to optimal realizations of 
the protocol at fixed and B^, either with input state H^q) (first five cases) or with 
input state (latter three cases). As expected, the average concurrence falls 

with the size of the sampling ensemble. However, some slightly surprising features 
are observed, according to the choice of input state and the strength of the applied 
field. 

If the sample is prepared in a single- magnon state (input state \iPq)), the likelihood 
of detecting an entangled pair is strongly dependent on the value of B^. If B^ 
is zero, C" is a slowly-varying function of t - in fact, at certain values of r the 
C vs. t curve shows an almost linear trend (see figure 16.191) . Consequently, the 
average concurrence remains finite for an appreciable range of t, and shows signs 
of converging to a lower limit as t A^. As the field is increased, the relationship 
between C and t changes rather dramatically, approaching a near-exponential decay 
in the intermediate field regime [figure 6.20(b)| . By the time B^ reaches B2+ the 



curve has somewhat flattened out, but values of C are extremely low, and fall to 
zero by t = 5. Best results are therefore achieved at zero or small fields. 

Conversely, if the sample is prepared in a fully polarized state, it makes little sense 
to consider an average concurrence unless a field is applied. This is due to the 
fact that, in zero field, the concurrence is extremely low for all m and n. The 
optimal operating mode for the system now seems to be the weak field regime, in 
which, for an optimal choice of r, the average concurrence falls as t~^. As above, 
the situation deteriorates as the field is raised, and in the intermediate and strong 
field limits the behaviour of the system is similar to that observed for input state 
\iPq). We conclude that, for either input state, the probability of a measurement 
on two randomly chosen chosen neutrons yielding an entangled state is substantial, 
subject to a set of experimental conditions. First, the neutron interaction time must 
be chosen appropriately. Second, the field must remain below the threshold -B^-. 
Third, if input state \ipo) is chosen, the neutron spin must be antiparallel to the 
quantization axis. 
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(a) (b)^ 




Fiffure 6.19: Input state |?/'^). (a) The average concurrence of any neutron pair within the range 
{m,n} G [1,10] for = 10 and = 0, sampled at times r ~ 0.8 and 0.6 (blue, red and 

green curves, respectively), (b) The average concurrence of any neutron pair within the range 
{m,n} G [1,20] for = 20 and ~ 0, sampled at times t = 1 and 0.8. All quantities are 

expressed in natural units, with A = 1. 
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Figure 6.20: Input state iV'c?)- The average concurrence of any neutron pair within the range 
{m, n} e [1, 20] for A^ = 20. Values of and r are as follows: (a) B^ = 0.3A, t = 0.9 and 

0.7 (blue, red and green curves, respectively), (b) B^ ~ 0.9A, r — Tq^, 1.7 and 0.8. (c) B^ = 3A, 
T = T0/2, 0.6 and 0.4. All quantities are expressed in natural units, with A = 1. 
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Figure 6.21: Input state \tpQ). The average concurrence of any neutron pair within the range 
{m,n} G [1, 10] for TV = 10. Values of and r are as follows: (a) Bz = A/3, t = 17 and 13.6 
(blue and red curves, respectively), (b) B^ = 0.9A, r = 14.41. (c) Bz = 3A, r = 6.4. The black 
dashed line in figure (a) corresponds to the function F{t) — at^^ + h, with a — 1.63 and b ~ —0.44. 
All quantities are expressed in natural units, with A = 1. 



6.5 Conclusions 



It was shown in chapter [5] of this thesis that two neutrons consecutively scattering 
from a sample prepared in a specific eigenstate of its internal Hamiltonian could 
come to share a substantial degree of entanglement, subject to an appropriate choice 
of the experimental parameters. Extending this protocol to include the — more 
realistic — possibility of multiple scattering events, it emerges that any two neutrons 
scattered within a certain time frame of each other can be measurably entangled 
in a similar fashion. The absolute value of this entanglement is generally low, and 
falls well short of the threshold Cs = 0.67 defined in section 15. 51 However, the 
scheme remains of great interest from a fundamental viewpoint, providing as it does 
a simple, practical means to answer the as-yet unresolved question: can two distinct 
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neutrons be entangled? 

In general, the outcome of the protocol is determined by the number of scattered 
neutrons relative to the number of scatterers. Provided this ratio remains small, 
there is a strong possibility that measurements made on any two scattered neutrons 
will yield an entangled state. Optimal performance can therefore be achieved by 
periodically resetting the sample to its initial state. In these circumstances, the 
average concurrence of any neutron pair is of order 0.2, corresponding to just under 
0.1 ebits of entanglement. This can be boosted either by specifying which neutron 
pair to measure, or by increasing the frequency of the resetting pulses. Distillation 
may also be a possibility, but only if one could ensure detection of the same neutron 
pair subsequent to each re-initialization of the sample state. 

A more detailed analysis of the distribution of entanglement within the scattered 
neutron bunch further shows the concurrence is finite with respect to several bi- 
partitions. Provided the number of scattered neutrons is small, one therefore creates 
genuine multi-partite entanglement. However, this particular aspect of the protocol 
has not been explored in great detail. 

The experimental requirements of the scheme were assessed in section 15.51 where 
it was concluded that measurement of the witness of equation fl5.99p may be ex- 
tremely challenging, but not a priori impossible. This is reinforced by the results of 
the present chapter. Therefore, the only outstanding problems remain the neutron 
velocity and the neutron coherence length. The former issue was discussed in detail 
in the previous chapter. As regards the neutron coherence volume, we recall this 
quantity depends on the degree to which a beam can be monochromated. A perfect 
crystal monochromator can achieve velocity resolutions of Av/v = 10~^ [125] . but 
for our purposes this is not enough. Higher resolutions may be achievable by select- 
ing neutrons from the beam with the aid of an extremely high frequency chopper. 
This has not been reported to date, but should not be impossible, as there exists no 
fundamental lower limit to Av/v. Hence, it is reasonable to conclude that meeting 
the final requirement of the protocol may simply be a matter of time. 



Chapter 7 



The Problem of Scalability and 
the Need for Communication 

In this chapter, I introduce the topic of quantum state transfer with spin chains. The 
use of a spin chain as a communication channel was first suggested by Bose I126f . in 
a bid to find an efficient means of transferring quantum states over short distances in 
solid-state devices. It was found that unmodulated spin chains are natural conduits 
of quantum information, but do not always allow high quality communication. Many 
improvements to the basic communication protocol have since been suggested, and 
these will be summarized. I will then explain the original communication scheme 
in detail, and review its performance. Possible experimental implementations will 
also be discussed, ending with a general overview of spin chain communication as a 
quantum computational technique. 



7.1 Introduction 



At the very beginning of this thesis, it was observed that communication represents 
the sharing of a correlation. In practice, this correlation arises when information 
flows from one party to another via some kind of mediating channel. There exist 
a wide variety of channels capable of transmitting classical information, but those 
same channels are not in general adequate for quantum state distribution. Indeed, it 
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has been shown that transferring quantum information through a classical channel 
involves a loss of at least a third of the original information |127] . This is un- 
acceptable for protocols such as teleportation or quantum key distribution, which 
require faithful exchange of qubits. Currently, the most common means of trans- 
ferring quantum information is the optical fibre, a choice born of necessity as most 
large-scale quantum computational tasks attempted to date have employed photons 
as qubits. This is down to the photon's low interaction cross-section, which enable 
it to remain in a coherent state over large distances both in air and in optical fibres. 
However, in recent years the problem of quantum state transfer in solid-state devices 
has attracted a great deal of attention; in fact, it has been suggested this may be 
the key to building quantum computers large enough to rival existing processors. 

Just as the computing power of a classical machine is directly related to the number 
of transistors on its Central Processing Unit (CPU), the computing power of a 
quantum computer is directly proportional to the number of qubits it can access. It 
has been calculated that a 10^-qubit quantum computer running with a clock speed 
of 100 MHz could, in just one hour, solve a factorization problem that would keep 
the most powerful existing supercomputer busy for a time comparable to the age of 
the universe. Unfortunately, this is at least a million times larger than the machines 
currently within our ability to build [128[ 1129] . 

Scalability is therefore a key requirement for a feasible quantum computer architec- 
ture. Up-scaling current designs is not a straightforward task, for two main reasons. 
First, many-qubit systems acquire extra degrees of freedom, which are difficult to 
control and provide a path for spurious interactions with the environment. As a 
result, the system becomes more vulnerable to decoherence and, consequently, a 
less reliable quantum computer. A second problem relates to the implementation of 
two-qubit gates between distant qubits. To achieve this, the qubits must interact 
either directly or by a common bus mode, but neither approach is without its prob- 
lems. If we aim for direct couplings, the qubits in question (or, alternatively, the 
information they encode) must be shuttled around the system until they are close 
enough to interact, and then returned to their original positions. This is clearly 
inefficient, and would slow the calculation down considerably. On the other hand, 
if one employs a common bus mode there is a physical limitation to the number 
of qubits that can be coupled. These issues become increasingly problematic as 
the size of the computer increases, thus rendering the construction of a many-qubit 
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'processor unit' an extremely challenging task. 

7.2 A Review of Quantum Communication with 
Spin Chains 

How else could quantum processors be scaled up, if not by simply adding extra 
qubits? The past few years have witnessed the advent of a 'divide and conquer' 
philosophy. This entails building a modular processor composed of several intercon- 
nected few-qubit registers, which could either act as parallel processors, or perhaps 
have different functions (one register could be the CPU, another the memory etc.). 

Such a 'collaborative' approach clearly requires fast, efficient communication be- 
tween the registers. The best known and most widely used carriers of quantum 
information are photons, because they travel naturally and are easily isolated from 
their environment. Photons can be successfully propagated over distances of hun- 
dreds of kilometers through optical fibres, free space - even through outer space, 
as recently demonstrated by a team of researchers led by Zeilinger |13m 11311 1132] . 
One might then suggest constructing a simple quantum computer by connecting a 
series of solid-state processors by a mesh of optical fibres, but until very recently this 
approach seemed impractical, owing to the difficulty of interfacing photonic qubits 
with solid-state devices. It has now been shown that strong coupling can be achieved 
in a variety of systems |133[ 1134] , which may open up the possibility of exploring 
this avenue. The option of a quantum computer with 'all optical' components has 
been looked into, but such designs are complex and often involve non- linear pro- 
cesses, which are very challenging to implement experimentally (see |135l I136[ 1137] 
and references therein). 

Recently, much effort has been dedicated to studying systems in which quantum 
information is encoded in stationary qubits and propagated from one part of the 
system to another by the natural interaction between the system's components. One 
of the simplest geometries in which this can be achieved is a one-dimensional chain 
of interacting particles, where the qubit is encoded in the particles' spin degree of 
freedom. This is usually referred to as a spin chain. The original proposal for 
quantum communication through spin chains was put forward by Bose |126] . who 
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studied a ferromagnetic, exchange-coupled chain of spins with a constant nearest- 
neighbour interaction, in an external magnetic field. It was found that information 
could be transferred with a fidelity exceeding the maximum classical value in a time 
that grows polynomially with the length of the chain. The importance of this result 
lies in the fact that the transfer protocol requires no external control whatsoever. 
All the experimentalist must do is encode a quantum state into the first spin of the 
chain, wait for a well-defined amount of time and measure the last spin of the chain. 
Depending on the chain length, the fortunate experimentalist will have recovered 
a more or less perfect copy of the initial state without having to intervene in the 
transfer process in any way. 

It was subsequently shown that slightly more complex systems allowed better per- 
formance. Osborne and Linden |138j demonstrated the viability of transmission 
fidelities close to unity if the qubit is taken to be a carefully designed 'wave packet', 
provided the sending and receiving parties can access a sufficiently large portion 
of the chain. In the absence of structural imperfections, an XY Hamiltonian on a 
hyper-cubic lattice allows perfect state transfer [1391 11401 1141j . as do parallel spin 
chains |142j . Arbitrarily perfect state transfer can also be achieved by applying a 
sequence of two-qubit gates at the receiving end of a quantum chain [143] , or by em- 
ploying auxiliary local memories [144j . Plenio et al. [145j have studied the situation 
for chains of harmonic oscillators (i.e. where each particle on the lattice possesses a 
continuous, rather than a discrete, degree of freedom), while Hartmann et al. |146j 
have recently found that quantum information can be made to propagate with ar- 
bitrarily high fidelity through both oscillator and spin chains near a quantum phase 
transition. This transfer is exponentially slow; more rapid transmission is possible 
at the quantum critical point, but at some cost to the fidelity. High fidelities can 
also be attained by engineering the strength and the nature of the interactions be- 
tween the spins |1471I148] . or by performing well-defined projective measurements on 
parts of the chain [1491 115U] . Perfect state transfer is even achievable in disordered 
or imperfect systems, and without state initialization, although extra resources or 
error correction may be required [1491 11511 1152] . 

The vast majority of these protocols assume only short-range couplings between 
the spins in the system. More recently, Kay [153] has made a detailed study of 
finite spin chains in which the total Hamiltonian accounts for a long-range magnetic 
dipole potential. It was found that state transfer can be optimized by fine-tuning 
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the inter-particle spacings and applying local magnetic fields at each spin site. 

Some effort has also been dedicated to examining the effect of thermal fluctuations 
on the transfer process |154] . For short chains, there exist analytical expressions 
which demonstrate that thermal fluctuations affect the fidelity of the transfer, but 
not its speed. For longer chains, the problem rapidly becomes intractable because 
the dimensions of the Hilbert space of the system blow up exponentially as the 
number of qubits increases. In most state transfer protocols suggested to date, this 
difficulty is avoided by choosing Hamiltonians with built-in symmetries, which allow 
one to work in a Hilbert space that grows linearly rather than exponentially in the 
number of spins (see section 17.3. ip . These symmetries cannot be exploited if dissipa- 
tive processes such as thermally induced excitations are present, therefore thermal 
effects are often neglected. This is not necessarily an unrealistic approximation, as 
it equates to considering thermal energy scales kT which are negligible compared 
to the strength of, say, an external magnetic field or the inter-qubit couplings. One 
could certainly engineer such a situation in a laboratory environment. 

7.3 The Communication Protocol 

We now review the communication protocol originally suggested by Bose |126] . The 
present discussion deals with a linear chain geometry, however the scheme is valid 
for an arbitrary graph of qubits. 

A sending party, who in the name of tradition we call Alice, wishes to transmit 
an unknown quantum state to a receiving party, surprisingly named Bob, who is 
located at the opposite end of the chain. The transfer can be achieved with minimal 
effort following four very simple steps: 

1. The chain is initialized to its ground state. 

2. Alice prepares a qubit in the state she wishes to transmit and replaces the first 
spin in the chain with this qubit. 

3. The chain is left to evolve for a finite amount of time under the effect of its 
internal Hamiltonian. 
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4. Bob measures the last spin in the chain and recovers a copy of the state 
originally prepared by Alice. 

The key stage of this protocol is point 3: the information encoded by Alice is 
transferred to Bob spontaneously by the inter-spin couplings. As we will see in the 
following sections, information has flowed from Alice's end of the chain to Bob's. 
The chain therefore acts as a channel for quantum information, albeit not a perfect 
one, as the state recovered by Bob is rarely an exact copy of that sent by Alice. This 
outcome is easily understood if one imagines Alice's quantum state as a wavepacket 
of finite width, initially centred on the first spin of the chain. The wavepacket spreads 
as it travels, hence the information it contains is no longer localized but dispersed 
over several spins. To some extent, this effect is related to the length of the chain: 
shorter chains often perform better as the initial wavepacket has less "room" to 
spread out. As a result of dispersion, when Bob measures the end spin he will only 
recover part of the original wavepacket, consequently receiving an incomplete sample 
of the information Alice sent. 

To optimize the quality of the information he receives. Bob must measure his spin 
at the optimal moment. One might expect this to be the time at which the peak 
of the initial wavepacket reaches Bob's end of the chain; however, this is not always 
the case. Indeed, it is often advantageous to wait a little longer, thus allowing the 
wavepacket to propagate back and forth through the chain several times. As it is 
reflected from the ends, it interferes with itself and may experience a partial revival 
at the site of Bob's spin. If Bob is quick, he can measure at this time, but the 
occurrence of such revivals is unfortunately somewhat difficult to predict. 

The realization of the protocol rests on three assumptions. First, it must be possible 
to initialize the chain to a known state, which in this case corresponds to the ground 
state. In a ferromagnetic chain, this is simply achieved by applying a strong magnetic 
field to polarize the spins, and is therefore not a great challenge. Second, Alice 
and Bob must be able to perform measurements on their spins in order to encode 
and extract the information. This is considerably trickier, owing to the difficulty 
of performing spin-selective measurements on the distance scales involved, though 
possible test systems have been suggested [126] . Finally, these measurements must 
take place on a timescale shorter than the spin- lattice relaxation time Ti, as both 
the encoding and the decoding operations rely on the chain being in a specific state. 
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7.3.1 The Hamiltonian 



The great strength of the protocol we have just described hes in the fact that Ahce 
and Bob do not need to intervene in the transfer process in any way: the state is 
transferred spontaneously by the interaction between the spins in the chain. In this 
respect, the transfer process can be thought of as the solid-state analogue of photon 
transmission through a vacuum. The protocol therefore exploits a natural property 
of the system, as opposed to the vast gamut of quantum computational tasks which 
require the experimentalist to keep a tight rein on any degrees of freedom their 
system may exhibit. One feels almost as if the much-fabled - and so far elusive - 
free lunch were within close reach. 

The characteristics of the transfer process are governed by the Hamiltonian of the 
spin chain. Bose's original proposal sees the spins coupled by a ferromagnetic 
nearest-neighbour exchange interaction and immersed in an external magnetic field, 
such that the overall Hamiltonian of the system reads 



where the first sum runs over all nearest-neighbouring pairs (ij) and the second over 
all N spins in the chain. Both Jij and Bj are chosen to be positive numbers, hence 
the ground state of this Hamiltonian is fully polarized in the positive z direction. 
Excitations from this ground state are represented by de-localized spin-flips, often 
described in terms of spin waves |105j . 

The choice of Hamiltonian is determined by two main factors. First, it is physically 
realistic, as it accurately describes the energy level structure of many ferromagnets 
|155] . Second, it conserves the number of excitations in the system; indeed, if the 
total magnetization is proportional to the number of spin flips, one finds 



TV 




(7.1) 



[Ho, S,] = 0, 



(7.2) 



Chapter 7. The Problem of Scalabihty and the Need for Communication 



151 



where is the sum of the individual spin components in the ^-direction 

N 

As a result of (17. 2p . the dimensions of the Hilbert space needed to to describe the 
system is drastically reduced, as dimTi = (^) , where m is the number of excitations. 
In the simplest case of m = 1, the dimensions of Ti grow linearly with the length of 
the chain. It is therefore possible to solve even very large systems analytically, an 
impossible task if one were obliged to work with the full Hilbert space. 

7.3.2 Notation 

Before proceeding with a more quantitative description of the communication pro- 
tocol, a comment on notation is required. In keeping with the Bloch sphere repre- 
sentation of eigenkets of the Pauli spin operators, it has so far been assumed that 
states |0) and |1) represent the spin up and down eigenstates of a^, respectively. 
The convention adopted in |126j reverses this labeling. Therefore, to avoid possible 
confusion, I will henceforth label the eigenstates of cTz explicitly as | |) and | |). 
The ground state of the chain can then be represented as | iii2i3 ••• Itv), where the 
subscripts refer to the positions of the spins; I abbreviate this state as | JJ-). The ket 
|j) will indicate a state with a spin flip on site j, the ket a state with spin flips 
at sites j and /c, and so on. It will be assumed Alice encodes her information in the 
spin at position s (standing for sender), and Bob chooses to measure the spin at 
position r (standing for receiver). In all but very few cases it will be assumed that 
s = 1 and r = N. 

7.3.3 The Fidelity as a Measure of Success 

With these tools in place, the steps of the protocol can be reiterated as follows: 

1. The chain is initialized by applying a magnetic field and waiting until the 
ground state | is attained. 
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2. Alice removes the first spin of the chain and replaces it with her encoded qubit 

= c^l i) + P\ T) to create the state \ipi) 

m = \a)i\i2... U) = «|^)+/3|l). 

3. The chain is left to evolve for time t to output state IV'/) under the effect of 
dZH): 

\^f) = a\ 4) +/3e-^^o*|l) = a\ 4) + /5E7=i(ie-^''"*|l)|j)- 

4. Bob measures the spin at his end of the chain and obtains the state \b{t)): 

m) = -jL=[a\ i) + /5/(t)|T)], 



with 

f{t) = (iV|e^*^"*|l) 

P{t) = \a\' + mf{t)\'- 

The function f{t) is a propagator, and represents the transition amplitude from 
state |1) to state |N). If the protocol is successful, at some time to Bob will find 
/(to) = 1, hence |fe(to)) = \a). The value of the propagator is therefore related to 
the quality of the transfer, and can be employed to measure of the success of the 
protocol via a quantity known as the fidelity 



(7.4) 



where p and p' are the initial and final density operators. The fidelity evaluates how 
faithfully a quantum channel maps an input state p to an output state p'. It is a 
real- valued function spanning the interval [0,1], which is non-zero provided there 
is some overlap between the input and output states, and unity if p = p'. The 
definition f l7.4p applies to both pure and mixed states, but when dealing with pure 
states it is more convenient to introduce a second quantity, known as the averaged 
fidelity 

\fit)\ , \m\' , 1 



6 2' 



(7.5) 
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where the average is calculated over all pure input states of the Bloch sphere. There 
are several other measures of the performance of a quantum channel, but these are 
more complex and will not be described here. As the present discussion is limited to 
pure states, the term 'fidelity' should hereafter be understood as referring to (17.51) : 
definition (17.41) is not used anywhere in this work. 

If nearest neighbour couplings are constant and the applied field is uniform, it is 
possible to derive simple expressions both for the spectrum of (17.11) and for the 
fidelity. Setting Jij = (-f)^ and Bj = B, we quote the result obtained by Bose 
in [T26] 



N 



\m) = a^^cos ^— (m- l)(2j- 1)J |j) 

Err,. = 2i? + 



1 1 — COS 




— (m — 1) 







(7.6) 
(7.7) 



where m = 1...N, ai = am^ti = y |f) ^-^d the ground state energy Eq = 
_ j{N~i) + 2BN zero. The propagator between sites 1 and of the chain can 



then be expressed as 



N 



/(t) = 5^(i\^m)(m|l)e- 



■iEmt 



(7.8) 



m=l 



To calculate the fidelity, this is numerically evaluated and substituted into equation 



7.3.4 The Performance of the Protocol 

For arbitrary chain length (with the exception of A^ = 2), the evolution of the fidelity 
as a function of time is disordered and largely unpredictable (fig. 17.11) . For A^ < 30, 
with few exceptions, time-optimized fidelities exceed 0.8, but the channel becomes 
progressively worse as the length of the chain increases. Nevertheless, the optimal 
fidelity at A^ = 80 is still greater than that achievable with a classical channel 
|126[ 1127] . The presence of several high 'spikes' in the fidelity is encouraging, as it 
suggests there are several times at which Bob could measure his spin and obtain a 
reasonably good copy of Alice's input state. However, it also implies Bob has a very 
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short window of opportunity in which to make his measurement. For the example 
shown in figure 17.11 for instance, if Bob measures his spin at time ti he will find 
state I a) with over 75% fidelity; if, on the other hand, he measures at the not too 
distant time ^2 he will obtain a much poorer result. 

One might expect the optimal measurement time to coincide with the time to taken 
by the wavepacket to reach the end of the chain. Roughly speaking, this can be 
approximated as the timescale of a single nearest-neighbour coupling multiplied by 
the length of the chain. From the uncertainty principle 

N-1 , , 

t„ « — . (7.9) 

where 1/2 J is the timescale of a single coupling, iV — 1 is the length of a chain of 
spins with unit spacing, and h is set to 1. 

The same result is obtained from equations (17.61) and (17.71) . recalling that a localized 
spin-fiip can be written as a superposition of spin waves with different momenta 

N 

\j) = am COS kra\km), (7.10) 

m=l 

where I have replaced the argument of the cosine with k^, to make explicit the 
momentum dependence of the j-states. Given the dispersion relation of equation 
(17.71) . the group velocity of the individual spin waves 

= ^ = 2Jsmkm, (7.11) 

is itself a function of momentum. It then becomes clear that wavepacket of equation 
(I7.10p spreads because each component of the superposition propagates at a different 
velocity. As observed in [138j . spreading the spin fiip over several sites counters the 
effect by reducing the number of momenta involved in the summation. The fastest- 
propagating component of (17.101) will reach the end of a chain of length L in time 

L N-1 

*» = 27 " — • <"2' 

which is the result of equation (17.91) . 
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The value of to for a chain of 27 spins with J = | is shown in figure 17. 1[ We see 
that, although a local maximum is attained at this time, there are several instances 
with t > to at which the fidelity of state transfer is greater than at to- 




There have been many interesting studies regarding the relationship between the 
optimal fidelity and the length of the chain. It was already noted in [126] that the 
fidelity seemed lower for chains whose length was divisible by 3. Subsequently, it 
was shown that chains of length > 4 cannot achieve unit fidelity |14Uj . and also 
that if is a prime number the fidelity approaches unity in exponential time |156] . 
This may sound enticing, but ultimately one cannot afford to wait so long to achieve 
perfect state transfer. We will return to this point in section 17.61 

7.4 Decoherence 

The discussion so far has been centred on the performance of an isolated spin chain. 
However, to obtain a realistic description of the performance of the protocol in a 
realistic physical setting, interactions of the spin chain with the environment should 
be considered. With few exceptions |157l 11581 1159j . the vast majority of existing 
spin chain protocols neglect the problem of decoherence, and the work presented 
in this thesis is no exception. However, recent studies regarding the influence of 
dissipation and decoherence on the transfer process |160l 1161] allow us to estimate 



Chapter 7. The Problem of Scalabihty and the Need for Communication 



156 



how several of the protocols described above may perform in less idealized settings. 

The decoherence mechanisms most relevant to quantum information processing are 
dephasing and damping [19]. Dephasing takes place when interaction with the en- 
vironment destroys the phase relationships between the components of a quantum 
state. It is a conservative interaction, because no energy is exchanged between the 
system and the environment, and the timescale on which it occurs is usually referred 
to as T2. Damping, on the other hand, affects the populations of the eigenstates 
of the local Hamiltonian. Typical examples might be the relaxation to the ground 
state of a sample prepared in some high-energy configuration, or the appearance of 
thermal spin excitations in a perfectly ordered ferromagnetic state when the tem- 
perature is raised above zero. Damping is therefore a dissipative process, because 
energy is exchanged with the environment. The timescale associated with damping 
is usually referred to as Ti. 

A detailed study of the effect of Ti and T2 processes on spin chain communication 
has recently been carried out in |160j . Here, the authors assume two possible modes 
of interaction of the spin chain with its environment. The 'global' mode involves 
each spin in the chain coupling to a single, 'block' environment, whereas in the 'local' 
mode each spin interacts with a local independent environment. In both cases, effects 
on the fidelity of state transfer are strongly dependent on the internal Hamiltonian 
of the spin chain. For a Heisenberg Hamiltonian, the fidelity in the presence of global 
noise is determined by the transfer time and the so called-environment parameter u, 
which quantifies the relative strength of spin-spin and spin-environment coupling. 
For a chain of given length, the quality of the transfer becomes steadily worse as u 
increases. For a given u, the quality of the transfer is inversely proportional to the 
length of the chain, because the longer the chain, the longer it takes for the state 
to be transferred, hence the longer the coupling time to the environment. Similar 
results are found for the local environment model. 

Interestingly, the situation changes if one considers the mirror-periodic Hamiltonian 
of |139j . In this model, in the absence of decoherence, the time at which perfect 
state transfer is achieved is constant. For a global dephasing or local damping 
environment, this yields an average state transfer fidelity which is independent of 
chain length, but decays with increasing environment parameter. On the other hand, 
in a local dephasing environment, the transfer fidelity falls both with increasing 1/ 
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and with increasing chain length. 

The effect of a pure local damping environment is analogous to that of a spin bath 
at zero temperature [IHl 1162] . Recent work has been dedicated to studying the 
effects of thermal fluctuations in more detail. For short spin chains, the problem 
has been solved analytically |154j . It emerges that finite temperature effects can 
influence the fidelity of state transfer, but not the timescale. Recently, more detailed 
results have been obtained for longer chains globally coupled to a heat bath via 
dissipative exchange interactions |161j . The authors consider state transfer of a 
spinless fermion in a tightly bound lattice, which is analogous to an XY spin model. 
At zero temperature, the fidelity of state transfer is found to decay exponentially 
with time. At fixed temperature, the same exponential decrease is observed; at fixed 
time, however, the performance of the system shows an interesting dependence on 
the coefficients of the state one sends. Specifically, for a state a|0)+/3|l), the transfer 
fidelity decreases with temperature if |a| > \(3\, but increases with temperature if 
|a| < However, this trend is not particularly strong, and becomes negligible at 
large times. 

The above discussion allows us to conclude that, in all realistic systems, decoher- 
ence can seriously impair the efficiency of state transfer. The effect is especially 
pronounced when the transfer time is long, which is increasingly the the 
length of the chain grows. Strategies to combat decoherence might then include 
more robust state-encoding mechanisms or faster transfer times. The latter objec- 
tive could perhaps be achieved by engineering nearest-neighbour couplings. In some 
systems, this might only be possible on a global scale; in others, however, locally 
tunable interactions may be a realistic possibility (see Section 17151) . 

7.5 Experimental Quantum Communication with 
Spin Chains 

In the years immediately following Bose's proposal, research in quantum communi- 
cation with spin chains was largely biased towards the theoretical. Practical achieve- 
ments are now emerging, as progress in instrumentation begins to offer increased 
freedom in the systems we use as qubits. The majority of proposals for experimental 
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implementation of existing communication schemes are centred on superconducting 
qubits. It has been shown that both charge and flux qubits allow high-fidelity quan- 
tum state transfer, even in the presence of manufacturing defects and decoherence 
|163[ 1164] . though no experimental evidence has yet been presented. Quantum state 
mirror-inversion in a chain of three spins has recently been demonstrated in the field 
of NMR, which is well-suited to proof-of-principle trials, but fails for large systems, 
owing to the limitations of the state-preparation procedure |165j . Promising can- 
didates are also 'carbon peapods', single- walled carbon nanotubes containing spin 
qubits trapped in fullerene cages |166j . The properties of the peapod are strongly 
dependent on the nature of the trapped spin. ScQCeo peapods, for example, can 
be accurately modeled as Heisenberg antiferromagnetic spin chains with tunable 
nearest-neighbour couplings, as the spin qubit resides mainly on the surface of the 
fullerene cage |167] . Conversely, in systems such as N@C60, the spin is well-screened. 
Therefore, exchange couplings are more limited, but other mechanisms, such as dipo- 
lar couplings can come into play. In any case, the spin qubit is largely unaffected by 
the presence of the nanotube, which only serves as structural support, and a buffer 
against environment-induced decoherence. Atoms in coupled optical cavities have 
also been considered; recent results show these can simulate anisotropic XXZ or 
isotropic XXX chains of spin qubits with s > | [168]. G iven an optimal cavity, all 
necessary qubit control is provided by laser fields; however, obstacles to practical im- 
plementation include decoherence effects and the requirement of strong atom-cavity 
coupling. Heisenberg and Ising Hamiltonians, as well as more exotic couplings, can 
also be realized in optical lattices [169] . The idea here is that atoms can be made 
to interact controllably by adjusting the trapping potential. Several possible imple- 
mentations have been suggested, however these are also not without their difficulties 
|17U] . Finally, several proposals centred on quantum dot architectures have been ad- 
vanced. These hold much promise, as it is now possible to grow regular two- and 
three- dimensional quantum dot arrays of controllable size and density, potentially 
allowing single-qubit addressability |171] . Interactions between dots are provided 
by the Forster effect, through which an excited donor dot transfers an exciton to an 
un-excited acceptor via an electric dipole coupling |172] . Previous work identified 
the main difficulties with the quantum dot approach in short decoherence times and 
manufacturing defects |173j : however, recent findings suggest these difficulties can 
be largely overcome |174] , and with some luck, practical tests will soon bear out this 
result. 
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7.6 Conclusions 



Quantum communication with spin chains is currently one of the most active ar- 
eas of research in quantum information processing. It has been widely shown that 
an isolated spin chain can afford versatile, high quality distribution of arbitrary 
quantum states with minimal control requirements, though optimal performance 
sometimes requires less flexibility over the operating conditions. A general survey of 
the communication protocols suggested to date indicates that perfect state transfer 
is achievable in a variety of different ways, provided one can engineer the geometry 
and the energetics of the system, and has access to ancillary registers. It is gener- 
ally observed that spin chains allow more robust state transfer than sequences of 
SWAP operations, because the use of a permanently-coupled system minimizes the 
possibility of accumulating errors. The issue of state-transfer time has also been 
raised, because an optimal protocol must allow for both quality and speed. Transfer 
times do not always bear an obvious relationship to the size of the system, though 
in several cases a rising trend is observed. In absolute terms, the transfer time is 
determined by the inter-spin couplings, which in a typical ferromagnet may be of 
order 10 meV. We will see in the following chapter that for the original commu- 
nication protocol of |126j carried out with a chain of 20 spins, this might yield an 
optimal measurement time of order 10~^^ s [this is shown in the following chapter, 
in figure IH^ a)]. Compared to the timescale of a single SWAP gate, which hovers 
about the pico-second mark in a double quantum dot structure |174l 1175] , this fig- 
ure is clearly advantageous; however, other protocols may not perform as efficiently. 
Practical demonstrations of spin chain communication are few, but recent advance- 
ments in experimental methods and instrumentation suggest the situation is poised 
to change. 



Chapter 8 



Quantum State Transfer with 
Long-Range Interactions 



This chapter examines the quality and efficiency of quantum state transfer through 
closed- and open-ended chains of spin-^ fermions coupled by long-range magnetic 
dipole interactions. The two geometries are treated separately because they are not 
physically equivalent; indeed, an open chain has end points, but a closed chain does 
not. The structural difference is reflected in the dynamics, which show interesting 
and perhaps slightly counter-intuitive features. We will see that the fidelity of state 
transfer in a chain has a clear periodicity and a predictable evolution in time, but 
the same cannot be said for a ring. This might seem unusual, given that a ring is a 
periodic structure and a chain is not; however, analysis of the respective eigenspectra 
reveals the origin of this behaviour. After considering a basic communication scheme 
analogous to that described in \12b^ . I will explore possible optimization strategies, 
concluding with some remarks on the robustness of the protocol in the face of exper- 
imental uncertainties, fabrication defects and possible decoherence channels. 



8.1 Introduction 



The majority of spin-chain communication protocols are designed for systems ex- 
hibiting only nearest-neighbour interactions. However, in recent years new effort has 
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been dedicated to studying systems with anisotropic, long-range couplings. The pos- 
sibility of implementing quantum- computational tasks in one- and two-dimensional 
spin lattices with tunable dipolar interactions has been known for some years |176] ; 
however, until the work of Kay |153j . the performance of a spin chain with long-range 
couplings as a communication channel had never been investigated. The study de- 
tailed in this chapter represents the first rigorous attempt to analyze the dynamics 
of quantum state transfer in a chain of spins coupled by the magnetic dipole inter- 
action. Thanks to recent advances in ion- and electron- trapping technology, several 
possible implementations of the present scheme have come to light. These, together 
with subsequent newly published extensions to the scheme, will be discussed in 
section 18.61 

The protocol described here is, in essence, a mapping of Bose's original communi- 
cation scheme to a system in which all the spins are coupled to each other via a 
long-range magnetic interaction that scales as r~^. I examine a simplified system 
in which any external magnetic field is constant and parallel to the axis joining the 
dipoles. This axis is chosen to coincide with the ^-direction, so that the z-component 
of the total magnetization M is conserved (cfr. equation 17. 2p . It is therefore possible 
to work in the single excitation sub-space, which grows linearly in the size of the 
system (see Section [7.3. ip . Within this subspace, the effect of the magnetic field 
is to add a constant to the energies. This constant will hereafter be omitted. For 
simplicity, it will also be assumed that the system is at zero temperature, hence 
contains no thermally excited spin-flips. Unless otherwise stated natural units will 
be used, such that /iq = /is = ^ = me = 1. 



8.2 The General Dipolar Hamiltonian 

The magnetic dipole coupling between two particles separated by a distance VijZ 
can be expressed as follows 

i/., = ^[5,.5,-35^5^], (8.1) 

I "ij I 



where Si and Sj are the total spin operators at sites and r,-, and and Sj are 
their respective z-components. The value of C is determined by the type of particle 
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one considers; for a system of spin- 2 fermions such as electrons, C is given by 



C 



(8.2) 



where Qe is the electron g'-factor. In natural units, the value of C is then C = As 
the majority of the following discussion will be concerned with uniform systems, it 
is useful to define the interaction strength between anti-aligned nearest neighbours 



where a = |?"j.j±i| is the nearest-neighbour spacing. Length, energy and time units 
are hereafter defined by setting a and e to unity. This does not affect the qualitative 
behaviour of the system, as equation (18. ip has an overall scaling factor of ; hence, 
a uniform compression or expansion of the system simply re-scales all energies by 
a constant. Consequently, for a fixed number of component spins, the energy and 
performance of a system of any size can be extrapolated by adjusting the value of a 
as necessary. 

The above considerations are entirely general. The specific form of the energy spec- 
trum and Hamiltonian of the system are determined by the geometry. I will consider 
two cases: the closed and the open chain. These two systems differ in that the latter 
has end points, but the former does not. In other words, the closed chain is a linear 
system obeying periodic boundary conditions, whereas the open chain is a linear 
system not obeying periodic boundary conditions. For brevity, I will refer to the 
closed chain as a ring, however it must be noted this does not in any way assume a 
circular geometry. 
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8.3 Dipolar Rings: the Energy Spectrum and the 
Hamiltonian 



Let us begin by calculating the Hamiltonian of a dipolar ring. From equation (18. II) . 
it is easily found that the off-diagonal matrix elements are 



{j'lm 



\3-3'\ 



•4) 



where e represents unit energy. For each pair {j, j'}, equation (18.41) must be evalu- 
ated in accordance with the minimum image convention, whereby if 



\J-J'\> 



N/2 



for even N 



(N - l)/2 for odd N, 



one sets 



{j'\H\j) = {\N~f\ \H\j). 



1.5) 



1.6) 



The diagonal elements correspond to the average energy of a single-spin-flip state. 
Given the ground-state energies for even and odd 



j=i J \ 

N-l 
^ 1 



it is possible to show that the diagonal matrix elements of the corresponding Hamil- 
tonians are related to the Euler Gamma function T{z) and the Reimann Zeta func- 
tion C{s), as 



{jmjr 
{j\H\jr 



N-4 
2 

iV-4 



Inr(z) +2C(3) 



(8.9) 
(8.10) 



with z = Y and z = ^^j^, for even and odd A^, respectively. Equations (18.91) and 
(I8.10p depend only on A^; therefore, as required by symmetry, the energy needed to 
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flip a spin with respect to the ground state is independent of position along the ring. 

To calculate the eigenspectrum, we recall that the single excitation eigenvectors of 
a periodic system are single-magnon Bloch states. Labeling these states with |m), 
one has 

1 ^ 

|m) = ^ Ve'*^-"^|j), (8.11) 

where km is the magnitude of the Bloch wavevector, r,- is the position of the spin at 
site j, and \j ) is the single spin-flip state. 

If the nearest neighbour separation is set to unity, one has |r, | = j — 1, as the origin 
Vj = coincides with the first site of the ring. Since the excitations travel along the 
dipole axis, it follows that k^, ■ rj = km{j — 1)- The eigenenergies of the ring can 
therefore be expressed as follows 

Em = {m\H\m) = l-J2e"'-^^~^'Hj'\H\j), (8.12) 

j, j ' 

where the sum runs over all pairs j ^ j'. 

According to condition (18. 5p . for a given j there exist y (even A^) or ^ (odd N) 
independent values of \j The sum over j ' contains twice this number of terms, 

because it must account for positive and negative values of (j — j')- It follows that, 
for a given j and even Al] 



2 



J'J' (i-j')=-f 



where we have used the fact that interaction energies of distant neighbours {j, j '} 
are obtained by scaling the nearest-neighbour coupling by a factor \j — j '\^. Finally, 



^The term in j ~ j ' is neglected. 
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considering all possible values of j 

JV 

COS k^{j - j') 



JV 
2 



^^cosKa^ (8.14) 



.=1 ^ 

as e is set to unity. The corresponding expression for odd is obtained by simply 
changing the upper limit of the summation to j = {N — l)/2. 

Equations (17. 5p . (17.81) . (18.111) and (I8.14p are now sufficient to provide a full charac- 
terization of the time-evolution of the fidelity of state transfer in a dipolar ring. 



8.4 The Performance of a Ring 

Thanks to the symmetries of the ring geometry, the fidelity of state transfer can be 
expressed in analytical form. To this end, it is convenient to write the propagator 
of equation (17.81) in more general terms as 

N 

fit)r,s = J](r|m)(m|s)e-^^™*, (8.15) 

m=l 

where | s) and | r) are the sender and receiver sites, respectively. Equation (18. lip can 
be substituted into (18.151) . to obtain 

^ N N 

f{t)r,s = ^ E E {r\j){j'\s)e^'-^^-^\-^^-'. (8.16) 
»^=i j, j '=1 

For a ring of circumference L = Na, Bloch's theorem requires the phase factor 
^ikmNa unity. This yields the allowed values of km'- 



271 _ 2Tin 



(8.17) 
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with a = 1. As the \ j) states form an orthonormal set, the propagator between sites 
r and s is found to be 



N-l 



. 27m(r — s) 



N 

n=0 



with 



N 



^ cos^ 



— (8.19) 



.=1 ^ 



To obtain the fidehty of state transfer, equations fl8.18p and fl8.19p are then substi- 

^ or ^ 
2 2 



tuted into (17.51) . taking (r — s) = ^ ^y^, for even and odd respectively. 



The evolution of the fidelity in time for = 2 — 7 is illustrated in figure 18.11 Figure 
18.21 (a) represents instead the maximum fidelity F^ax that can be achieved within 
an arbitrary time period T = for both a dipolar and a Heisenberg ringj. The 
definition of a cutoff time stems from the practical consideration that one cannot 
afford to wait indefinitely for the state to be transferred. To give a more complete 
description of the efficiency of the protocol, the time r required to achieve fidelity 
Fmax is also shown, in Figure [HI2] (b). 

For A^ > 5, the fidelity is an unpredictable function of time. The maximal fidelity 
falls rapidly with increasing A^, though for even A^ the fidelity is systematically 
higher. A comparison with a nearest-neighbour-coupled ring shows their average 
behaviour is very evenly matched in both quality and speed, although at large 
(odd) A^ a Heisenberg ring has higher fidelity. The results illustrated in figure 18.21 
are complementary to those published in figure 1 of |177j . though not completely 
identical. In the majority of cases, this is owed to the fact that in [177] optimization 
of the fidelity was carried out over the much shorter timescale of However, 
there exist additional numerical inconsistencies between the results presented in 
figure IHTST b) and those of figure 1(b) of |177j . which occur at A^ = 6, 8, 9 and 16 
for the dipolar ring, and A^ = 8, 12, and 24 for the Heisenberg ring. At A^ = 8, 
9, and 12 these arise from the fact that, upon repeating the analysis, it was found 
that the value of the fidelity at the (longer) time of figure 1(b) of [177] exceeds that 
at the (shorter) time of figure IH^ b) by O (10^'^). Given such a small advantage. 



^For a Heisenberg ring, the fidelity is easily obtained from equations ()7.5|) and (jS.lSp by setting 
the eigenenergies = — 2Jcos^^ (see Appendix EJ. 
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it is more efficient overall to measure the output state at this shorter time, hence 
the result presented here. The discrepancies at = 6, 16, and 24 are unfortunately 
owed to a recently discovered indexing oversight, but this does not yield qualitative 
disagreement with previous findings except in the case = 16. Fortunately, this is 
not sufficient to affect the conclusions drawn so far, both in |177j and in the present 
thesis. 



N=2 



N=3 



1.0F 

0.9 
0.8 

0.7 



20 40 60 80 100 



N=4 




0.80 
0.75 
0.70 
0.65 
0.60 
0.55 
0.50 



I! 11 





20 40 60 80 100 

t 

N=5 



0.75 






0.70 






e 0.65 

Li_ 






0.60 






: 0.55 






- 0.50 


- 





20 40 60 80 100 

t 

N=7 



20 40 60 80 100 

t 



Figure 8.1: The evolution of the fidehty at short times t < — for N 
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Figure 8.2: The maximal fidelity (a) and the timcscale (b) of state transfer between diametrically 
opposite sites for a dipole-couplcd (thick red dashed curve) and nearest-neighbour coupled (blue 
dotted curve) ring. The fidelity is optimized over a time period T = i^^. In SI units, the value 
of T depends on the inter-spin distance one chooses. For a chain of particles spaced at 3 A, say, 
the dipolar interaction between nearest neighbours is approximately 4/xeV, roughly corresponding 
to a timescale of 10~^° s. The chosen cutoff time is then of order 100 ns, which is well within 
the decoherence time of, for example, NQCeo |1161 1117] , The black dot-dashed line indicates the 
maximum fidelity for classical transmission of a quantum state. 



8.5 State Transfer in a Dipolar Chain 

The following section is dedicated to the analysis of quantum state transfer in a 
dipolar chain. The lack of periodicity resulting from the finite structure yields 
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important qualitative differences between the chain and ring geometries, which from 
many viewpoints render the chain the more efficient of the two quantum channels. 
As we will see, the main drawback of the chain configuration is the length of time 
required for a quantum state to be optimally transferred, which remains essentially 
impervious to most optimization techniques suggested so far. However, drawing also 
on more recent work published earlier this year |178j . the performance of the dipolar 
chain will be critically evaluated, to establish whether such a system will ever be 
useful for anything other than proof-of-principle applications. 

8.5.1 The Hamiltonian 

It is instructive to begin our analysis by calculating the Hamiltonian of the uniform 
dipolar chain. It is easily shown that the off-diagonal matrix elements are given 
by equation (18. 4p . but this time without recourse to condition (18.51) . The quantity 
\j — j'\ is therefore the linear distance between sites j and j', irrespective of the 
values of these indices. The diagonal elements are calculated by evaluating the 
energy gap between the ground state and the spin- flipped state \j ) . Given the 
ground-state energy 



Equation (I8.2ip can be compared to (18.91) and (18.101) to highlight an extremely 
important characteristic of the chain geometry: the energy required to flip a spin 
from the ground-state configuration is dependent on the site of the spin flip, and 
reaches a minimum at sites 1 and (Fig. 18.31) . The first and final sites of the chain 
are therefore the most favourable environments for a spin flip to exist, which bodes 
extremely well for the prospect of transferring a spin flip between these sites. The 
full implications of this result will be illustrated in greater detail in section I8.5.2[ 




(8.20) 



it follows that 




(8.21) 
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Figure 8.3: The on-site energy {j \H\j) as a function of j for a dipolar chain of 15 spins. 

8.5.2 The Performance of a Chain 

Equations (17.51) . (17. 8p . (18.41) and (18.211) provide all the information required to cal- 
culate the fidelity of state transfer in a dipolar chain. The time evolution of the 
fidelity of state transfer between end sites is shown in figure 18.51 for N = 2 — 7. 

It is immediately evident that, provided a large enough sampling time is chosen, the 
behaviour of a dipolar chain is substantially different from that of both a ring and 
a Heisenberg chain. The first thing we observe is the presence of a regular - though 
slightly noisy - oscillation, whose amplitude is close to unity. This shows there exist 
predictable times at which the transfer process is close to being perfect. Second, the 
period of the oscillation is long and increases with A^, which suggests transfer time 
grows with the length of the chain. In some respects, this is a disadvantage; in others, 
one could argue that a slow variation of the fidelity in time provides the receiving 
party with a reasonably large window in which to measure their spin. Indeed, the 
fidelity remains at least above the classical threshold for roughly a quarter of a 
period either side of the peak. 

The period of oscillation of F{t) is uniquely defined by the energy splitting AA 
between the two lowest eigenvalues of the Hamiltonian. The transfer process is 
therefore dominated by the beating of two nearly degenerate states localized near 
the ends of the chain. As noted in the previous section, this behaviour is explained 
by the variation of the on-site energies of the spins as a function of j; it is clear 
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from fig. 18.31 tliat the most favourable positions for a spin to flip are sites 1 and 
N. Consequently, states |1) and |A^) are the most strongly coupled to the system's 
(two) bound states, which are shown in fig. 18. 4[ In a dipolar chain, this phenomenon 
is a natural consequence of the geometry, but systems in which the spin flip energy 
is specifically chosen site by site have also been studied |179j . 
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Figure 8.4: The energy splitting AEm between the ground state | ij- ) and states with a single 
flipped spin, which shows the evolution of the two bound states. Values of AEm of the same index 
m, counting from the bottom of the spectrum, are shown in the same colour. 



The period of F{t) is related to AA by 



. - ^ (8.22) 

Consequently, for each the time at which F(t) first peaks is 

t,,(N) = I = ^. (8,23) 
This time grows with chain length, as the splitting AA decreases with increasing N. 
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Figure 8.5: The evolution of the fidelity for = 2 — 7. The range of abscissa values varies with TV 
in order to highlight the characteristic oscillation of the fidelity with time. 

Fig. 18.61 compares the performance of dipolar and nearest-neighbour-coupled chains. 
We note that, in addition to = 2, = 3 and = 4 also give perfect transfer, 
and in general F^ax > 0.9. This is a marked improvement on a Heisenberg chain; 
specifically, it is no longer the case that the fidelity is poor when is a multiple of 3 
|126] . These high fidelities are achieved at the cost of long transfer times; indeed, at 
large A^, tQ{N) is proportional to the cube of the chain length (fig. 18.61) . One is forced 
to conclude that, unless optimization strategies can be found, the time required to 
complete the protocol will become impractical as the chain length increases. 
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Figure 8.6: The maximal fidelity (a) and the timescale (b) of state transfer between sites 1 and N 
of a dipole-coupled (thick red dashed curve) and nearest-neighbour coupled (blue dotted curve) 
chain. The black dot-dashed line in (a) indicates the maximum fidelity for classical transmission 
of a quantum state. Fig. (b) shows the times of optimal fidelity on a logj^Q scale for both types 
of chain. Here, the black dot-dashed line represents the function y = {N — 1)^. At large N, this 
is parallel to the red dotted curve, indicating that the transfer time in a dipolar chain scales as 
the cube of the chain length. In absolute terms, for a chain of particles spaced by 3 A, say, unit 
time corresponds to 0(10"^°) s. Transfer times for chains longer than 20 spins therefore exceed 
the micro-second mark. 



Following the work of Kay |153] and Osborne and Linden |138j . two main optimiza- 
tion methods will be investigated. The first consists of changing the encoding and 
measurement convention so that the sending and receiving parties have access to 
more than one spin at the ends of the chain. I refer to this as input and measure- 
ment optimization. The second involves changing the geometric properties of the 
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chain, such as the number of spins or their placement. I refer to this as structural 
optimization. 

8.5.3 Input and Measurement Optimization 

Let us begin by asking whether the performance of the protocol can be boosted by 
allowing Alice and Bob manipulate the first and last two spins of the chain. For 
example, suppose Alice wishes to transmit the state |1). Rather than encoding this 
state in the first spin only, Alice encodes it in the first two spins, generating the 
initial state 

|s) = ci|i) + C2|^). (8.24) 

Following the discussion of the previous section, it is advantageous to choose the 
coefficients Ci and C2 so as to maximize the overlap of \s) with one of the system's 
bound states, as these dominate the transfer. The eigenvectors of the system are 
symmetrized and antisymmetrized combinations of |j)-states, having the form |m) = 
aj^lj), with |a7^| = |a^_]^^j|, and G For either of the bound states (m = 1 
and m = 2), the amphtudes {a^ ■ ■ ■ a^„2} decrease steadily with N; hence, most of 
the amplitude is concentrated on the first and last two sites. Suppose Alice takes 
the first of these bound states as her 'reference' state; she should then choose the 
coefficients ci and C2 to be as close as possible to a\ and a^, while still obeying the 
normalization condition + c\ = 1. 

Bob must now establish which output state |r) to aim for. Knowing as he does 
the criterion Alice has adopted to fix her ci and C2, and also that la™] = la^v-i+jl' 
he chooses |r) = C2\N-1) + ci\N). The propagator of equation (17.81) can then be 
expressed as 

N 

fit) = Yl + (C2a:^_i + cia^) e'^-\ (8.25) 

m=l 

The outcome of this optimization procedure for a chain of 7 spins is shown in figure 
18. 7[ Comparing with figure 18. 5^ it is evident that the transfer time has remained 
unaffected. However, we have been successful in smoothing out the oscillation and 
boosting its amplitude even closer to unity. This agrees completely with the result 
of |138j : indeed, the quality of the transfer can be further improved if Alice and Bob 
are have access to an even larger portion of the chain, as shown in figure [821(b) . 
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Figure 8.7: The evolution of the fidelity for N ~7, when the initial and final states are distributed 
over (a) the first and last two, or (b) the first and last three spins of the chain. 



It is also interesting to ask how well the system performs if the initial and final states 
are chosen at random; this is illustrated in figure We find that the characteristic 
oscillation of F{t) is lost, unless either |s) = \2) and |r) = |iV), or |s) = |1) and 
\r) = \N - 1) [figure 18.8( b)]. However, in both cases the signal is considerably 
noisier, and the maximum fidelity is greatly reduced. This is a result of the lesser 
efficiency of coupling to the bound states as one moves away from the ends of the 
chain (cfr. fig. 18.41) . 

Conversely, fidelities in excess of 0.8 can be attained if sites s and r are symmetrijfl 
(Figure I8.8p . In this case also the oscillatory pattern is lost, as the eigenstates 
that dominate the transfer are well-separated in energy. In return, however, for 
< 23 the fidelity can peak at times shorter than r/2. The system is therefore 
more efficient, though less predictable. Furthermore, the choice of symmetric sites 
implies Alice and Bob must be capable of site-specific encoding and measurement 
at arbitrary positions in the chain; in other words, they must be able to act upon 
a more or less extended portion of the chain containing sites s (Alice) or r (Bob). 
From a practical viewpoint, this is much more demanding than end control only. 

We conclude that by optimizing the encoding and decoding operations whilst re- 
maining in control of the end spins only, Alice and Bob can ensure perfect trans- 
mission of quantum information. However, the problem of excessively long transfer 
times still remains. This is not surprising, as the transfer time is determined by the 

■^A mathematical justification of this result was recently published in |178| . 
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Hamiltonian rather than the details of the encoding. It is therefore clear that to 
act on the transfer time one must change the Hamiltonian, and as the Hamiltonian 
is proportional to the spacing of the spins, it can only be modified by changing the 
structural properties of the chain. 



(a): s=3, r=9 (b): s=1, r=14 




400 800 



Figure 8.8: The evolution of the fidehty for = 15, when the initial and final states are (a) chosen 
at random, (b) taken as \s) = |1) and |r) = |14), and (c) symmetric (|s) = |3) and |r) = \13)). We 
note that in (c) the fidelity peaks within the time limit T = imposed for the ring geometry. 



8.5.4 Structural Optimization: Non-Uniform Chains 

It has already been suggested in |153j that non-uniform chains may offer some ad- 
vantages over uniform systems. To continue this line of investigation, we study a 
chain in which the placement of the spins is uneven, though not random. There are 
several ways to obtain such a structure; here, we imagine compressing or expanding 
some part of a uniform chain by a fixed amount. The system can then be described 
in terms of four parameters: (a) the number of spins in the chain, N; (b) the number 
of evenly spaced spins before the compressed (expanded) region, i; (c) the number 
of evenly spaced spins after the compressed (expanded) region, /; (d) the spacing 
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of the spins in the compressed (expanded) region, 6. The Hamiltonian remains that 
of equation (18.11) . except the overall scaling factor of is now lost. 

Figure IHlQl shows the performance of a chain of 10 spins for different {N,i,f,S} 
combinations with i f. It has been assumed that the 'default' spin spacing is a, 
hence 6 is expressed in units of a. We observe that when the chain is not mirror- 
symmetric the fidelity is poor irrespective of the choice of parameters. Conversely, 
if the chain has mirror symmetry [i = f) it is again possible to achieve high-quality 
state transfer between all symmetric sites, as illustrated in figure [KTUl The times of 
peak fidelity vary considerably; however, for a given s and r, these times are shorter 
if the central section of the chain has been compressed (i.e. 6 < 1). This is not 
surprising, as compressing the central section shortens the chain. 

If the chain is non uniform but mirror-symmetric, the time evolution of the fidelity 
of state transfer between the end sites exhibits a regular oscillatory behaviour, which 
is different from that of a uniform chain. The details of the oscillation depend on 
whether the chain has been expanded or compressed in the central region. If the 
chain has been expanded, the transfer is no longer dominated by the bound states 
but by the four upper-lying eigenstates of the system. This is a consequence of the 
geometry, as the end spins are now closer together (hence more strongly coupled) 
relative to the spins in the mid-section of the chain. As a result, there are now at 
least three clearly discernable frequencies beating against each other. The first is 
the frequency of the slowly varying envelope function, which is determined by the 
splitting of the third and fourth highest eigenenergies of the system. The second 
characterizes the rapid oscillation within the envelope, and is determined by the 
energy splitting of the highest and second-highest eigenstates. In turn, this serves 
as an envelope for an even faster oscillation, whose period is defined by the energy 
gap of the highest and third-highest eigenstates. These details are shown in figure 
18.111 The peak fidelity is again close to unity, but the transfer now takes even longer 
than in a uniform chain, because expanding the central section of the chain increases 
its length. 

Similar results are found for compressed chains with the exception that, as in uniform 
chains, we are concerned with the bottom of the energy spectrum rather than the 
top. The frequency of the outer envelope is therefore determined by the energy 
splitting of the two lowest-lying eigenstates, the frequency within this envelope is 
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determined by the energy splitting of the third and fourth lowest-lying eigenstates, 
and so on. The transfer is quicker than in a uniform chain, as compressing the centre 
of the chain makes it shorter. 
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Figure 8.9; The fidelity of state transfer between sites 1 and 10 of a ten-spin chain. Values of 
{N, i, /, 6} are indicated on each figure, with the corresponding chain structure inset. 
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Figure 8.10: The fidelity of state transfer between symmetric sites of a non-imiform ten-spin chain. 
The Icft-liand figures take s — 3 and r ~ 8, the right-hand figures take s = 4 and r — 7. Values of 
{N, i, /, i5} are indicated on each figure, with the corresponding chain structure inset. 
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Figure 8.11: The fidelity of state transfer between sites 1 and 10 of an expanded ten-spin chain, 
illustrating the different frequencies involved in the transfer. The plot to the right of each figure 
is a 'zoom-in' of the boxed region. 
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Figure 8.12: The fidelity of state transfer between sites 1 and 10 of a compressed ten-spin chain, 
illustrating the different frequencies involved in the transfer. The plot to the right of each figure 
is a 'zoom-in' of the boxed region. 

In conclusion, we find that expanding or compressing the chain while keeping the 
number of spins fixed can offer some degree of speed-up over a uniform geometry 
provided the transfer takes place between symmetric spins not located at the ends. 
However, this comes at the price of losing the regular oscillation which allows us to 
predict when the fidelity will peak. It would therefore be more desirable to find a 
way to minimize the transfer time while preserving the fidelity's periodic behaviour. 
This brings us to our last attempt at optimization, which will investigate whether 
the performance of a chain of a given length depends on the number of component 
spins. 

8.5.5 Structural Optimization: Fixed length Chains 

In order to analyze the relationship between the efficiency of state transfer and the 
number of spins for a fixed chain length L, we define t^ as the transmission time 
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giving maximum fidelity at unit chain length. Hence 

*; = (8.26) 

with L = a{N — 1). Fig. 18.131 shows the development of with increasing A^, and 
reveals two interesting features. The first is the asymptotic behaviour of at large 
TV; this suggests there exists a threshold length above which the central part of the 
chain ceases to play a role in the transfer. In other words, the evolution of the 
system is determined almost exclusively by the coupling between the first and last 
i spins, irrespective of the number of spins that separate them. 
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Figure 8.13: The behaviour of as a function of the number of spins in the chain. Note the 
minimum at = 4 and the flatness of the curve for > 15. Units are as specified in Section [5TT1 



To explore this hypothesis, and determine the behaviour of tg for large A^, we work 
with states \B) and \E) localized on the first and last i spins of the chain, respectively. 
These are the bound-state eigenf unctions of a semi-infinite chain extending to the 
right and the left, and can be expressed as follows 

i 

\B) = J]a„|n), (8.27) 

n=l 

i 

\E) = ^a„|iV + 1 - n). (8.28) 

n=l 

Let us now assume the coefficients are non-zero only if < i oi > 

N + 1 — i. The states \B) and \E) can then be used as an alternative basis to 



Chapter 8. Quantum State Transfer with Long-Range Interactions 



183 



represent the Hamiltonian as a 2 x 2 matrix, which takes the form 



n 



{B\H\B) {B\H\E) 
{B\H\E) {B\H\B) 



since by definition {B\H\B) = {E\H\E) and {B\H\E) = {E\H\B). The difference 
between the two bound-state energies of (18.1 p therefore corresponds to the difference 
between the eigenvalues of Ti. This is given by 



From (18.41) . we have 



Hence 



AA = 2{B\n\E). 



{j'\H\j) = H{\j-j'\). 



(8.29) 



(8.30) 



{B\n\E) = J2 



\H\N + 1 - m) = ^ alamH{\N + 1 - m - n|) (8.31) 



n,m=l 



n,m=l 



We adopt a 'dummy' variable X = |j — j '|, so that 

C 



H{X) 



It follows that 



2a3X3 
dH{X) 3C 



(8.32) 



(8.33) 



dX 2a^X^' 

Using (I8.3ip and the fact that L = a{N — 1), we can expand H{\N + 1 — m — n\) as 
a Taylor series to first order inx = m + n — 2. Then 



H{\L-x\) = H{L)-x^ 



C ^ 3Ca(m + n - 2) 



X=L 



2L3 



2L4 



(8.34) 



Therefore 



{B\n\E) 



~2 

C 
~2 



^ i i 



n,m=l 

Q aR 

1?^U 



n,m=l 



(8.35) 
(8.36) 
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with 

i 

Q = ^ (^*n(^m, (8.37) 
n,m=l 
i 

R = 3<a^(m + n-2). (8.38) 

n,m=l 
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Figure 8.14: Comparison between the predictions of a 14-spin model (blue points) and the data 
calculated from the treatment of the system in its entirety (red points). Note the model is only 
accurate at > 13. 



It emerges that the asymptotic behaviour of the system can be modeled very ac- 
curately if the coefficients taken to be the amplitudes of the ground state 
eigenvector of the i x i sub-matrix of the original dipolar Hamiltonian. The values 
of an used in the fit of figure 18.141 are thus obtained from the dipolar Hamiltonian 
of a 14-spin chain, with i = 4. The corresponding values of Q and R are Q ~ 0.325 
and R ~ —0.957. These quantities show only a very weak dependence on N, so they 
have been treated as constants. The fact that Q < 1 indicates that, at large A^, 
the transfer rate is always less than that attained between two completely isolated 
spins; equation 08.371) indicates this is a result of interference between positive and 
negative components in the localized states \B) and \E). 

At small A^, however, this is clearly not the shown by the minimum in 

fig. 18.131 at A^ = 4. The concavity of the curve between N = 2 and N = 7 indicates 
that only chains with less than 7 spins can improve on the performance of a simple 
dipole pair; specifically, in a uniform chain of a given length, the best compromise 
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between the quality and the speed of the communication is obtained with 4 spins. 
This occurs because for small the bound states at the ends of the chain have 
a large overlap, i.e. there exist terms in eqn. (18.31 1) which simultaneously have 
significant positive values of a*am and small values of|A^ + l — m — ?t,|. 

The uniform 4-spin chain can be optimized still further by symmetrically moving the 
inner spins slightly closer to the ends, so that ri^2 = 1^3,4 ~ 0.312L and r2,3 ~ 0.375L. 
The resulting transfer rate is then ^ 0.512, which is a 10% improvement on the 
value one finds for a uniform system, ^ 0.568. Considering the extra difficulty 
involved in such high precision placement of the spins, it is probably not useful to 
adopt this optimization technique as a rule. However, viewed from another angle, 
our result also indicates that if a symmetric error were made in the placement of 
the central spins, the protocol would not necessarily be compromised. 

To explore in more detail the robustness of the protocol against structural imper- 
fections in the chain, we can study the effect of introducing a random error on the 
placement of each spin. To evaluate whether there is a significant drop in perfor- 
mance, I have counted over 1000 iterations the number of times the fidelity of a 4-spin 
chain of arbitrary length, measured at the optimal time to = Lt^, falls below 2/3. 
It has been assumed the magnitude of the error lies in the interval [— ^ < A < ^], 
with p G [0, 1] so that the maximum allowed deviation of a spin is the distance to the 
next (or previous) site. As one might expect, the performance of the system is found 
to depend on the magnitude of A. If A < O.OIL, the fidelity at time to rarely drops 
below the classical threshold. However, as A increases the number of instances of 
poor performance grow, and by the time A ^ 0.03L the system is failing more than 
50% of the time. Let us now be more demanding and raise the lower limit of the 
fidelity to 0.8; in this case, an error of A < O.OIL yields a 25% failure rate, and 
at A ^ 0.07L the system can essentially be written off. These statistics may seem 
discouraging, but we must remember that as L becomes large the absolute value of 
A also increases. Hence, in terms of precision requirements, long chains are more 
versatile than short ones because they allow more leeway for experimental error in 
the placement of the spins. 
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8.6 Extensions and Experimental Implementations 

In section 17.51 several methods for implementing quantum state transfer with short 
range interactions were reviewed. Among them, endohedral fullerenes are the leading 
candidates for the protocol discussed here, though natural magnetic dipolar systems, 
such as LiHoF4, should also be taken into account [180] . Trapped ions are currently 
emerging as ideal subjects for one- and two-dimensional quantum manipulations. 
Indeed, as widely demonstrated by a variety of impressive experiments [36l 11811 11821 
11831 1184] , ions can be efficiently cooled, initialized, entangled, read out and, most 
importantly in the present context, stored at fixed positions. Ion microtraps allow 
the highest degree of control over the qubits' placement, as each ion is individually 
confined. Linear traps can be used to create linear arrays, whereas confinement in 
two dimensions is achieved with Penning traps. In these systems, however, the single 
ionic positions cannot be fixed to the same degree of precision. In all three cases, 
the trapping potential is provided by a series of orthogonal electric and magnetic 
fields, whereas the inter-qubit interactions are mediated by common vibrational 
modes, usually tuned with directed laser pulses. It has recently been shown that by 
adjusting the trapping potential and the laser frequencies, it is possible to simulate a 
dipolar Hamiltonian similar to equation (18. ip . This involves accessing the so-called 
'stiff' vibrational modes, in which the Coulomb repulsion of the trapped qubits can 
be considered a perturbation on the trapping potential. 

Motivated by these findings, new effort is being dedicated to optimizing quantum 
state transmission through dipolar chains |178] . An interesting proposal, advanced 
earlier this year, involves a system known as the 'double hole (DH) chain'. This is 
a linear chain of uniformly spaced, magnetic-dipole coupled spins, from which the 
particles at sites 2 and — 1 have been removed, isolating the end spins from the 
middle of the chain. As a result, the transfer process is faster than in a 'complete' 
chain, and unit fidelities are achieved even in chains of over 80 spins. One can 
show that, by continuing to symmetrically remove spins, the performance of the DH 
chain can be further improved; however, no matter how many times the process is 
repeated, for a given chain length the lower bound for the transfer time is still that 
provided by a uniform four-spin system (figure [S.lSp . Nevertheless, a 'many hole' 
chain may have other advantages; indeed, it is argued in |178] that by splitting off 
the sending and receiving spins from the mid-section of the chain, perturbations in 
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the channel are less likely to impair the protocol, making the DH chain more robust 
than its complete counterpart against possible decoherence. 




Figure 8.15: The evolution of the fidehty of state transfer over a distance L = 9 in (a) a complete, 
double-hole and four-hole chain (blue, red and cyan curves respectively), and (b) a six- hole chain, 
a uniform four-spin system, and two single spins (green, pink and orange curves, respectively). 
We note that a uniform chain of length L = 9 is composed of ten spins, therefore a six-hole chain 
corresponds to a four-spin configuration, in which the two central spins have unit spacing and are 
separated from the end spins by a distance L = 4. It is evident from figure (b) that the minimum 
transfer time is achieved in a uniform four-spin chain. 

A dipolar Hamiltonian can also be simulated in a solid-state spin system, using 
nuclear magnetic resonance |185l 1186] . Most recent work in this area relates to 
the appearance of multiple correlations between different parts of the spin chain, 
which can generate entanglement between distant spins. Interestingly, it is observed 
that long-distance entanglement is more likely to persist if the spins are coupled by 
nearest-neighbour interactions only. 

8.7 Conclusions 

Unmodulated spin chains with long-range dipolar interactions allow high quality 
transfer of quantum information. Periodic systems do not perform as well as finite 
chains, particularly if the component number of spins is even. Nevertheless, for 
systems of up to 30 spins, a dipolar ring often outdoes its nearest-neighbour coupled 
counterpart. In general, the fidelity of the channel is a decreasing function of chain 
length. If the chain is periodic, the fidelity falls quite rapidly with A^; however. 
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in a finite chain the optimal fidehty remains well above the classical threshold for 
chains of up to 100 spins |178j . Furthermore, in a linear system the fidelity peaks at 
predictable times, and its dynamics are such to ensure the receiver has a reasonably 
ample window of opportunity in which to retrieve their spin. Both these factors 
are owed to the fact that the transfer process is dominated by the two end spins. 
The presence of the mid-section of the chain registers merely as noise; indeed, the 
more one separates the end spins from the mid-section, the better the fidelity. For 
a uniform chain of given length, the number of component spins is relevant to the 
overall performance of the system, which peaks when = 4. 

The main drawback of the linear chain geometry is the time required for the state 
to be transferred, which grows as the cube of the chain length. Recent work has 
managed to reduce this figure by a factor of three, however in absolute terms this 
is not a sizeable advantage, because the transfer time in chains longer than L = 20 
still remains of micro-second order. Therefore, if all factors are taken into account, 
a dipolar chain is a less efficient communication channel than a nearest-neighbour 
coupled chain. However, evaluating the effect of dipolar couplings on quantum state 
transfer is a valuable exercise, because any practical implementation of spin-chain 
communication would most likely be carried out in systems exhibiting both types 
of couplings, though perhaps not in equal measure. One could therefore imagine 
engineering the relative strengths of the nearest-neighbour and dipolar interactions 
so as to generate constructive interference, thus yielding a fidelity higher than that 
achievable by either effect taken separately. 
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Summary and Future Work 



The past 50 years have witnessed the evolution of Quantum Information Technology 
from a theoretical construct to a tangible reality. Most recently, the field has been 
developing at a truly impressive rate, and seems on the verge of a leap of quality 
despite several outstanding issues, mainly related to preserving quantum properties 
for long enough to perform useful tasks. This thesis has explored means to achieve 
the extremely desirable goals of quantum communication over short distances in 
solid-state devices, and experimental entanglement generation. Novel techniques to 
attain these goals have been suggested, and related practical difficulties have been 
discussed. 

Chapter [2] was dedicated to an overview of quantum entanglement and its basic 
properties. Methods to quantify and measure entanglement were also analyzed. 

Chapter [3] put forward a proposal for a scattering experiment designed to cre- 
ate entanglement between distinct neutrons. The protocol was illustrated in detail, 
together with the approximations underlying the mathematical model of the ex- 
periment. Given the possibility of describing the system in terms of pure states, 
it was decided that the scattering event should be treated using a time-dependent 
quantum mechanical analysis. In Chapter HI it was verified that if no momentum 
is transferred between the incident neutron and the sample, the neutron interaction 
time can be identified with its free-flight time through the sample. 

In Chapter [5l the performance of the scattering protocol at zero momentum trans- 
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fer was assessed. It was found that the first pair of neutrons to scatter from a sample 
having the properties of a ferromagnetic insulator could come to share a substantial 
degree of measurable entanglement, subject to the possibility of applying a static 
external magnetic field of tuneable strength. The experimental feasibility of the 
protocol was discussed, and it was concluded that the only true obstacle to a practi- 
cal realization stemmed from the coherence properties of the neutron beam. These 
properties depend only on beam preparation apparatus and technique. Therefore, 
the protocol may indeed become feasible as neutron facilities become better equipped 
and more refined. 

Chapter [6] extended the scope of the entangling proposal to account for the possi- 
bility of many neutrons scattering from the sample. The aim of the chapter was to 
establish whether a measurement on any pair of neutrons detected within a certain 
time of each other would yield an entangled state. It was shown that the average 
entanglement of any pair is at best of order 0.1 ebits. This falls with the size of the 
sampling ensemble of outgoing neutrons. Consequently, to produce highly entangled 
states it is necessary to detect one of the very first scattered pairs. The emerging 
trends are interesting from a more fundamental viewpoint also, as they provide a 
practical means to establish whether or not two distinct neutrons can in fact be 
entangled. 

In Chapter [71 the focus shifted from entanglement generation to quantum informa- 
tion transfer in the solid state. The chapter was dedicated to a review of quantum 
communication through spin chains, and a detailed analysis of the pioneering pro- 
posal of Bose [126] . Recent experimental implementations of this proposal were also 
discussed. 

Finally, Chapter [H] analyzed the performance of dipolar chains as conduits of quan- 
tum information. It was found that near-perfect transfer can be achieved without 
any special engineering of the system, but at the cost of long transfer times. Possible 
optimization methods were explored, and it was concluded that for a given distance 
the optimal balance between quality and speed was obtained with a uniform chain 
of four spins. Recent theoretical and experimental work on dipolar systems was also 
reported. 

Much of this work has been concerned with simulating possible experimental con- 
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ditions. As with many theoretical models, this leaves a great deal of scope for 
refining underlying approximations to attain a more faithful description of reality. 
Consequently, there are several directions in which future work might proceed. 

In the context of neutron entanglement, a first question to ask might be how the 
protocol is affected when the temperature is raised above zero. Simple preliminary 
calculations suggest that if scattering takes place from a thermal state of the sample 
the outcome of the protocol is determined by the relative magnitude of the thermal 
energy to all other energy scales in the problem. Due to heavy computational 
demands, only small systems have been considered so far. It would be interesting to 
look at larger systems, and in greater detail. One could also examine the effect of 
changing the properties of the sample, and in this respect the possibilities are clearly 
boundless. Samples with different electrical properties (metals, semi-conductors, 
superconductors...) or more realistic anisotropies could be considered. 

As regards quantum state transfer through dipolar spin chains, it may be interesting 
to re-asses the performance of the protocol in the presence of mixed short- and long- 
range interactions. A slightly more ambitious goal might be to scale up from a single 
spin chain, and study the flow of information around some form of 'dipolar quantum 
network' . 

The problem of time in scattering theory also offers huge scope for further investiga- 
tion. From a very brief survey of the field, it would seem the matter has been 'on the 
table' for the past fifty years at least, and still no unique resolution exists. As com- 
mented in chapter m this may be due to a dependence on the context. Indeed, much 
of the literature tends to explore very specific scenarios. One could therefore begin 
with an extensive survey of work done to date, to identify results of general validity 
and gain insight on the most efficient way to tackle the problem. Some research into 
the dynamical theory of neutron scattering may also prove enlightening. 

As a final word, I would like to say that the past three years have been an 
incredible experience, which has taught me a great many things about myself in 
both a professional and a personal capacity. I have often felt the only thing that 
could have prepared me to do a PhD would have been... to do a PhD! Paraphrasing 
slightly the words of Cervantes' legendary Don Quixote: I have struggled, I have 
made mistakes, but I have done this work as best I can according to the world as I 
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see it. 

One of the most precious things I have learnt is what it takes for me to truly under- 
stand a concept and make it mine. I discovered this rather late in my PhD - perhaps 
I was not in a position to see it sooner, but that is another story. What matters is 
that I have rediscovered the joy of questioning, and the freedom to question what I 
wish. I hope to carry that feeling of wonder and adventure into any work I may do 
in the future, conveying it to those around me always. Thank you for listening. 



Appendix A 

The 5-Matrix Method 



The S'-matrix formalism was originally proposed by Heisenberg in the early 1940s, 
as an attempt to reformulate quantum field theory in terms of observable quantities 
[187] . This project, which became known as the S'-matrix program, was abandoned 
soon after, but laid the foundations for the future developments of scattering theory 
at the hand of, among others, Gell-Mann and Goldberger |109] . and Lippmann 
and Schwinger. The treatment I present here is that of |lU9j . which describes the 
emergence of an operator akin to the Heisenberg S'-matrix from the interaction 
representation. 

Consider a scattering system represented by a Hamiltonian H = K + V , where 
K represents the kinetic energy of the system, and V represents the scattering 
potential. Let us denote with $j (t) = 0je~*'^'* the stationary state solutions to the 
Schrodinger equation (SE) in the absence of the potential, where the (pi are solutions 
of the time-independent SE. The full SE takes eigenstates \E'(t), such that 

i^-^ = {K + V)^{t). (A.l) 

We work in the approximation that interaction timescales are negligible with respect 
to measurement timescales. Therefore, the initial and measured states of the system 
coincide with eigenstates of K. Suppose, then, we wished to calculate the cross- 
section for scattering from initial state $j to final state $j at some convenient time 
t = 0. To describe the history of the system up until the moment t = 0, one can 
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represent the initial state ^1/^(0) as a superposition of $j states, i.e. as a 'train' of 
waves released at some time T in the distant past and fed into the sample over a 
prolonged period of time. Hence 

where r = will eventually be taken to the limit of +00, and the factor of e^^ is 
included to ensure the integral converges in this limit. 

Let us now define an interaction representation wavefunction by 

qj'(t) = e*^**(t), (A.3) 
such that ^'(0) = ^(0). The SE for ^'(t) reads 



dt 

with 



V{t)^'{t), (A.4) 



(t) = e*^ Ve"'^*. (A.5) 

can be written in terms of \l/'(to) at some time to with the help of the time- 
evolution operator: 

^\t) = U{t,to)^\to). (A.6) 



Knowing that = e *-'^*\l/(to), equations flA.Sp and flA.6p then give 



U{t,to) = e^^*e-^(^+^)(*-*°)e-^^*«. (A.7) 



By differentiating U {t, to) with respect to t and substituting back the form of (1A.7P 
and flA.Sp . one finds a SE for the time-evolution operator: 

^^^^ = V{t)U{t,to)■ (A.8) 
Integrating both sides from to to t then gives 

U{t,to) = l-t [ V {t')U (t',to)dt'. (A.9) 
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By repeating the same process, but with respect to t^, we obtain 



to 



U{t,to) = l + i U (t, t') V it') dt'. 



(A.IO) 



With the help of Dyson's ordering operation |188j , the solutions to equations (1A.9P 
and ( lA.lOl) can be expressed as 



U (t, to) = 1 exp 
U {t, to) = ( exp 



-t [ V{t')dt' 

Jto 

I / V{t')&t' 



(A.ll) 
(A.12) 



where the + (— ) subscripts indicate that terms in the series expansion should be 
written from left to right in order of decreasing (increasing) times. The Heisenberg 
S'-matrix is often written as U (oo, — oo). Therefore, we must now show that for the 
scattering process yielding the state (lA.2p . it is meaningful to extrapolate equation 



(1A.7P to the limit of t, to = ±00. 

Let us transform (1A.2P with the help of (1A.3p . letting r ^ 00: 



which becomes, from (1A.7P 

mt)=e / e^^f/(t,T)0,dT. 
J —00 



)jdT, 



(A.13) 



(A.14) 



We take 



U (t, -00) = lim,^o+e / e'^U (t, T) dT, 

J —00 

POO 

U (00, t) = lim,^o+e / e-'^U (T, t) dT. 
Jo 



(A.15) 
(A.16) 



By thus defining our limits, it is possible to show that equations (1A.9P and (lA.lOp 
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do indeed reduce to 



U{t,-oo) = l-i f V {t')U {t\-oo)dt', (A.17) 

J —oo 

/CO 

U {t,t')V {t')dt'. (A. 18) 



We therefore conclude that 



/oo 
V {t')U {t',-oo)dt'. (A.19) 
-oo 



Now to justify the form of equation flA.7p as to ~^ — oo. We know the scattered 



states of the system in the asymptotic regime must satisfy the Lippmann-Schwinger 
equation (13 .34^ . such that 

'il:j = (Pj+G{Ej)Vijj. (A.20) 
where the ip are eigenstates of the total Hamiltonian, satisfying 

Hipj = Ejipj, (A.21) 

and G [Ej) is the retarded Green Function of the free system evaluated at energy 

Ef 

G (EA = . (A.22) 

From equations (1A.7P and (lA.lSp 



U (t, -oo) = e-^*e-(^+^)* lim,^o+ . ^ . , ^ ^ , \<I>,){U (A.23) 
with \(f)j){(f)j\ = 1. One can show that 

U {0,-oo)(f)j = ijj, (A.24) 

meaning U (0, — oo) takes an incident wave (pj, which is a stationary state of the 
interaction-free Hamiltonian, and produces an initial state ipj which is an eigenstate 
of the total Hamiltonian. Similarly, U (0, +oo) takes (pi to an eigenstate of the total 
Hamiltonian ipi corresponding to an outgoing wave (pi. These operators are found 
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to obey the unitarity condition 

U (0, -oo) U (-CX), 0) = 1, (A.25) 
f/(0,+oo)f/(+cx),0) = 1. (A.26) 

Then, using f EOSl) and f[04l) 

[/(t,-oo) =e^^* J]e-^^^*|^,)(0,|, (A.27) 

3 

giving 

\/(t)f/(t,-oo) = 5^|0.)(0.|e'(^-^^-)V|^,)(0,|. (A.28) 
«,i 

Finally, substituting into (1A.19I) : 

S = l-2mY, \<P^) {<Pj\^ {Ei - E,) (A.29) 

with 



% = (0, 1^1 0, + g^). (A.30) 

To first order in the scattering potential, the S'-matrix is then 

S^^^ = I - 2niV5 {Ei - Ej) , (A.31) 

where Ej and Ei are the initial and final eigenenergies of the free Hamiltonian, 
respectively. 
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The PF2 Source of UCN 



The PF2 source produces UCNs by further coohng neutrons that emerge from the 
Vertical Cold Source (VCS), currently operating at the Insitut Laue-Langevin in 
Grenoble. The structure of the VCS comprises a fuel element and liquid deuterium 
moderator immersed in a pool of heavy water, all encased in concrete shielding. 
As the neutrons interact with the moderator, the deuterium is heated to boiling 
point. The vapour is then recondensed in a heat exchanger, and flows back into 
the moderator vessel. The moderated neutrons are directed into a nickel guide, 
the lower half of which is filled with helium to increase cooling. They emerge in a 
turbine house with five exit ports, part of which is occupied by a device known as 
a Garching turbine. Half the beam feeds the turbine, which slows the neutrons to 
the meV regime by total reflection along semi-circular nickel blades. These neutrons 
emerge as ultra-cold beams, whose cross section is determined by the exit port to 
which they are distributed. The remaining part of the beam bypasses the turbine, 
and emerges through the upper of the five exit ports at a velocity of approximately 
40 ms^^. A schematic of the reactor and of the PF2 source is shown in figure iBTl 
below. 
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Figure B.l: A schematic of the Vertical Cold Source (left) and of the PF2 instru- 
ment at the ILL, Grenoble. The image of the VCS is reproduced from the web- 
site of the Technischen Universitat Miinchen, which operates the cold neutron source 
FRM II, whose design is extremely similar to that of the VCS. Image available at 
http: / / www.£rm2 . tum.de/en / technik/secondary-sources / cold-source /index. html 
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Appendix C 

The Energy Spectrum of a 
Heisenberg Ring 



The Heisenberg nearest-neighbour interaction is expressed by the Hamiltonian of 
equation fl7.ip : 

The eigenenergies of this Hamihonian in a periodic system are obtained by calcu- 
lating the expectation values {m\Hf)\m) between eigenstates of the form 

1 ^ 

where j represents the location of the spin flip. From equations (IC.ip and (lC.2p . 
one therefore obtains 

1 ^ 

= {m\Ho\m) = - ^ e''-^^-^'\j'\Ho\j) = 
T ^ 

T ^ 

= 'nY1 (^'''" + ^"*") = -2^cosfc^. (C.3) 
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This is the result used to compute the fidehties and timescales shown in figure I8.2[ 
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